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ABSTRACT 


Rosner,  John  Christopher.  Ph.D.,  Purdue  University, 
August  1969.  Theoretical  Study  of  Landing  Mat  Behavior. 

Major  Professor:  Dr.  Milton  E.  Harr. 

Mechanistic  models  are  developed  to  help  predict  the 
behavior  of  landing  mat  systems.  The  first  model,  which  is 
based  upon  elastic  theory,  is  capable  of  duplicating  the 
action  of  such  systems  under  static  loads.  The  associated 
assumptions  are:  1)  that  an  infinite  beam  is  the  structural 

equivalent  of  the  mat;  2)  that  the  subgrade  is  homogeneous; 

3)  that  horizontal  displacements  within  the  subgrade  are 
negligible;  and  4)  that  the  mat  always  remains  in  contact 
with  the  subgrade.  The  model  parameters  are  established  from 
simulations  of  full-scale  experimental  tests.  These  param- 
eters are  also  correlated  with  prototype  test  variables. 

Results  from  prototype  tests  indicate  that  the  model 
parameter  k,  the  subgrade  modulus,  decreases  as  trafficking 
of  the  section  progresses.  Contrary  to  common  belief  the 
model  behavior  is  found  to  be  extremely  sensitive  to  the 
magnitude  of  the  subgrade  modulus.  Results  also  indicate 
that  the  performance  of  dual-wheel  prototype  tests  can  be 
predicted  with  reasonable  confidence  by  the  procedure  devel- 
oped herein:  use  for  the  single-wheel  tests  appears  to  be 
somewhat  limited. 

xii 


'An  investigation  of  the  effect  of  end  joint  connections 
indicates  that  some  increase  in  service  life  can  be  attained 
by  strengthening  the  standard  end  joint  connectors. 

A second  mechanistic  model  is  developed  to  simulate 
actual  loading  sequences  and  to  provide  a means  of  esti- 
mating the  residual  deformations  of  the  mat  surface.  The 
applicability  of  this  model  is  demonstrated  for  both  single- 
wheel and  dual-wheel  tests.  Complete  evaluation  of  the 
potentialities  of  this  model  is  not  possible  from  the  avail- 
able data. 
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INTRODUCTION 


Present  indications  are  that  landing  mats  will  be  used 
extensively  in  the  forward  areas  of  any  future  air  opera- 
tions. This  is  particularly  true  in  theaters  of  operations 
where  time  and  the  availability  of  native  construction 
materials  will  not  permit  a more  permanent  type  of  installa- 
tion. Since  the  efforts  and  cost  required  to  produce  and 
deliver  landing  mats  to  theaters  of  operations  are  so  great, 
the  efficient  use  of  the  mat  is  essential.  Inefficient  use 
of  mats  may  prevent  the  timely  accomplishment  of  assigned 
missions . 

Existing  criteria  for  designing  runways  to  be  surfaced 
with  landing  mats  have  developed,  by  and  large,  from  full- 
scale  experimental  tests.  Prototype  tests  were  conducted 
by  the  Corps  of  Engineers  [10,  12r  on  individual  mat  types 
under  various  conditions  of  load,  number  of  coverages, 
wheel  spacing,  tire  pressure,  tire  ply  rating,  base  thick- 
ness, and  California  Bearing  Ratio  (CBR)  of  the  subbase 
material.  Most  of  these  tests  were  conducted  until  failure 
occurred.  Failure  was  judged  on  the  bas:'s  of  the  roughness 
of  the  mat  surface  or  excessive  breakage  in  the  mat.  In 
regards  to  the.  roughness  requirement,  the  test  section  was 
considered  to  have  failed  when  permanent  deviations  of  3 
inches  or  more  occurred  in  a 10  foot  length  in  any  direction 
within  the  traffic  lane.  When  breakage  developed  in  10  per- 
cent or  more  of  the  mat  panels  within  the  traffic  lane,  the 
test  section  was  considered  to  have  failed.  In  most  tests 
both  failure  criteria  were  met  simultaneously. 

The  results  of  each  test  when  failure  occurred  were 
used  in  connection  with  established  design  curves  developed 
originally  for  flexible  pavements.  Such  design  curves  are 
given  by  Yoder  [68]  and  the  procedure  for  using  the  obtained 
data  is  given  in  Corps  of  Engineers  publications  [8,  10]  . 

An  effective  thickness  of  flexible  pavement  was  assigned  to 
the  mat  using  the  design  curves.  The  performance  of  a land- 
ing mat  could  then  be  estimated  by  interpolation  ana  extra- 
polation. 

T The  numbers  in"brackets  refer  to  references  given  in  the 
List  of  References. 
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In  tests  where  failure  was  not  in  evidence  after  a 
maximum  established  number  of  coverages,  the  effective 
thickness  of  the  mat  was  determined  using  the  same  procedure 
given  above.  For  these,  the  effective  thickness  assigned  to 
the  mat  was  less  than  optimal  for  the  established  failure 
criteria. 

For  each  type  of  mat  tested  a set  of  design  curves  has 
been  developed.  In  line  with  this  procedure,  each  altera- 
tion in  the  structure  or  configuration  of  the  individual  mat 
elements  has  necessitated  a new  series  of  performance  tests. 

The  objective  of  this  study  was  to  develop  a reliable 
procedure  or  method  for  predicting  the  performance  of  land- 
ing mats.  It  was  anticipated  that  the  success  of  this  study 
would  serve  to  reduce  greatly  the  need  for  costly  perform- 
ance testing  and  the  procedure  could  be  used  for  comparisons 
of  the  relative  effectiveness  of  existing  mat  types.  In 
addition,  it  was  hoped  that  the  study  would  provide  a means 
whereby  more  efficient  utilization  of  landing  mats  could  be 
achieved . 

For  this  study  the  results  of  the  prototype  tests  as 
reported  by  the  Corp  of  Engineers  [10,  12]  were  used. 


REVIEW  OF  LITERATURE 


^ Existing  Theories  of  Surface  Displacements 

The  benefits  to  be  accrued  from  the  use  of  selected 
surface  materials  as  agents  for  the  transmission  of  wheel 
loads  to  the  existing  sub-soils  have  long  been  realized. 

The  behavior  of  these  materials,  be  they  concrete,  asphaltic 
mixtures,  chemically  stabilized  soils,  or  landing  mats,  has 
been  of  considerable  interest  to  the  engineer.  Many  theo- 
ries for  the  analysis  of  such  systems  have  been  proposed; 
yet,  no  general  theory  exists  today. 


Winkler  Concept 

One  of  the  initial  analyses  of  load  transfer  between 
surface  materials  and  a subgrade  was  provided  by  Winkler 
[67],  In  his  work  he  assumed  that  the  action  of  the  sub- 
grade was  analogous  to  a dense  fluid  or  a system  of  linear 
independent  springs  wherein  the  reactive  pressures  developed 
in  the  subgrade  were  proportional  to  the  deflection  of  the 
surface.  Winkler's  hypothesis  stated  that 


p » kw  (1) 

V 

where  p is  the  pressure  exerted  by  the  subgrade 
k is  a proportionality  constant 
w is  the  deflection  of  the  subgrade  surface. 

However  he  did  not  indicate  what  value  or  values  should  be 
associated  with  the  proportionality  constant,  k. 

In  his  analyses  of  railway  ties  and  rails,  Zimmermann 
[70]  utilized  the  Winkler  concept  and  concluded  that  the 
proportionality  constant  was  dependent  upon  the  type  of  sub- 
grade. Experimental  studies  conducted  by  Foppl  [ 19]  indi- 
cated that  the  response  of  such  systems  could  be  approxi- 
mated by  using  the  Winkler  hypothesis. 
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In  a mathematical  treatment  of  the  subject,  Hayashi 
[ 25]  employed  the  tinkler  hypothesis  to  analyse  beams  rest- 
ing on  soil.  He  presented  the  solution  for  the  condition 
where  shearing  forces  existed  at  the  interface  of  the  beam 
and  the  soil.  The  governing  differential  ecuation  took  the 
form 


A 

El 


ox 


c(x) 

El 


(2) 


wherein  El  is  the  rigidity  of  the  beam 
A is  a constant 
q(x)  is  the  loading  function. 


If  interfacial  shear  can  be  neglected,  the  second  term  in 
Eqn.  (2)  may  be  omitted.  From  a later  investigation  which 
considered  the  presence  of  interfacial  shear  forces,  Florin 
[18  ] concluded  that  these  frictional  forces  had  negligible 
influence  on  the  distribution  of  the  vertical  pressure 
exerted  by  the  soil. 

In  1926,  Westergaard  [64]  used  the  Winkler  hypothesis 
to  compute  the  stresses  in  concrete  highway  pavements. 

Later  [65]  he  extended  his  analysis  to  airport  pavements. 

In  his  initial  study,  he  observed  that  an  increase  of  the 
subgrade  modulus  from  50  lbs./cu.  in.  to  200  lbs./cu.  in. 
produced  only  minor  variations  in  pavement  stresses.  There- 
fore it  was  surmised  that  an  approximate  single  value  of  the 
subgrade  modulus,  k,  should  be  sufficient  for  determination 
of  the  stresses  in  pavements.  It  was  also  suggested  that  a 
standard  procedure  be  developed  whereby  the  value  of  k could 
be  determined . 


Prior  use  of  the  Winkler  hypothesis  assumed  that  the 
subgrade  followed  the  structural  member.  This  necessitated 
the  development  of  normal  tensile  stresses  at  the  interface 
of  the  load  transfer  mechanism  and  the  subgrade.  Murphy 
[39]  investigated  the  stresses  and  deflections  of  plates 
whose  edges  were  free  of  foundation  support.  In  addition 
to  observing  that  the  size  of  the  plate  had  a considerable 
effect  upon  the  developed  stresses,  calculations  indicated 
that  an  increase  in  the  stiffness  of  the  foundation  reduced 
appreciably  the  stresses  and  deflections  of  the  place.  In 
an  analogous  manner,  Harr  and  Leonoras  [24  ],  employing  the 
Winkler  hypothesis,  demonstrated  that  appreciable  stresses 
and  deflections  coulu  be  produced  in  concrete  slabs  when 
moisture  ano/or  temperature  gradients  across  the  depth  were 
present. 
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In  1946,  Hctenyi  [26]  presented  a collection  of  solu- 
tions for  beams  on  './inkier  foundations;  finite  as  well  as 
infinite  length  beams.  To  remove  the  discontinuity  in  the 
deflection  pattern  that  occurs  at  the  end  of  a finite  length 
beam  using  the -Winkler  hypothesis,  ho  suggested  that  the 
finite  beam  be  placed  over  an  infinite  beam  which  in  turn 
wan  embocuoci  in  the  system  of  springs.  Hotenyi  did  not  give 
numerical  values  for  the  subgradc  modulus,  nor  did  he  sug- 
gest what  factors  rnay  influcnco  this  measure. 

In  a study  conducted  in  1955,  Drapkin  [16]  obtained  the 
solution  for  finite  length  beams  utilizing  the  principle  of 
superposition.  Ho  noted  that  contrary  to  prevailing  prac- 
tice increasing  the  length  of  a beam  did  not  reduce  sub- 
stantially the  maximum  vertical  foundation  pressures. 

In  1955,  Terzaghi  [62  ] presented  a critical  review  of 
the  history  and  development  of  existing  theory  based  upon 
the  subgrade  modulus.  He  pointed  out  that  the  magnitude  of 
the  subgrade  modulus  was  dependent  upon  the  dimensions  of 
the  loaded  area  as  well  as  the  elastic  properties  of  the 
soil.  He  also  established  a procedure  whereby  consistent 
values  for  the  subgrade  modulus  could  be  obtained.  Terzaghi 
noted  that  intelligent  use  of  the  Winkler  hypothesis  would 
produce  reliable  stresses  and  bending  moments  but  the  theory 
failed  to  provide  reasonable  estimates  for  settlements. 

In  1962,  Lenczer  [36]  experimentally  investigated  the 
effect  of  soil  depth  on  the  magnitude  of  the  subgrade 
modulus.  He  observed  that  for  shallow  depths,  less  than  12*. 
inches,  the  variation  in  the  numerical  values  of  k was 
appreciable.  He,  as  did  Terzaghi  before  him,  found  that  the 
subgrade  modulus  was  dependent  upon  the  size  of  loaded  area. 
In  addition,  he  presented  an  empirical  relationship  wherein 
the  subgrade  modulus  could  be  taken  as  a function  of  deposit 
depth . 

Many  other  investigators  [21,  29,  38,  51]  have  employed 
the  Winkler  hypothesis  or  minor  variations  thereof  to  the: 
analysis  of  load  transfer  systems. 


Elastic  Solid  Models 

In  the  1930 's  with  the  advance  of  soil  mechanics,  ques- 
tions were  directed  concerning  the  validity  of  the  Winkler 
hypothesis.  In  1937,  Biot  [4]  presented  the  means  of  com- 
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puting  the  contact  pressures  on  the  base  of  an  infinite 
elastic  beam  resting  on  a semi-infinite  elastic  solid.  This 
solution  made  it  possible  to  determine  an  equivalent  value 
for  the  subgrade  modulus,  k.  Subsequent  investigations  led 
to  the  conclusion  that  the  subgrade  modulus  was  a compli- 
cated function  dependent  not  only  on  the  elastic  modulus  of 
the  soil  and  the  width  of  the  beam  but  also  upon  the  beam's 
flexural  rigidity.  Biot  concluded  that  no  unique  value  of 
the  subgrade  modulus  could  be  assigned  to  a given  subgrade. 

In  1943,  Burmister  [ 6 ] obtained  the  solution  for  a two 
layered  system  acted  upon  by  a uniform  circular  load.  Each 
layer  was  assumed  to  be  a homogeneous,  isotropic,  and 
linearly  elastic  solid.  A perfectly  rough  interface  condi- 
tion was  assumed  to  exist  between  layers  with  the  lower 
layer  providing  continuous  support  to  the  upper  layer. 

Lemcoe  [35  ] confirmed  the  results  given  by  Burmister  and 
also  investigated  the  condition  of  a frictionless  interface. 
He  reported  that  only  a very  slight  change  occurred  in  the 
vertical  normal  stresses  as  the  interface  condition  changed 
from  perfectly  rough  to  perfectly  smooth. 

The  semi- infinite  elastic  solid  model  and  an  elastic 
layer  of  finite  thickness  were  employed  by  Pickett  [43]  in 
his  investigation  of  the  behavior  of  concrete  pavements. 
Influence  charts  for  deflections  and  moments  at  different 
points  in  the  concrete  slab  were  presented  for  both  an  elas- 
tic solid  subgrade  and  a dense  liquid  subgrade.  However, 
he  gave  no  suggestions  for  associating  the  model  parameters 
with  given  subgrade  conditions.  The  solution  for  an  elastic 
solid  of  finite  depth  was  also  presented  by  Sovinc  60  . As 
an  aid  to  obtaining  this  solution  he  approximated  the  dis- 
placements within  the  solid  by  a double  Fourier  Series. 

In  1961,  Davis  and  Taylor  [13  ] investigated  the  influ- 
ince  of  horizontal  loading  on  the  surface  displacements  of 
an  elastic  layer.  They  concluded  that  care  in  selecting 
the  values  of  the  elastic  constants  was  necessary  for  this 
type  of  loading  because  some  portions  of  the  soil  mass  ex- 
periences a stress  reduction.  In  1963,  Lee  [ 34  ] concluded 
that  for  horizontal  surface  loads  the  normal  stress 
distribution  along  the  contact  surface  of  a flexible  strip 
was  essentially  independent  of  the  magnitude  of  the  contact 
shearing  stresses.  Lee  also  expressed  the  belief  that  the 
assumption  of  a smooth  contact  surface  leads  to  conservative 
estimates  of  the  bending  moment  and  shearing  force  induced 
in  the  strip,  in  an  earlier  study,  Leonov [ 37]  also  found 
that  horizontal  forces  have  practically  no  influence  on  the 
distribution  of  the  vertical  reactive  pressure. 


6 


In  an  attempt  to  simplify  future  computations#  DeBeer 
[14]  using  an> iterative  procedure#  obtained  the  contact 
pressure  distribution  beneath  a beam.  He  suggested  that  a 
second  or  fourth  order  polynominal  could  be  used  to  approxi- 
mate the  contact  pressure  distribution.  Using  this  relation- 
ship he  obtained  reasonably  correct  deflections  and  bending 
moments.  From  the  results  of  an  experimental  investigation 
DeBeer  [ 15  ] concluded  that  the  second  degree  parabolic  dis- 
tribution reflected  the  results  of  the  tests  regardless  of 
the  beam  stiffness.  Comparisons  of  the  test  results  with 
computations  based  upon  the  Winkler  hypothesis  were  also 
made.  Barden  [ 2 ] in  his  studies  of  finite  beams  employed 
DeBeer 's  approximation  for  the  contact  pressure  and  con- 
cluded that  the  results  obtained  were  more  valid  than  those 
obtained  from  the  Winkler  analysis. 

Many  procedures  and  techniques  have  been  developed 
[ 57,  41,  23,  34,  48  ] to  approximate  the  distribution  of 
surface  contact  pressures.  A method  developed  by  Zhemochkin 
[69  ] assumed  that  the  load  transfer  element  was  connected 
to  the  elastic  solid  by  linkages  (two-force  members).  By 
treating  the  forces  in  these  linkages  as  unknowns  and  using 
the  "method  of  virtual  displacements"#  he  was  able  to 
approximate  the  contact  pressures  transmitted  by  the  elastic 
solid.  A high  degree  of  accuracy  could  be  achieved  by  using 
many  linkages  in  the  system.  This  procedure  gave  results 
which  agreed  well  with  the  more  rigorous  mathematical  solu- 
tion. Sinitsyn  [ 59  ] extended  the  application  of  Zhemochkin' s 
procedure  and  Barden  [ 1] later  presented  another  variation. 

In  his  discussion  to  Barden's  paper,  Chueng  [7]  indicated  a 
further  simplification  and  noted  that  this  method  had  been 
used  in  Asiatic  countries  in  design  for  many  years. 

A recent  theory  for  analyzing  structures  on  elastic 
foundations  based  upon  a general  variational  method  was, 
developed  by  Vlasov  and  Leont ' ev  [ 63 ] . This  theory  is  more 
complex  than  the  Winkler  theory  yet  simpler  than  the  elas- 
tic half-space  analysis.  The  form  of  the  basic  differential 
equation  given  to  describe  the  state  of  strain  in  a loaded 
single-layer  foundation  had  the  same  form  as  that  obtained 
earlier  by  Filonenko-Borodich  [ 17  ] and  by  Pasternak  [42]. 

The  method  possesses  great  flexibility.  By  judicious  selec- 
tion of  the" coefficients  in  the  differential  equation,  the 
equation  for  either  the  Winkler  hypothesis  or  the. elastic 
solid  model  may  be  obtained.  In  addition,  boundary  condi- 
tions are  easily  introduced  and  the  developed  procedure  can 
be  readily  extended  for  complex  three-dimensional  considera- 
tions. This  method  forms  the  background  of  a large  part  of 
the  analysis  developed  in  this  report. 
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Viscoelastic  Models 


Recently  attention  has  been  directed  toward  application 
of  viscoelastic  theory  to  the  response  of  soils  subjected 
to  static  and  cynawic  loading.  Many  rheological  models  have 
been  proposed  for  simulating  the  stress-strain-time  behavior 
of  soils  and  for  modeling  studies  of  creep  phenomena  and 
stress  relaxation  in  soil. 

Freudenthal  and  Lorsch  [20]  were  among  the  first  inves- 
tigators who  employed  a linear  viscoelastic  analysis  to 
infinite  beams.  In  their  study  the  soil  support  was  replaced 
by  a series  of  Kelvin,  Maxwell,  or  Standard  Linear  Solid 
elements.  For  each  of  these  elements  the  authors  were  able 
to  develop  relationships  for  the  deflection  of  an  infinite 
beam  subjected  to  time  invariant  concentrated  and  uniform 
loads. 

In  a brief  note  in  1958,  Reissner  [52]  presented  the 
solution  of  a thin  plate  resting  on  a viscoelastic  founda- 
tion which  possessed  shear  interaction;  a frictionless 
plate-foundation  interface  was  assumed.  Two  years  later, 
Pister  and  Williams  [ 47]  extended  the  work  of  Reissner  to 
incorporate  a rough  interface  condition.  In  addition,  they 
presented  response  curves  for  the  maximum  deflection  and 
moment  for  an  applied  step  force.  Later  Pister  [46] 
presented  the  solution  for  the  axisymmetric  bending  of  a 
viscoelastic  plate  of  finite  thickness  supported  by  a visco- 
elastic foundation  of  infinite  extent.  He,  as  did  Reissner, 
assumed  a frictionless  interface  condition. 

In  his  study  of  the  creep  behavior  of  snow  foundations, 
Kerr  [ 28]  employed  a linear  viscoelastic  media.  The  time 
invariant  load  was  applied  directly  to  the  foundation  rather 
than  by  a load  transfer  mechanism  (such  as  a beam  or  a plate) . 

Hoskin  and  Lee  [ 27]  used  linear  viscoelastic  models  for 
their  analysis  of  the  stress  and  deformation  characteristics 
of  a subgrade,  strengthened  by  a flexible  surface  plate  and"" 
subjected  to  a sudde  ily  applied  uniformly  distributed  invari- 
ant load.  They  concluded  that  the  Maxwell  model  gave  unreal- 
istic values  for  deflections,  subgrade  pressures,  and  plate 
bending  moments  for  long  time  loading;  however,  the  use  of 
the  Standaro  Linear  Solici  mooel  was  found  to  be  satisfactory. 
The^c  also  noted  that  under  certain  circumstances  (by  use  of 
transform  theory  the  time  dependence  of  the  system  could  be 
removed)  stress  and  deflection  analyses  of  a system  contain- 
i°9  linear  viscoelastic  components  could  be  treated  as  an 
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analogous  clastic  problem  having  tho  enmo  geometric  boundary 
conuitionc.  This  analogy  is  known  as  the  "correspondence 

principle".  '"l - - 

viscoc.  Ins  lie 


lhcy  concluded  that  the  response  of  a linear 
tlorial  to  stress  could  be  prooictcd  by  an 


elastic  .y>at o eve  1 wherein  the  elastic  constants  wore  replaced 
by  time-dependent  parameters. 


.ennpery  [5J(  54]  reported  that  the  solution  of  many 
non-linear  viscoelastic  problems  may  ' 


[5  3,  54] 

also  be  closely  approx- 
imated by  reducing  them  to  an  equivalent  elastic  condition 
wherein  the  parameters,  which  characterize  the  non-linear 
viscoelastic  material,  are  taken  as  functions  of  time.  The 
validity  of  this  procedure  (quasi-elastic  method  of  analysis) 
was  demonstrated  by  Schapery  [53]  for  a viscoelastic  canti- 
lever-beam subjected  to  a concentrated  load  at  the  free  end. 
Schapery  [54]  also  applied  this  procedure  to  the  analysis  of 
the  creep  of  glass  fiber  reinforced  resin  under  an  uniaxial 
stress  condition  and  demonstrated  that  the  procedure  could 
also  predict  the  creep  response  of  Polystyrene  under  a vary- 
ing stress  condition.  Three  different  samples  of  Polysty- 
rene were  subjected  to  a triple-step  stress  application. 
Predictions  of  the  magnitude  of  the  strains  were  within  9% 
of  those  observed  from  experimental  investigation  over  an 
elapsed  time  of  1512  hours. 


Konder  and  Krizek  [31]  investigated  the  creep  response 
of  a commercially  available  cohesive  soil,  Jordan  Buff, 
under  uniaxial  and  constant  compressive  stress.  They  were 
able  to  predict  creep  response  for  this  material  which 
exhibited  definite  non-linear  behavior  with  an  expression 
of  the  form 

2 

e = C(o)  + Cj  (o) (log  t - Cfc  log  t)  (3) 


where  C (a),  C*  (a)  are  constants  with  respect  to  time 

but  are  stress  dependent 
C2  is  a constant 
t is  the  time  variable. 

Singh  and  Mitchell  [58]  studied  the  creep  response  of 
a number  of  soil  types  subjected  to  various  test  conditions. 
They  developed  an  expression  for  the  strain  of  the  form 


c = C(o)  + C,  eoC('’)  (t)''m  (in  / 1)  (4) 


9 


where 


C(o)  is  constant  which  is  dependent  on  the 
strain  at  unit  time 
Cj  is  a constant 
a is  a dashpot  constant 
D (o  ) is  the  dcviator  stress 

M is  the  slope  of  logarithm  of  strain  rate  versus 
logarithm  of  time  (straight  line) 
t is  time. 

They  concluded  that  the  expression  was  applicable  irrespec- 
tive of  whether  the  soils  were  undisturbed  or  remolded,  wet 
or  dry,  normally  consolidated  or  over-consolidated,  or 
tested  in  a drained  or  undrained  condition.  They  stated 
* that  the  developed  creep  function  reflected  the  effect  of 
soil  composition,  soil  structure,  stress  history,  stress 
intensity  and  the  slope  of  the  strain  rate  vs.  time  relation- 
ship on  a log-log  plot.  They  indicated  that  further  research 
was  being  conducted  to  attest  the  validity  of  their  creep 
function  for  repetitive  load  conditions. 

In  his  analysis  of  viscoelastic  layered  pavement  sys- 
tems, Barksdale  [3]  was  able  to  develop  a creep  compliance 
response  for  an  asphalt  mixture  and  one  for  a clay  subgrade. 
Both  relationships  were  developed  for  materials  subjected  to 
repetitive  stationary  step  loadings. 

The  primary  limitation  of  applying  viscoelastic  models 
to  real  soil  behavior  has  been  the  difficulty  in  assigning 
representative  numerical  values  to  the  parameters  for  even 
the  simplest  model. 


Prototype  Tests  on  Landing  Mats 

Accelerated  traffic  tests  simulating  aircraft  taxiing 
operations  were  condi  :ted  by  the  Corps  of  Engineers  [ 10]  on 
test  sections  constructed  with  a range  of  subgrade  strengths 
that  represented  airstrips  surfaced  with  landing  mats.  The 
purpose  of  th  'se  tests  was  to  provide  data  on  the  service 
life  of  mat-surfaced  airstrips  under  various  conditions  of 
wheel  load,  tire  pressure,  and  subgrade  strength.  As  noted 
previously,  design  curves  were  developed  by  modification  of 
the  conventional  CBR  (California  Bearing  Ratio)  design 
curves  for  flexible  pavements. 
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The  traffic  simulated  the  operations  of  military  air- 
c'ni.t  with  both  single-  and  dual-wheel  assemblages.  Single- 
v.Yioal  loads  ranged  from  10,000  to  50,000  lbs  and  tire  pros- 
surer,  ranged  from  40  to  300  psi.  For  dual— wheel  assembly, 
Iona a of  50,000  lbs  were  used  with  a wheel  spacing  of  37.5 
in.  (center  to  center)  and  tire  pressures  ranged  from  100  to 
300  psi.  These  loading  arrangements  were  applied  by  a load 
cart  C 1 0 3 . bach  test  section  was  approximately  26  ft.  wide 
■and  30  ft.  long  with  the  traffic  applied  over  a width  of 
12  ft.-- 


Two  types  of  subgrades  were  used:  (1)  a fat  clay  (CH) 

with  an  average  liquid  limit  of  60  and  a plasticity  index 
of  39  and  (2)  a lean  clay  (CL)  with  an  average  liquid  limit 
of  36  and  a plasticity  index  of  15.  The  strength  of  the 
subgrade  was  classified  into  three  groups:  (1)  high  strength 

with  a CBR  above  20,  (2)  medium  strength  with  a CBR  between 

7 and  20,  ahd  (3)  low  strength  with  a CBR  between  3 and  7. 

The  tests  were  made  on  standard  M6,  M8,  and  M9  landing 
mats.  The  M6  mat  is  a pierced  steel  p,JLank  with  a moment  of 
inertia  per  foot  of  width  of  0.069  in.  ; the  M8  mat  is  a 
heavy,  deep-ribbed  steel  mat  with  a moment  of  inertia  per 
foot  of  width  of  0.269  in.  The  M9  mat  is  a deep-ribbed 
aluminum  mat  with  an  average  moment  of  inertia  per  foot  of 
width  of  0.618  in.  4 These  mats,  which  are  rectangular  in 
shape,  were  placed  on  the  cutegrado  in  a masonry  typo  arrange- 
ment and  connected  by  integral  locking  lugs  on  the  sides  arid 
hooked  connectors  on  the  ends. 

The  tests  were  continued  until  the  test  section  failed 
or  until  700  coverages  had  been  completed  by  the  load  cart. 
Failure  of  the  sections  was  judged  on  the  basis  of:  (1) 

development  of  roughness  of  the  mat  surface  to  the  point  of 
endangering'  aircraft  operations  and  (2)  excessive  mat  break- 
age . 

In  1967,  the  Corps  of  Engineers  conducted  an  extensive 
study  [12  ] to  develop  criteria  for  the  efficient  design  of 
aircraft  landing  gear  for  aircraft  required  to  operate  from 
mat-surfaced  airfields.  Traffic  tests  were  conducted  with 
numerous  combinations  of  wheel  configurations,  loads,  and 
tire  pressures.  The  wheel  configurations  varied  from  a 
single-wheel  up  to  12  wheels:  the  loading  varied  from  35,000 

to  273,000  lbs  and  the  tire  pressures  ranged  from  50  to  250 
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psi.  The  tv/o  soils,  which  formed  the  subgrade,  exhibited 
only  minor  differences  in  characteristics.  One  soil  was  a 
fat  clay  (CII)  with  a liquid  limit  of  50  and  a plasticity 
index  of  21  while  the  other  soil  was  a fat  clay  (CH)  with  a 
liquic.  limit  of  61  and  a plasticity  index  of  37.  The  in- 
place,  initial  strength  of  these  subgrades  as  indicated  by 
CBR  values  ranged  from  1.1  to  9.0. 

These  te_ts  were  made  on  the  M8  mat  (previously 
described)  and  the  modified  Til  mat  which  is  a lightweight, 
extruded-aluminum  panel  with  an  abrasive  surface  [ 11] . The 
moment  of  inertia  per  foot  of  width  of  the  Til  mat  is 
1.368  in.  **  The  method  of  placv’.ment  of  the  rectangular- 
shaped Til  mat  was  also  of  the  masonry  design.  The  behavior 
al  characteristics  and  performance  of  the  mat  surfaces, 
whether  loaded  or  unloaded,  are  well  documented  [ 12  ] for 
each  test  at  various  coverage  levels.  The  single-wheel  and 
dual-wheel  data  from  this  test  series  constitute  the  basis 
of  the  investigation  reported  herein. 


Application  of  Theories  to  Landing  Mats 

It  is  well  recognized  [9]  that  the  action  of  landing 
mats  in  distributing  wheel  loads  to  the  subgrade  is  compli- 
cated and  that  the  exact  mechanism  through  which  distribu- 
tion is  accomplished  has  not  been  defined.  The  consensus 
of  opinion  is,  however,  that  the  mat  distributes  loads  in  a 
manner  similar  to  that  of  a beam  or  flat  plate.  Consequent- 
ly, most  theoretical  studies  of  landing  mats  have  employed 
this  type  of  analysis. 

In  1951,  Pickett  [44]  conducted  a theoretical  investi- 
gation of  the  behavior  of  landing  mats.  In  this  investiga- 
tion, he  considered  the  landing  mat  capable  of:  (1)  tensile 

strength  only,  (2)  flexural  strength  only,  and  (3)  both 
flexural  and  tensile  strength.  Initially,  the  assumption 
was  made  that  the  composite  mat  and  subgrade  material  could 
be  represented  by  a model  composed  of  a thin  membrane  of 
infinite  extent  supported  on  a liquid  subgrade  incapable  of 
supporting  shear  stresses.  Vertical  loads  were  applied 
directly  to  the  membrane.  The  governing  differential  equa- 
tion for  this  model  was 


2 

T A w - kw  - -q  (5) 


12 


whege  T is  tho  membrane  tension  (assumed  to  be  constant) 
a is  the  Laplacian  operator 
w is  the  deflection  of  the  surface 
k is  the  density  of  the  liquid 
q is  the  intensity  of  pressure  in  tho  direction 
of  positive  deflection. 

In  earlier  studies,  Schiel  [56  ] and  Filonenko-Borodich  (17) 
arrived  at  this  identical  differential  equation  in  their 
analyses  of  beams  on  elastic  foundations.  After  obtaining 
the  solution  of  Eqn.  (5),  Pickett  developed  an  influence 
chart  by  which  the  deflection  of  the  surface  and  the  reac- 
tive pressure  could  be  determined  graphically  for  a uniform 
load  of  any  configuration. 

Next  Pickett  [44]  assumed  that  the  composite  mat  and 
subgrade  material  could  be  represented  by  a membrane  sup- 
ported  on  an  elastic  solid.  For  this  model,  influence 
charts  were  prepared  for  both  deflection  and  reactive  pres- 
sures for  two  finite  thicknesses  of  the  elastic  solid  in 
addition  to  the  case  of  infinite  thickness.  As  could  be 
expected,  the  presence  of  the  membrane  reduced  the  magnitude 
of  the  deflection  as  compared  to  the  simple  elastic  half- 
space model.  A model  composed  of  a thin  plate  type  of  mem- 
brane (wherein  the  membrane  tension  could  vary  with  direc- 
tion at  a point  and  vary  from  point  to  point)  and  a subgrade, 
assumed  to  be  a dense  liquid  with  respect  to  vertical  reac- 
tions but  providing  elastic  resistance  to  horizontal  dis- 
placements of  the  membrane,  was  investigated  by  the  finite 
difference  method.  Mathematical  difficulties  prevented  the 
completion  of  this  solution. 

In  a. later  study,  Pickett  [45]  extended  his  analytical 
work  to  include  considerations  of  orthotropic  mats  on  elastic 
subgrades.  In  these  investigations,  he  considered  the  mat 
capable  of  only  flexural  resistance  and  assigned  to  the  mat 
a wide  range  of  combinations  of  transverse,  longitudinal, 
and  torsional  stiffnesses.  He  concluded  that  a given  amount 
of  transverse  flexural  rigidity  was  several  times  as  effec- 
tive in  reducing  maximum  deflection  as  an  equal  amount  of 
torsional  rigidity  and  that  both  of  these  rigidities  were 
of  less  relative  importance  than  the  longitudinal  rigidity. 

It  was  further  noted  for  widely  distributed  loads  that  all 
mat  rigidities  were  of  less  relative  importance  than  the 
subgrade  rigidity. 


In  all  of  Pickett's  studies  it  was  assumed  that  the  mat 
and  subgrade  remained  in  contact  at  all  points  and  that  the 
shear  Developed  at  the  interface  could  be  neglected.  In 
none  of  the  studies  aid  he  actually  utilize  his  models  to 
predict  the  performance  of  landing  mats. 

In  1955,  the  Corps  of  Engineers  [9]  presented  the 
results  of  full-scale  landing  mat  tests  and  model  tests 
which  were  compared  with  the  analytical  studies  conducted 
by  Pickett  [44,  45  ].  It  was  concluded  that  the  beneficial 
effects  of  the  landing  mat  was  dependent  upon  its  longitu- 
dinal, transverse,  and  to  a very  minor  extent  its  torsional 
rigidity.  It  wrs  further  noted  that  because  of  slack  in 
the  end  joint  connections  and  practical  considerations  in 
laying  the  mat  little,  if  any,  over-all  membrane  action  was 
possible.  The  report  concluded  that  landing  mats  act  pri- 
marily "in  flexure"  and  that  additional  stiffness  at  the 
end  joints  add  considerably  to  the  smoothness  of  the  operat- 
ing surface  under  traffic  conditions,  approximately  doubling 
the  service  life  of  the  mat. 


MODEL  DEVELOPMENT  AND  SIMULATION 


Introduction 


A mechanistic  model  was  sought  whose  deflection  behav- 
ior reflected  the  action  of  landing  mats.  Initial  consider- 
ation was  given  to  the  selection  of  a representative  load 
transfer  element.  In  accordance  with  the  findings  of  the  i 
Corps  of  Engineers  [9],  a membrane  was  eliminated  as  a pos- 
sible element.  Also  as  indicated  in  Photograph  5 of  this 
report,  (particularly  for  a group  of  M8  type  mats)  there  is 
an  apparent  lack  of  resistance  to  lateral  bending  due  to  the 
presence  of  longitudinal  joints.  This  would  tend  to  exclude 
a thin  plate  from  consideration.  Thus  for  the  sake  of  expe- 
diency, a beam  was  selected  for  the  basic  load  transfer 
element.  .......  , 

The  second  consideration  was  the  idealization  of  the 
soil  media.  The  Winkler  hypothesis,  Eqn.  (1),  was  not  con- 
sidered suitable  due  to  its  shortcomings  as  pointed  out  by 
Terzaghi  [ 62  ] and  its  neglection  of  interfacial  shear 
stresses.  Initially,  a viscoelastic  model  was  not  used 
because  of  its  complexity  and  the  anticipated  difficulty  of 
correlating  the  model  parameters  with  the  given  soil  proper- 
ties. Since  it  was  shown  [12]  that  the  "average  deflection" 
increased  with  the  number  of  coverages,  a conventional 
elastic  solid  model  was  also  not  directly  applicable.  How- 
ever, this  behavioral  characteristic  could  be  accommodated 
by  an  elastic  solid  model  wherein  the  model  parameters  are 
made  coverage  dependent.  Such  a quasi-elastic  model  for 
the  soil  media  was  employed  initially  in  this  study. 


Development  of  Mat-Soil  Mod e 1 

Since  it  was  anticipated  that  the  mat-soil  model  param- 
eters would  have  to  be  established  by  analytical  simulation 
of  prototype  test  data,  careful  consideration  was  given  to 
the  method  of  analysis.  The  general  variational  method  of 
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analysis,  developed  by  Vlasov  and  Leont’ev  [63]  for  elastic 
foundations,  appeared  to  provide  the  desired  degree  of  flex- 
ibility. A detailed  presentation  of  this  method  for  a 
single  layer  foundation  is  given  in  Appendix  A. 

The  reaction  of  a beam  of  infinite  length  to  loading 
was  selected  as  the  mechanistic  equivalent  to  that  of  the 
mat  ur.d_r  similar  circumstances.  This  selection  was  predi- 
cated upon  several  prevailing  conditions.  First,  actual 
field  operations  a au  test  procedures  employed  by  the  Water- 
ways Experiment  Station  112  ] demonstrated  that  the  mat  ele- 
ments extend  laterally  for  an  appreciable  distance  outside 
the  normal  traffic  lane.  The  traffic  lane  widths  for  the 
single  wheel  tests  ranged  from  4.75  ft.  to  12.00  ft.  A 
minimum  extension  of  10.16  ft.  beyond  the  traffic  lane 
existed  for  all  cases.  Typically  the  width  over  which  the 
wheel  load  was  applied  was  only  14.7  inches.  Secondly,  the 
longitudinal  joint  used  for  the  mats  (necessitated  by  the 
construction  procedure)  provided  virtually  no  moment  trans- 
fer from  one  mat  element  to  the  next.  However,  as  attested 
by  Photograph  4 of  reference  9,  in  the  transverse  direction 
an  appreciable  amount  of  moment  transfer  was  afforded  by 
the  end  joint  connectors.  Tests  [ 9 ] on  the  M8  mat  indi- 
cated that  the  longitudinal  rigidity  was  approximately  150 
times  larger  than  the  lateral  rigidity. 

The  width  of  the  infinite  beam  was  taken  as  the  length 
of  a rectangle  whose  area  was  equivalent  to  the  tire  print 
area  and  whose  width  was  equal  to  the  maximum  width  of  the 
tire  print.  In  those  cases  where  this  equivalent  beam  width 
was  greater  than  the  width  of  the  actual  mat  element  the 
longitudinal  joint  was  neglected.  Initially,  it  was  assumed 
that  complete  continuity  existed  at  the  end  joints  of  the 
mat  elements;  that  is,  the  end  joint  connections  provided 
total  shear  and  moment  transfer  between  mat  elements. 

In  addition  to  the  above,  the  following  assumptions 
were  made: 


1. 

The  wheel  loads  can  be  represented 
distributed  loads. 

by  uniformly 

2. 

The  beam  obeys  Eulerian  conditions 
the  stress  level. 

regardless  of 

3. 

The  beam  and  the  soil  always  remain 

in  contact. 

4. 

Horizontal  displacements  within  the 
are  negligible. 

soil  media 
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Under  these  assumptions,  the  response  of  the  mat-soil  model 
is  given  by  the  differential  equation  (Ecn.  (B5)  , /.opendi;: 
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At  the  present  time  there  is  no  reliable  information 
available  as  to  the  actual  distribution  of  displacements 
with  depth  in  a soil  layer.  This  is  particularly  true  for 
the  selected  soil  model  wherein  horizontal  displacements  in 
the  soil  media  are  assumed  to  be  of  negligible  magnitude. 
Information  available  for  homogeneous  deposits  indicates 
that  the  distribution  of  displacements  with  depth  is  non- 
periodic and  has  a maximum  at  the  surface.  Tests  conducted 
by  the  Corps  of  Engineers  [9]  on  a rubber  subgrade  loaded 
both  without  a mat  and  through  a steel  mat  indicate  an 
asymptotic  attenuating  distribution  of  displacement  with 
depth-  (see  Plate  17  of  reference  9) . It  appears  reasonable 
to  assume  that  with  different  mat  rigidities  the  same 
general  shape  would  be  maintained  but  that  the  rate  of 
attenuation  would  vary  with  depth.  Functional  representa- 
tions of  this  type  of  distribution  may  take  many  forms. 

One  form,  suggested  by  Vlasov  and  Leont’ev  [63] , assumed  a 
ratio  of  hyperbolic  functions  as 

sinh  y 

IT-  (10) 

sinh  Y| 
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where  H is  the  thickness  of  the  soil  layer 

y is  the  distance  from  the  subgrade  surface 
Y is  the  dimensionless  parameter 


With  the  distribution  of  Displacements  described  by  Eon. 
(10)  and  after  applying  the  assumption  of  negligible  hori- 
zontal displacements,  the  stresses  can  be  expressed  as 
(from  Eqns.  (;.  16)  , Appendix  A) 
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A typical  distribution  of  the  above  stresses  is  given  in 
Figure  1. 


With  the  form  of  (y)  taken  as  in  Eqn.  (10),  any  vari- 
ation in  the  rate  of  attenuation  of  the  displacement  due  to 
different  mat  rigidities  can  be  incorporated  by  judicious 
selection  of  the  parameter y . As  indicated  in  Figure  2, 
fj  (y)  does  not  produce  a unique  distribution  until  the 
dimensionless  parameter  y is  established.  Since  the  value 
of  y could  not  be  established  from  any  previous  studies,  it 
was  necessary  to  examine  the  "average  deflection"  patterns 
for  the  cases  at  hand  with  the  hope  that  a simulation  pro- 
cedure could  be  developed  which  would  yield  reasonable 
measures  of  this  parameter. 


In  aodition  to  the  parameter  y , some  measure  was 
required  for  the  equivalent  layer  thickness,  K.  In  most  of 
the  test  sections,  the  natural  soil  was  excavated  to  a 
depth  of  6 ft.  and  the  excavation  was  backfilled  under  con- 
trolled conditions  to  produce  tne  required  soil  strengths. 

In  some  test  sections,  the  subgrade  was  controlled  only  for 
a depth  of  two  feet.  An  analytical  study  of  the  effect  of 
layer  thickness  on  the  deflection  pattern  was  made.  For 
this  investigation  the  layer  thickness  was  varied  in  incre- 
ments  from  12  in.  to  an  infinite  depth.  From  the  simulation 
of  the  average  deflection"  patterns,  results  (to  be  dis- 
cussed later,  see  Table.  3)  for  the  extreme  conditions  ir.ci- 
cate  that  the  layer  thickness  was  immaterial.  For  simplicity, 
the  soil  mecia  was  assumed  to  be  of  infinite  extent.  Under  ' 


this  assumption,  the  model  characteristics  k and 
(7)  and  (6) , become 


t in  Ecr.s. 


18 


FIGURE  I.  DISTRIBUTION  OF  STRESSES 

WITH  DEPTH 


FIGURE  2.  DISTRIBUTION  OF  DISPLACEMENT 

WITH  DEPTH 
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As  can  be  noted  from  the  above  relations,  the  two  model 
characteristics  are  functions  of  the  unknown  parameters  I 
(modulus  of  elasticity  of  the  soil),  p (Poisson's  ratio  lor 
the  soil) , and  y.  Therefore  to  use  the  developed  mat-soil 
model,  the  parameters  E and  p must  be  numerically  identi- 
fied in  addition  to  theparameter Y . The  identification  of 
these  parameters  could  be  accomplished  by  simulating  the 
"average  deflection"  patterns;  however,  to  consider  the 
effect  of  each  parameter  independently  would  necessitate 
extensive  computer  time. 

Previous  studies  of  beams  on  elastic  foundations  indi- 
cated that  v generally  has  a negligible  effect  on  the  re- 
sulting deflection  pattern.  To  assess  the  validity  of  this 
assertion  for  the  present  problem,  the  deflections  of  two 
points  in  the  mat,  one  directly  under  a wheel  and  one  36.75 
in.  from  the  centerline  of  the  wheel,  were  determined  for 
conditions  identical  to  those  of  Section  13,  Lane  28,  Item 
2 for  Poisson's  ratios  of  0.3  and  0.5.  The  resulting  deflec- 
tions are  tabulated  in  Table  1.  For  this  particular  test, 
a Poisson's  ratio  of  0.4  appears  adequate.  Hence  in  this 
study,  p Q is  assumed  to  be  0.4.  This  assumption  is  also 
in  line  with  that  made  previously  by  Pickett  [ 43,  44]  . In 
any  event,  any  appreciable  error  incurred  by  this  assumption 
would  be  compensated  for,  in  part,  by  the  other  parameters, 

E0  Y * obtained  from  the  simulation  procedure. 
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Table  1 


Effect  of  Poisson's  Ratio  on  Deflection 


Y 

Eo 

(psi) 

Deflection 
at  i of 
Wheel  (in.) 

Deflection 
36.75  in. 
from  i of 
Wheel  (in.) 

Vo  0-3 

3.00 

250 

1.098 

0.099 

Vo  = °-5 

3.00  • 

250 

0.900 

0.054 

Actual^ 

(from  "average 

1.0 

0.0 

deflection" 

pattern) 


. Parameter  Identification  Techniques 

Several  identification  methods  were  applied  to  the 
deflection  data  obtained  from  the  prototype  tests  to  deter- 
mine equivalent  mat-soil  model  parameters.  The  first  method, 
which  has  been  used  successfully  by  mechanical  engineers  for 
similar  systems,  utilizes  a state  variable  filter.  Unfor- 
tunately, for  the  present  problem  it  proved  to  be  totally 
unsatisfactory.  The  second  approach,  the  so-called  "steep 
descent"  method,  proved  to  be  adequate  for  the  identification 
of  the ‘model  characteristic  k,  Eqn.  (12a).  A simplification 
was  introduced  to  this  method  to  reduce  computer  time.  The 
"steep  descent"  method  failed  to  provide  a reliable  measure 
of  the  parameter  y.  Finally,  success  was  registered  in  this 
regard  by  employing  a trial  and  error  procedure. 

For  each  method  and  modification,  input  conditions  were 
imposed  upon  the  developed  mat-soil  model  identical  to  those 
of  the  prototype  tests.  The  criterion  followed  in  the  iden-  * 
tification  of  the  model  parameters  was  to  minimize  the 
response  differences  in  deflection  of  at  least  nine  discrete 
points  taken  from  the  "average  deflections"  patterns.  In 
cases  where  the  curvature  of  the  deflection  pattern  changed 
appreciably  from  point  to  point,  more  reference  points  were 
selected.  In  some  cases  as  many  as  sixteen  reference  points 
were  used. 
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State  Variable  Filter 


The  initial  procedure  utilized  a state  variable  filter 
in  conjunction  with  a steep  descent  linear  identifier  as 
outlined  b / Kohr  [ 30].  The  procedure  (see  Appendix  C)  was 
prcgrammoe'u.-ing  the  field  data;  however,  unrealistic  values 
for  the  parr.nOters  resulted.  Additional  modifications  also 
proved  unsuccessful  and  the  method  v/as  eventually  abandonee. 
Failure  of  this  procedure  to  produce  meaningful  values  was 
believed  to  be  primarily  due  to  the  selection  of  a step 
function  as  the  forcing  function  for  the  load.  The  l?c.<  of 
success  for  a step  load  had  been  noteo  previously  in  ether 
cases  by  Kohr 


step 

[sol . 


Steep  Descent  Method 

The  second  method  employed  to  provide  a measure  of  the 
parameters,  E and  y,  was  based  upon  the  "steep  descent" 
method  (40,  68].  A brief  explanation  of  the  procedural 
aspects  of  this  method  is  given  in  Appendix  D. 

This  procedure  was  initially  applied  to  data  from  the 
test  designated  as  Section  1,  Lane  2,  Item  3,  for  zero 
coverages.  The  initial  values  assumed  for  E and  y were 
100  psi.  and  1.55,  respectively.  Using  the  steep  descent 
procedure,  a minimum  of  the  error  functional  (defined  in 
Appendix  D)  of  0.0S0  was  obtained  when  E ~ 750  psi.  and 
Y = 1.59C?  this  is  indicated  as  trial  1 in  Figure  3.  To 
determine  whether  the  minimum  obtained  was  global  rather 
than  local  another  trial  was  performed.  Trial  2,  Figure  3, 
which  was  initiated  with  E = 200  psi.  and  y = 6.00,  pro- 
duced a minimum  of  0.079  wRen  Eq  = 280  psi.  and y = 6.006. 

With  these  results  it  was  apparent  that  the  surface  of  the 
error  functional  was  not  bowl-like  in  form.  Additional 
trials  were  made  as  indicated  in  Figure  3.  As  can  be  seen 
from  Figure  3,  the  error  functional  possessed  a curved 
valley  of  minimal  values  which  for  all  practical  purposes 
were  identical.  Unfortunately,  the  shape  of  the  error 
functional  proved  not  to  be  unique  as  shown  for  another  test 
ection  in  Figure  4. 

It  was  apparent  that  unique  values  of  E and  y could 
not  be  obtained  with  the  selected  form  of  thR  error  functional. 
Fortunately,  as  can  be  seen  in  Table  2,  the  values  of  the 
parameter  k varied  only  slightly  along  the  valley  of  the 
error  functional.  Therefore,  representative  values  of  the 
characteristic  k could  be  generated  regardless  of  the  value 
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FIGURE  3.  BEHAVIOR  OF  ERROR  FUNCTIONAL 
(SECTION  I,  LANE  2,  ITEM  3 ) 


GAMMA 

FIGURE  4.  BEHAVIOR  OF  ERROR  FUNCTIONAL 
(SECTION  9,  LANE  22  , ITEM  2 ) 


Table  2 


Values 

of  k Along  the  Valley 

of  the 

Error  Functional 

Section  1 

Lane  2 

Item  3 

Zero  Coverage 

' E0 

Y 

Error 

k 

750  psi. 

1.598 

0.060 

52.5  pci. 

530  psi. 

2.620 

0.095 

52.4  pci. 

400  psi. 

3. 608 

0.100 

51.0  pci. 

280  psi. 

6.006 

0.079 

52.8  pci. 

Since 

the  line  of  steepest 

descent 

for  all  trials, 

Figures  3 and  4,  was  essentially  parallel  to  the  Eq  axis, 
the  identification  procedure  was  modified  somev.’hat;  Values 
of  1.0,  2.5,  4.0,  5.5,  7.0,  and  8.5  were  assigned  to  the 
parameter  y an<3  f°r  each  of  these  values  Eq  was  incremented 
until  the  error  functional  was  minimized.  This  procedure 
was  subsequently  employed  for  all  relevant  test  sections  and 
at  all  coverage  levels. 

For  a few  of  the  coverage  levels,  where  the  actual 
deflection  pattern  was  relatively  unsymmetrical , the  magni- 
tude of  the  error  functional  was  found  to  be  greater  than 
1.0.  Nevertheless,  the  representative  value  of  k was  taken 
as  that  which  existed. when  the  error  functional  was  a mini- 
mum for  the  selected  values  of  y. 

With  the  parameter  k defined,  only  one  additional  param- 
eter had  to  be  identified,  As  noted  in  Table  2,  the  param- 
eter E showed  considerable  variation  for  the  range  of  values 
of  y investigated.  As  the  parameter  y appeared  more  stable 
it  was  selected  for  identification. 


Identification  of  the  Parameter  Y 

Due  to  the  small  variation  in  the  error  functional  for 
range  of  Y values  investigated , it  was  concluded  that, 
the  preceding  procedure  was  not  satisfactory  for  identifying 
the  parameter  y • The  procedure  developed  for  determining  y 
was  less  direct  than  that  used  previously  for  the  parameter 
k.  Preliminary  studies,  Figure  5,  indicated  that  the  value 
of  Y did  not  influence  greatly  the  magnitude  of  the  deflec- 
tions. However,  it  was  noted  that  as  y increased  deflections 
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in  tho  vicinity  of  the  loads  became  larger  and  wore  dissi- 
pntod more  rapidly  as  the  distanco  from  the  load  increased. 
This  indicate:;  that  the  value  of  y was  dependent  upon  the 
rigidity  of  the  mat?  that  is,  the  less  rigid  the  mat,  the 
larger  tho  cof lection  under  the  load  and  the  more  rapid  the 
dissipation  of  deflection.  From  comparison  of  deflection 
patterns  at  identical  coverage  levels,  it  was  also  observed 
that-,-  in  general,  the  convexity  (curvature)  of  the  deflec- 
tion patterns  for  the  MG  mat  was  greater  than  that  for  the 
Til  mat.  In  addition,  the  subgrade  strength,  as  represented 
by  the  CBR  values,  was  higher  for  the  M8  mat  tests  than  that 
for  the  Til  mat  tests.  It  was  therefore  concluded  that  the 
difference  in  the  convexity  of  the  deflection  patterns  resul- 
ted primarily  from  the  differences  in  mat  rigidity  and  could 
be  accounted  for  in  the  model  by  the  selection  of  the  param- 
eter y. 


The  parameter  y was  established  by  a trial  and  error 
procedure  utilizing  the  computed  model  deflections.  The 
previously  determined  value  of  k was  maintained  constant 
for  each  coverage  level  while  different  values  were  assigned 
to  y . The  "correct"  value  of  y was  established  by  comparing 
the  computed  model  deflection  configuration  to  the  prototype 
deflection  pattern.  After  the  parameter  y had  been  estab- 
lished slight  modifications  were  made  in  the  value  of  the 
parameter  k to  produce  even  better  correspondence  between 
deflection  patterns.  This  procedure  was  applied  to  all 
tests  to  provide  both  y and  k values. 


Correlation  of  Model  Parameters  to  Mat-Subgrade  Properties 

From  the  results  of  the  identification  procedure  (to  be 
discussed  later,  see  Figures  6 and  7) , it  was  found  that  the 
magnitude  of  k at  any  coverage  level  could  be  established  as 
a function  of  the  initial  value.  The  functional  relationship 
was  established  as 


kINT 
O'  0 4 85" 

N * 


(13) 


where  k is  the  k value  at  zero  coverage  and  k is  the  k 
va!SI  aiSSr  N number  of  coverages.  N 
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In  order  to  determine  KT  and  the  resulting- deflections 
at  any  coverage,  means  for  obtaining  a measure  of  k..:n  hao 
to  be  established.  To  avoic  developing  a new  soil  celt 
which  could  be  both  expensive  and  time  consuming  to  conduct, 
it  war-  decided  to  try  correlating  kTNT  with  the  stanoard 
soil  properties,  water  content,  dry^censity,  CBR,  obtained 
at  the  test  site.  Previous  work  by  Black  [5]  indicated  that 
there  exists  a strong  correlation  between  soil  strength, 
water  content,  and  dry  density.  Since  k is  in  a sense  a 
measure  of  soil  strength,  it  was  thought  unat  there  might 
exist  some  correlation  between  k , w,  and  Y . Representa- 
tive values  of  the  water  content1and  dry  density  were  taken 
as  the  average  of  the  respective  values  given  for  the  top 
18  inches  of  the  subgrade.  Attempts  to  develop  a relation- 
ship among  these  parameters  proved  to  be  fruitless. 

The  literature  [49,  50,  68  ] also  indicated  the  exist- 
ence of  an  empirical  relationship  between  the  subgrade 
modulus  and  CBR.  However,  this  relation  was  established 
for  subgrade  moduli  greater  than  100  pci.  Extension  of 
these  relations  to  prototype  test  conditions  proved  unreli- 
able. Computations  did  disclose  a reasonable  correlation 
among  CBR,  wy  . (weight  of  water  per  unit  volume  of  soil)  and 
the  parameter  R.  It  was  found  for  these  prototype  test  sub- 
grades that  k could  be  reliably  established  from  the 
relation 


kJNT  = 164.0  + 3.0  CBR  - 5.45  w^d  (14) 

where  CBR  is  the  average  CBR  for  the  upper  eighteen  inches 
of  subgrade  and  w ^ is  in  pounds  per  cubic  foot. 

Observations  of  the  "average  deflection"  patterns  dur- 
ing the  prototype  tests  indicate  that  the  curvature,  in 
general,  increased  with  increasing  number  of  coverages. 
Hence,  the  parameter  y must  also  increase  with  coverages. 
The  identification  procedure  demonstrated  this  response: 
the  value  of  Y did  increase  with  coverages.  In  addition  it 
became  slightly  larger  with  decreasing  mat  rigidity.  The 
latter  was  found  to  be  more  pronounced  for  low  coverage 
levels  than  for  high  coverage  levels.  At  all  coverage 
levels,  the  variation  with  mat  rigidity  was  small.  From 
the  simulation  of  the  "average  deflection"  patterns  the 
relationship  for  yN  was  established  as 
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(N  > 1) 


(15) 


* . — ^"ET  + 


(EI)Tu  ^ E ^ MAT 


'L\ 


W) 


MAT 


where  N is  the  specific  number  of  coverages,  (EI)_u  is  the 
r.'.giuity  per  foot  of  width  of  the  Til  aluminum  mat,  and 
(hi),  is  the  rigidity  per  foot  of  width  of  the  mat  being 
investigate--; . The  second  term  of  Eqn.  (15)  has  greater 
influence  on  y for  low  coverages  than  for  high  coverages. 

V.’ith  the  par 
by  Kens.  (14)  end 
the  characteristi 
obtained  the  characteristics  k and  t,  the  mat-soil  model  may 
be  used  to  estimate  the  expected  deflections  for  the  mat  at 
different  coverage  levels. 


ameters  k . and  y defined,  respectively, 
(15),  Eqn.  (12a)  may  be  solved  for  E and 
c t may  be  found  from  Eqn.  (12b) . Having  t 


Results 


‘ Layer  Thickness  Investigation 

Data  from  six  different  test  sections  [12]  were  used 
for  this  investigation.  Identical  conditions  of  wheel  load, 
wheel  spacing,  width  of  tire  print,  and  mat  rigidity  as  used 
in  the  prototype  tests  were  imposed  upon  the  mat-soil  model. 
The  subgrade  thickness  was  varied  from  12  in.  to  an  infinite 
depth.  For  each  layer  thickness  the  simulation  of  the 
"average  deflection"  pattern  was  achieved  by  minimizing  the 
error  functional,  Eqn.  (D2) . The  modified  steep  descent 
method  was  used  with  the  value  of  2.50  assigned  to^the  paranv 
etery  . Values  of  the  error  functional  and  the  model  charac- 
teristic k for  the  extreme  layer  thicknesses  are  given  in 
Table  3. 


Correlation  of  Model  Parameters  to  Mat-Subgrade  Properties 

The  identification  procedure  was  applied  to  all  test 
sections.  Values  of  the  parameters  y,  E , k,  and  the  error 
functional  for  a few  typical  test  sections  are  given  in  < 
Table  4.  For  the  remaining  sections,  the  k values  producing 
the  minimum  error  functional  for  the  selected  values  ofy 
are  given  in  Table  5. 

Plotting  the  magnitude  of  k against  the  number  of  cover- 
ages, Figures  6 and  7,  it  can  be  observed  that  with  only 
rare  exception  the  magnitude  of  k decreases  with  coverage 
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Table  3 


Comparison  of  Simulation  Error  for 
Different  Soil  Layer  Thicknesses 


H = 12 

inches 

H = infinity 

Sect 

Lane 

Item 

Coverages 

Error 

k-pci. 

Error 

k-pci 

1 

1 

3 

0 

.164 

51. 

.167 

51. 

20 

.182 

39. 

.186 

39. 

200 

.278 

40. 

.281 

39. 

300 

.435 

45. 

.439 

45. 

1 

2 

3 

0 

.072 

55. 

.068 

54. 

20 

.064 

45. 

.062 

45. 

40 

.284 

41. 

.278 

41. 

68 

.286 

40. 

.279 

40. 

3 

5 

1 

0 

.105 

17. 

.122 

17. 

30 

1.982 

21. 

.758 

15. 

3 

5 

2 

0 

.176 

23. 

.214 

24. 

30 

.347 

19. 

.349 

19. 

6 

11 

1 

0 

.040 

36 . 

.040 

37. 

600 

.069 

28. 

.070 

28. 

6 

11 

2 

0 

.020 

53. 

.020 

53. 

600 

.017 

47. 

.017 

47. 

30 


Table  4 


Comparison  of  the  Variability  of  the  Characteristic  k 

Section  _1  Lane  1 Item  2, 

Coverage  zero  Coverage  20 


Y 

Kc-pei. 

k-pci . 

Error 

E -psi. 
o 

k-pci . 

Error 

1.0 

1030. 

50. 

0.218 

830. 

37. 

0.i229 

2.5 

530. 

51. 

0.167 

430. 

39.  ‘ 

0 .'<186 

4.0 

380. 

53. 

0.161 

310. 

40. 

0.193 

5.5 

300. 

53. 

0.157 

260. 

44. 

0.368 

7.0 

260. 

56. 

0.181 

210. 

42. 

0.214 

8.5 

220. 

54. 

0.156 

180. 

41. 

0.184 

Coverage  200 

Coverage  300 

Y 

EQ-psi. 

k-pci . 

Error 

E -psi. 
o 

k-pci . 

Error 

1.0 

840. 

38. 

0.334 

940. 

44. 

0.513 

2 . 5 

' 430. 

39. 

C . 281 

480. 

45. 

0.439 

4.0 

310. 

40. 

0.271 

340. 

46. 

0.429 

5 5 

260. 

43. 

0.404 

270. 

46. 

0.425 

7 0 

210. 

42. 

0.278 

230. 

47. 

0.421 

8.5 

180. 

41. 

0.265 

210. 

51. 

0.512 

Section  1 

Lane  2 

Item  3 

Coverage  zero 

Coverage  20 

Y 

EQ-psi. 

k-pci. 

Error 

E -psi. 
o 

k-pci. 

Error 

1 0 

1070. 

52. 

0.052 

930. 

44. 

0.057 

2 5 

550. 

54. 

0.068 

480. 

45 . 

0.061 

4 0 

390 . 

55 . 

0.075 

340  - 

46. 

0 . 067 

hr  • v 

310. 

55. 

0.077 

270. 

46. 

0.070 

7 O 

260. 

56. 

0.079 

230. 

47. 

0.072 

8.5 

230. 

57. 

0.079 

210. 

51. 

0 . 167* 
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Table  4 (Cont'd.) 


Section 

1 

Lane  2 

Item  3 

Coverage  40 

Coverage  68 

Y 

vpsi 

k-pci. 

Error 

E -psi. 
0 

k-pci. 

Error 

1.0 

870. 

40. 

0.262 

850. 

39. 

0.242 

2.5 

450. 

41. 

0.278 

440. 

40. 

0.279 

4.0 

320. 

42. 

0.286 

310. 

40. 

0.294 

5.5 

260. 

44. 

0.297 

260. 

44. 

0.369 

7.0 

220. 

45. 

0.305 

210. 

42. 

0.299 

8.5 

210. 

51. 

0.682 

210. 

51. 

0.971 

Section  3 

Lane  5 

Item  1 

Coverage  zero 

Coverage  30 

Y 

EQ-psi.  k-pci. 

Error 

E -psi. 
o 

k-pci. 

Error 

1.0 

430.  16. 

0.116 

380. 

13. 

0.433 

2.5 

230.  17. 

0.122 

210. 

15. 

0.758 

4.0 

160.  17. 

0.135 

200. 

22. 

1.693 

5.5 

200.  31. 

3.189 

120. 

16. 

0.672 

7.0 

110.  18. 

0.107 

150. 

27. 

1.950 

8.5 

150.  32. 

1.453 

150. 

32. 

6.370 

Section  3 

Lane  _5 

Item  2 

Coverage  zero 

Coverage  30 

Y 

EQ-psi.  k-pci. 

Error 

E0-psi 

. k-pci. 

Error 

1.0 

480.  22. 

0.222 

413. 

18. 

0.390 

2.5 

260.  24. 

0.214 

22\ 

19. 

0.349 

4.0 

180.  24. 

0.174 

160. 

20. 

0.379 

5.5 

160.  27. 

0.601 

150. 

25. 

2.065 

7.0 

120.  24. 

0.196 

110. 

21. 

0.433 

8 . 5 

110.  26. 

0.194 

150. 

39. 

2.374 

Section  6 

Lane  11 

Item  1 

Coverage  zero 

Coverage  600 

Y 

EQ-psi.  k-pci. 

Error 

Eo-psi 

. k-pci. 

Error 

1.0 

770.  34. 

0.038 

640. 

26. 

0.077 

2 . 5 

410.  37. 

0.040 

340. 

28. 

0.070 

4.0 

290.  37. 

0.041 

240. 

29. 

0.069 

5 . 5 

230.  37. 

0.041 

190. 

29. 

0.069 

7.0 

190.  37. 

0.043 

160. 

29. 

0.069 

8 . 5 

170.  38. 

0.041 

140. 

30. 

0.069 
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Table  5 


1;  Valuer  for  .'lininmi.i  Simulation  Error 


Sect 

Lane 

Item 

Coverage 

k-pci. 

Error 

2 

3 

1 

0 

38.4 

0.102 

200 

28.8 

0/023 

600 

28.7 

0.040 

2 

3 

2 

0 

42.6 

0.172 

66 

32.2 

0.872 

120 

57.4 

o.oe6 

200 

30.3 

0.144 

2 

4 

1 

0 

19.9 

0.177 

20 

16.5 

0.373 

2 

4 

2 

0 

33.0 

0.130 

20 

32.2 

2.836* 

3 

6 

1 

0 

12.3 

0.174 

20 

10.5 

0.211 

76 

15.0 

0.661 

156 

13.9 

0.298 

3 

6 

2 

0 

21.8 

0.308 

20 

15.9 

0.841 

76 

18.4 

5.818* 

6 

11 

2 

0 

50.2 

0.014 

600 

48.5 

0.016 

6 

11A 

1 

0 

32.1 

0.137 

20 

27.0 

0.259 

130 

26.4 

0.168 

6 

11A 

2 

0 

48.8 

0.425 

20 

38.0 

0.623 

130 

32.6 

1.951* 

6 

12 

1 

0 

21.6 

3.689* 

20 

18.5 

0.646 

90 

21.6 

0.203 

6 

12 

2 

0 

60.0 

0.152 

20 

34.1 

0.253 

44 

34.3 

0.669 
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Table  5 (Cont'd.) 


Sect 

Lane 

Ite:n 

Coverage 

h-pci. 

Erroi 

Q 

21 

1 

0 

18.9 

0.070 

V 

20 

15.6 

0.159 

200 

12.6 

0.131 

300 

11.3 

0.541 

600 

10.4 

0.868 

a 

21 

2 

0 

29.1 

0.053 

V 

20 

26.7 

0.041 

200 

25.9 

0.500 

300 

21.8 

0.126 

9 

22 

1 

0 

13.9 

0.511 

20 

15.3 

1.266* 

400 

9.6 

7.937* 

9 

22 

2 

0 

16.8 

0.573 

20 

16.8 

0.425 

100 

16.8 

0.927 

10 

23A 

1 

0 

42.9 

0.222 

32 

32.9 

0.160 

10 

23^ 

0 

*- 

2 

32.2 

3.825* 

10 

23B 

1 

4 

32.1 

0.586 

13 

28 

1 

0 

14.0 

0.687 

200 

8.9 

0.992 

550 

5.6 

2.890* 

13 

28 

2 

0 

25.3 

0.145 

200 

21.1 

0.336 

550 

20.5 

0.575 

13 

29 

1 

0 

22.7 

0.395 

42 

7.8 

2.20^* 

140 

5.6 

7.7°4* 

13 

29 

2 

0 

26.1 

0.384 

42 

19.9 

0.542 

140 

19.9 

0.432 

200 

19.9 

1.714* 

* 

Indicates  that 

the  deflection  pattern  as  given 

1 12  1 is 

significantly  unsymmctrical . 
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100 


COVERAGES 

FIGURE  6.  VARIATION  OF  PARAMETER  k WITH  COVERAGE 
FOR  SINGLE-WHEEL  TESTS. 


for  both  the  single-wheel  and  dual-wheel  tests.  In  these 
illustrations  to  accommodate  the  log  scale,  the  value  ob- 
tained by  simulating  the  "average  deflection"  pattern  for 
zero  coverages  has  been  plotted  as  occurring  at  one  coverage. 
Also  shown,  is  the  slope  used  for  establishing  the  functional 
relationship  for  kVT,  Eqn.  (13) . For  clarity  of  presentation, 
results  for  seven  Available  sections  were  omitted  from  these 
figures.  As  can  be  observed  from  Table  5,  they  exhibited 
similar  behavior  to  those  shown  plotted. 

Deflection  patterns  were  calculated  incorporating  the 
refined  values  of  k and  Eqns.  (13)  and  (15).  Some  typical 
deflection  patterns^are  presented  in  Table  6 along  with  the 
actual  deflection  patterns.  The  distances  given  in  this 
table  are  referenced  from  a point  ten  feet  left  of  the  cen- 
terline Of  the  traffic  lane. 

Table  7 provides  a comparison  of  the  initial  subgrade 
modulus  as  given  by  Eqn.  (14)  with  simulated  values.  It  is 
noted  that  the  largest  discrepancies  occur  for  those  sections 
where  the  simulated  values  are  generally  high.  In  Figure  8, 
it  is  observed  that  small  variations  in  the  subgrade  strength, 
as  reflected  by  the  subgrade  modulus,  have  appreciable  influ- 
ence on  the  deflection  characteristics  of  the  load  transfer 
element.  Consequently,  Eqn.  (14)  was  developed  with  a bias 
directed  toward  the  better  reproduction  of  low  values  of 

kINT* 


Summary  and  Discussion  of  Model 

As  indicated  by  the  results  presented  in  Table  3,  any 
choice  of  layer  thickness  greater  than  twelve  inches  appears 
to  have  negligible  effect  upon  the  simulation  of  the  average 
deflection"  patterns.  Similar  behavior  was  also  observed  by 
Lenczer  [36  f in  his  study  of  the  influence  of  layer  thickness. 
The  assumption  of  an  infinite  depth  for  the  subgrade  was  both 
expeditious  and  justifiable. 

indications  were  that  a Poisson’s  ratio  of  0.4  for  the 
subgrade  material  was  reasonable;  selection  of  this  valu 
was  based  upon. the  work  of  Pickett  [ 43,  44 J and  the  results 
of  the  preliminary  study,  Table  1. 
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Table  6 


Comparison  of  Deflection  Patterns 


Section  1 Lane  1 

Item  3 

k..._ 

= 50.5pci. 

YT1l.m  = 

1.00 

IET 

INT 

aro  Coverages  - Error  = 0.059 

20  Coverages  - Error  = 0.059 

Actual 

Estimated 

Actual 

Estimated 

Deflection 

Deflection 

Deflection 

Deflection 

[stance 

(in) 

(in) 

Distance  (in) 

(in) 

(in) 

(in) 

47. 

0.00 

0.07 

44. 

0.00 

0.01 

51. 

0.10 

0.14 

51. 

0.20 

0.14 

56. 

0.22 

0.26 

56. 

0.35 

0.29 

66 . 

0.75 

0.66 

66. 

0.80 

0.78 

76. 

1.10 

1.22 

69. 

1.00 

0.96 

86. 

1.55 

1.65 

76. 

1.40 

1.43 

96. 

1.79 

1.75 

86. 

1.85 

1.94 

106. 

1.59 

1.65 

96. 

2.20 

2.08 

116. 

1.18 

1.22 

106. 

1.90 

1.94 

119. 

1.00 

1.05 

116. 

1.30 

1.43 

126. 

0.75 

0.66 

123. 

1.00 

0.96 

136. 

0.31 

0.26 

126. 

0.85 

0.78 

146 . 

0.00 

0.05 

136 . 

0.31 

0.29 

148. 

0.00 

0.01 

Coverages  - Error 

= 0.313 

68  Coverages  - Error 

= 0.348 

Distance 
(in) 

41. 
51. 
56. 
67. 
76. 
84. 
96. 
106. 
116. 
126. 
131. 
141. 
151. 


Actual  Estimated 
Deflection  Deflection 


(in) 


0.00 

0.31 

0.38 

1.00 

1.50 

2.00 

2.50 
2.00 
1.20 
0.92 
0.50 
0.35 
0.00 


(in) 


-0.03 

0.15 

0.31 

0.87 

1.48 

1.93 

2.15 

2.01 

1.48 

0.81 

0.53 

0.15 

•0.03 


Actual 
Deflection 
Distance  (in) 


(in) 

49. 

56. 

66. 

76. 

87. 

96. 

109. 

116. 

121. 

126. 

136. 

141. 

151. 


0.00 

0.48 

1.10 

1.40 

2.00 

2.50 

2.00 

1.48 

1.38 

0.90 

0.50 

0.35 

0.00 


Estimated 

Deflection 

(in) 


0.10 

0.32 

0.83 

1.52 

2.09 

2.20 

1.94 

1.52 

1.17 

0.83 

0.32 

0.15 

-0.03 
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Table  6 (Cont'd.) 


Section  3 Lane  5 

Item  2 

kINT  * pci 

* \liT 

' 1.75 

Zeie  Cove 

raj on  - Error 

•■=  0.621 

30  Coverages  - Error  **  0.619 

Actual 

Estimated 

Actual 

Estimated 

Deflection 

Deflection 

Deflection  Deflection 

Distance 

(in) 

(in) 

Distance  (in) 

(in) 

(in) 

(in) 

99. 

0.00 

0.47 

96. 

0.00 

0.42 

104. 

0.50 

0.73 

104. 

0.79 

0.88 

109. 

1.00 

1.04 

114. 

1.72 

1.64 

119. 

1.71 

1.73 

119. 

2.00 

2.03 

124. 

2.16 

2.01 

124. 

2 . 66 

2.35 

134. 

2.10 

2.20 

135. 

2.80 

2.60 

144* 

1.99 

2.17 

144. 

2.60 

2.58 

154. 

2.11 

2.00 

154. 

2.60 

2.60 

166. 

2.00 

1.91 

164. 

2.35 

2.35 

174. 

1.37 

1.39 

169. 

1.67 

2.03 

183. 

0.50 

0.79 

174. 

1.60 

1.64 

190. 

0.00 

0.42 

184. 

0.68 

0.88 

194. 

0.00 

0.34 

Section 

13  Lane  28 

Item  2 

kTNT 

25.0  pci. 

yint 

1.75 

Zero  Coverages  - Error 

= 0.157 

550  Coverages  - Error  = 0.788 

Actual 

Estimated 

Actual 

Estimated 

Deflection 

Deflection 

Deflection 

Deflection 

Distance 

(in) 

(in) 

Distance 

(in) 

(in) 

(in) 

(in) 

» 

32 

-0.16 

0.05 

32. 

-0.16 

0.09 

36. 

0.00 

0.10 

42. 

0.15 

0.32 

42. 

0.13 

0.22 

52. 

0.52 

0.70 

62. 

1.00 

0.86 

57. 

1.00 

0.92 

72. 

1.00 

1.06 

72. 

1.35 

1.38 

82. 

1.00 

0.93 

82. 

1.30 

1.24 

92. 

0.73 

0.64 

94. 

1.00 

0.88 

102. 

0.51 

0.52 

102. 

0.82 

0.79 

112. 

0.68 

0.65 

114. 

1.00 

0.98 

122. 

1.00 

0.92 

122. 

1.00 

1.23 

132. 

1.02 

1.06 

132. 

0.83 

1.38 

142. 

1.00 

0.86 

142. 

0.81 

1.14 

152. 

0.40 

0.51 

152. 

0.30 

0.70 

162. 

0.12 

0.22 

170. 

0.00 

0.12 

167. 

a. oo 

0.12  39 

Table  7 


Comparison  of  Simulated  to  Calculated  Values  of  k^ 


Test 

k 

pci. 

Sect 

Lane 

Item 

Simulated 

Calculated 

1 

1 

3 

48.5 

59.3 

1 

2 

3 

50.5 

66.1 

2 

3 

1 

35.5 

18.5 

2 

3 

2 

41.0 

35.1 

4 

1 

19.0 

15.9 

2 

4 

2 

31.5 

38.3 

3 

5 

1 

16.0 

17.0 

3 

5 

2 

21.0 

32.7 

3 

6 

1 

12.5 

20.4 

3 

6 

2 

18.0 

29.2 

6 

11 

1 

34.0 

20.6 

6 

11 

2 

50.5 

41.8 

6 

12 

1 

22.0 

29.4 

6 

12 

2 

36.0 

36.7 

9 

21 

1 

14.0 

11.9 

9 

21 

2 

25.5 

20.2 

9 

22 

1 

10.0 

11.2 

9 

22 

2 

16.0 

13.5 

10. 

23A 

1 

37.5 

27.1 

10 

2 3A 

2 

27.5 

31.7 

10 

23B 

1 

32.5 

27.1 

13 

28 

1 

11.0 

11.6 

13 

28 

2 

25.0 

25.3 

13 

29 

1 

18.0 

11.7 

13 

29 

2 

28.0 

24.2 
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FIGURE  8.  INFLUENCE  OF  k|NT  ON  MAXIMUM  DEFLECTION 


It  is  believed  that  failure  of  the  state  variable  filter 
and  its  modifications  to  identify  the  model  characteristics 
k and  t vas  due  to  the  intrinsic  representation  of  the 
applied  loads  as  step  functions.  Cimilar  difficulties  with 
step  loads  had  been  noted  earlier  by  Kchr  [30] . 

The  model  characteristics  k and  t,  Eqns.  (12a)  and  (12b), 
are  functions  of  y in  addition  to  the  conventional  Eq  and 
parameters.  The  "steep  descent"  method,  in  which  the  deflec- 
tion difference  between  the  model  and  the  prototype  for  at 
least  nine  discrete  points  was  minimized,  also  failed  to 
produce  unicrne  values  of  the  parameters  Eq  and  y , Figures  3 
and  4.  However,  the  error  functional  Eqn.  (D2) , was  found 
to  possess  a valley  of  minimums  along  which  the  value  of  the 
characteristic  k was  essentially  constant. 

The  simulation  of  the  "average  deflection"  patterns, 

Table  4,  indicated  that  the  value  of  the  model  parameter  k 
increased  slightly  with  increases  in  Y.  For  a specific  value 
of  y (as  shown  in  Table  4)  the  magnitude  of  k was  found  to 
decrease  with  the  number  of  coverages.  Computations  indi- 
cated that  the  magnitude  of  k was  more  sensitive  co  the 
number  of  coverages  than  to  the  value  of  the  parameter  Y . 

The  evident  validity  of  the  developed  expression  for  k , 

Eqn.  (13)  is  amply  demonstrated  in  Tables  4 and  5 and  Figures 
6 and  7. 

As  noted  in  Figure  5,  variations  in  the  parameter  Y were 
reflected  primarily  as  alterations  of  the  deflection  pattern 
curvature.  This  type  of  response  coulc  not  be  identified  by 
the  selected  form  of  the  error  functional,  Eqn.  (D2) . The 
representative  value  of  Y was  established  from  the  similarity 
of  model  deflection  curvature  for  various  values  of  y with 
prototype  deflection  curvature.  From  this  comparison  it  was 
noted  that  the  curvature,  thus  Y,  increased  with  coverages 
and  decreased  with  increasing  mat  rigidity.  This  behavior 
is  expressed  by  Eqn.  (15). 

After  the  parameter  Y was  established  (by  the  above 
procedure) , slight  modifications  were  made  in  the  magnitude 
of  k to  provide  even  better  correspondence  between  the  deflec- 
tion patterns.  The  simulated  values  of  the  characteristic  k 
at  zero  coverages,  Table  7,  were  found  to  be  less  than  51  pci. 
In  this  range,  the  model  deflections  were  found  to  be  cruite 
sensitive  to  the  magnitude  of  k . For  an  incremental 
change  in  k_  , the  difference  Irrdef lection  characteristics 
was  found  to  increase  as  k decreased,  Figure  8. 
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An,- empirical  relationship,  Eon.  (14),  was  established 
to  relate  k m to  tho  given  coil  properties  of  CBR,  w,  and 
y . Due  to  J"Uhe  above  described  model  sensitivity,  this 
relationship  was  developed  with  a bias  toward  the  smaller 
valuer;  of  . 

In  all  cases  where  tho  "average  deflection"  patterns 
were  fairly  symmetrical,  the  developed  mat-soil  model  was 
able  to  duplicate  tho  mat  behavior  under  static  loads  with 
a reasonably  lev;  degree  of  error.  The  degree  of  reliability 
of  the  simulation  of  prototype  behavior,  as  reflected  by  the 
magnitude  of  the  error  functional,  can  be  inferred  from 
Table  6.  The  magnitude  of  the  error  reflects  the  number  of 
points  used  for  simulation.  For  example,  assuming  that  the 
total  error  is  uniformly  distributed  among  all  points, 
errors"of  0.5625  and  0.6250  are  equivalent  to  a quarter  inch 
difference  at  each  point  on  deflection  patterns  represented 
by  nine  and  ten  points,  respectively.  It  is  further  noted 
that  only  a few  deflection  patterns  were  described  by  the 
minimum  of  nine  points.  The  results,  Table  5 and  particu- 
larly Table  6,  demonstrate  the  general  reliability  of  the 
chosen  model  and  the  latter  table  confirms  Eqns.  (13)  and 
(15). 


For  those  coverage  levels  where  the  resulting  error  was 
greater  than  unity,  the  prototype  deflection  pattern  was 
qrossly  non-symmetrical  [12].  Behavior  of  this  type  could 
not  be  simulated  with  the  chosen  mat-soil  model  because 
imposed  conditions  assumed:  1)  that  the  applied  load  was 

symmetrical  and  2)  that  the  subgrade  was  homogeneous. 


DEVELOPMENT  OF  FAILURE  CRITERION 


Introduction 


A failure  criterion  had  to  be  established  to  use  the 
developed  mat-soil  model  as  a means  of  predicting  mat  per- 
formance. The  search  for  a failure  criterion  which  would 
relate  well  to  performance  produced  many  blind  alleys. 
Finally,  success  was  had  through  an  educated  "trial  and 
error"  procedure.  The  full  chronicle  of  procedures  and 
methods  examined  will  be  given  below  in  the  hopes  that 
their  display  will  prove  valuable  to  others  on  similar 
expeditions . 

Initially,  the  failure  criterion  employed  was  the 
roughness  criterion  established  by  the  Corps  of  Engineers 
[12].  The  basis  of  this  failure  criterion  was  a deviation 
in  the  mat  sucface  of  3 inches  or  more  within  a 10  feet 
length  in  any  direction  within  the  traffic  lane.  Measure- 
ments were  always  made  on  the  unloaded  mat  and  the  con- 
figuration of  the  mat  surface  at  intermediate  coverage 
levels  and  at  failure  were  recorded  as  "average  cross- 
sectional  deformations"  [12].  These  configurations  are 
referred  to  as  "deformation"  patterns  in  the  present  study 
and  they  represent  the  residual  or  irrecoverable  displace- 
ments of  the  mat. 

The  established  roughness  criterion  depended  in  some 
complicated  way  upon  the  residual  deformation  of  the  sub- 
grade and  the  permanent  set  that  occurred  in  the  mat  element. 
As  the  selected  mat-soil  model  was  capable  of  estimating 
deflections  only  when  loads  were  applied,  a fictitious  load 
was  assumed  to  act  on  the  mat.  The  distribution  of  this 
load  was  taken  to  be  similar  to  the  distribution  of  the 
traffic  imposeo  during  each  test.  The  magnitude  of  this 
fictitious  load  was  selected  to  be  that  necessary  to 
simulate  the  maximum  deformation.  Investigation  of  the 
model  cieflection  curvature  revealed  that  the  roughness 
criterion  could  not  be  satisfied  at  the  failure  coverage; 
thus  this  approach  was  deemed  not  successful  and  another 
procedure  was  sought. 
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A dimensional  model  was  selected  as  the  second  means 
of  establishing  a failure  criterion  for  the  mat'or  of 
indicating  the  dominate  factors.  Attempts  to  develop  a 
functional  relationship  between  the  failure  coverage  and 
tnc  tost  parameters  indicated  that  the  developed  functional 
i elationshi.p  was  not  a reliable  measure  of  performance. 

Next,  the  mat-soil  model  deflection  at  the  indicated 
failure  coverage  was  expressed  as  a product  function  of  the 
fail  lure  coverage  ana  test  parameters.  A stepwise  regression 
analysis  of  the  data  in  this  form  was  undertaken.  The  data 
were^ divided  into  two  categories:  single-wheel  tests  and 

dual— wheel  tests.  Unfortunately,  for  both  categories,  the 
regression  analysis  indicated  that  the  number  of  coverages 
to  failure  was  an  insignificant  factor. 


In  light  of  the  above  study,  various  parametric  com- 
binations were  analyzed  in  an  effort  to  obtain  an  ordering 
which  would  be  more  dependent  upon  the  number  of  coverages 
to  failure,  Nf.  By  adjusting  the  exponents  occurring  in 
the  product  function,  a combination  was  found  for  the  dualr- 
wheel  tests  which  proved  reasonably  reliable  as  a failure 
criterion.  Success  for  the  single-wheel  tests  was  limited. 


Residual  Load  Concept 

Comparison  of  the  observed  "deformation"  patterns  with 
the  observed  "average  deflection"  patterns  revealed  that, 
although  the  magnitude  of  displacements  in  both  patterns’ 
were  Relatively  the  same,  the  curvature  was  much  less  pro- 
nounced in  the  "deformation"  pattern.  The  primary  differ- 
ence between  the  "deformation"  pattern  and  the  "average 
deflection"  pattern  was  assumed  to  be  attributed  solely  to 
the  manner  of  loading.  On  this  basis  the  model  parameters 
obtained  from  simulation  of  the  "average  deflection"  pattern 
were  equally  applicable  for  the  simulation  of  the  "deforma- 
tion" patterns. 

Two  types  of  traffic  distribution,  uniform  and  nonuni- 
form, were  used  during  the  prototype  tests.  For  the  tests 
with  uniform  distribution,  the  fictitious  load  on  the  model 
was  applied  uniformly  over  the  traffic  lane.  For  the 
nonuniforrn  traffic  distribution,  the  distribution  of  the 
fictitious  load,  q^,  was  taken  as  shown  in  Figure  9. 
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PERCENT  OF  qf 


FIGURE  9.  TRAFFIC  DISTRIBUTION 

ACROSS  TRAFFIC  LANE 


determined  by  simulating  tho 
prototype  tost.  Tho  range  of 
large,  in  addition,  ae  noted 
tho  deflection  patterns, 
much  less  than  that  of  tho 


The  magnitude  of  a,  './os 

maximum  ch' formation  for**ach 

va.lu-.r  of  r v;ac  found  to  bo 

proviouj:  ly , ‘ tV.  j curvature  o l 
obtained  fro  .\  tho  model,  v/aa 

"do formation”  jx-tferno.  This  was  particularly  noted  in 
tho  Inner  which  were  uniformly  loaded.  Consequently,  the 
three  inr'.x  roujhnevs  criterion  could  not  be  met  at  a 
coverage  level  which  could  be  considered  to  produce  failure. 
Thus,  this  approach  for  establishing  a failure  criterion 
was  abandoned . 


Dimensional  Analysis  and  Taylor  Expansion 

Another  approach  followed  to  define  failure  was  to 
obtain  a functional  relationship  among  the  mat,  soil*  loader 
ing  properties,  and  the  number  of  coverages  which  caused 
failure.  It  was  thought  that  a functional  relationship 
of  the  form 


N = f(WL,  CA,  El,  w^,  CBR)  (16) 

v? 

where  WL  is  the  single  wheel  load  in  kips 
CA  is  the  contact  area  in  sq.  in. 
and  the  other  symbols  are  as  previously  defined 

would  provide  some  insight  to  the  development  of  a failure 
criterion  or  could  possibly  even  serve  as  a means  of  predicting 
performance.  To  Reduce  the  number  of  variables,  the 
Buckingham  Pi  theorem  [32]  was  applied  to  the  system,  ■ 

allowing  Eqn.  (16)  to  be  expressed  as 


N * f (*j , *2  , w3  ) (17) 

where  the  it's  are  dimensionless  ratios  or  numbers.  Letting 
ttj  be  a dimensionless  ratio  of  wy^,  WL,  and  El, 

Ci  q, 

*i  = (wyd)  (WL)  (El)  (1G) 

where  Cj  and  are  constants,  it  was  found  that 
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2/3 


(EI) 


*1  = 


(::l)5/2 


(19) 


Similarly,  letting  »2  be  a dimensionless  ratio  of  v;  ■£,  T..X, 
and  CA, 

. <'-v 2/3 


(CA) 


2 ~ 


(V?L)2/3 

as  CBR  is  a dimensionless  quantity, 


(20) 


*3  = CBR  (21) 

With  the  * terms  so  designated,  the  pursuit  of  a reasonable 
functional  relationship  followed . The  function  was  expanced 
into  a Taylor  series  about  the  fixed  points  Wj  , *2  , and  »3 , 
which  were  assigned  typical  values  for  the  tj,  *2,  and  n3 
factors.  The  expansion  of  Eqn.  (17)  lead  to  the  following 
relationship 


f (*i  , "3)  = f (ir°  , u®,*3  ) 


f(*i'  v'j1 
3 "I2 

ffw — J f(ltl  ' V 


0 0 0 

ir,  »ir  <n, 
1 2 J 


+ (,2-wi>ST-+(»  -O: 


iT/ir  -M  JL 
1 ‘’h, 

57  [(’l  -'I’Sp'1’'!  -’'2)3^~+‘,3  - ’?tr]  1 ' ’3> 


3 

+ . . . 


(22) 


0 , 0.0 
1 2 3 


0 0 0 
*1  '*2  ' *3 


where  f ( ir°,  n°)  represents  the  value  of  the  function  at 

the  fixed  points3  *°  , and*0  and  hence  is  a constant. 

1 2 3 
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The  expression 


f (ifj 


i ^ 2 1 


,)l 


0 

*1 


0 

' *2 


0 

*3 


(i  * 1,  2,  3) 


represents  the  partial  derivative  of  the  functiog  with  0 
respect  to  it.  evaluated  at  the  fixed  points  *i  , > and  w3  . 

All  of  the  partial  derivatives  are  unknown  but  once  fixed 
points  are  selected  they  become  constants. 

Data  were  used  from  the  prototype  tests  for  the  single- 
wheel  load  assemblages  wherein  failure  occurred  before  700 
coverages.  In  all,  a total  of  26  different  sets  of  data’ 
were  utilized.  V7ith  this  amount  of  available  data  as  many 
terms  of  the  Taylor  series  were  used  as  possible.  Considering 
only  the  first  4 terms  of  the  series,  20  different  coefficients 
or  constants  had  to  be  evaluated  to  obtain  the  required 
functional  relationship.  If  the  functional  relationship 
was  to  apply,  of  necessity,  it  had  to  satisfy  each  set  of  , J: 
data.  Thus  in  order  to  determine  the  coefficients,  the 
data  for  20  different  tests  were  used  to  produce  20  differ- 
ent relationships  which  were  then  solved  simultaneously  to 
produce  the  20  required  coefficients. 


All  of  the  dual-wheel  tests  which  produced  failure  were 
also  investigated  in  the  same  manner  with  the  inclusion  of 
the  wheel  spacing,  WS,  as  an  additional  variable.  This 
additional  variable  necessitated  an  additional  T term. 
Letting  ir4  be  the  dimensionless  ratio  between  wY  , WL,  and 
WS,  it  was  found  that 


(wY  J 


V3 


(ws) 


4“ 


(WL) 


T7T 


(23) 


Hence,  for  dual-wheel  tests  the  functional  relationship 
becomes 


N = f (w  , *2  , *3,*4  ) (24) 

where  * , tt  , and  tt.  are  as  defined  previously.  Due  to  the 
limited  1 amount  of  data  available,  only  the  first  three  terms 
of  the  series  expansion  could  be  considered.  This  necessi- 
tated the  evaluation  of  15  coefficients. 
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Data  from  20  different  tests  were  used  to  develop  the 
function?!  relationship  for  the  single-wheel  tests,  Den.  (22). 
Tre  functional  relationship  obtained  from  one  particular 
combination  of  20  data  sots  was 


N = - 62.200 

- 0 . 154 ( *2-  nj 

- 0.001(  tt®) 

+ 2.042  ( it®)  (w  - 

* 22  3 3 

t 0.020  ( *2-  hJ) 

- 0.006  ( *j-  *®)  (*2-*2)  (n3 

2 

- 0.003  3 -*3) 

- 0.016  («2- (ir3-»J) 

2 

- 3.602  (»2-  \)  (»3-*J) 

- 0.000(if2-  it®? 


+ 0.265  ( Uj-  w?) 

- 85 . 261  ( t<3  - tt®) 

+ 1.070  ( Tij-n®  ) (tt3  -tt®) 

+ 0.002(  itj-ir®  ) 

o2 

+ 115 .730  (ir3  - ir°)  (2 

-n®)  - 0.000(  7 (7r2-ir2) 

- 0.000(  ffj-lfj  ) (t72  - 77j) 

+ 1.036  ( 77  j -*®  ) (773-7*3) 

- 0.000(  77^®  ) 

+ 8 . 466  ( Tfj-if®  )3 


Thin  relationship  was  applied  to  all  26  available  field 
tests.  Table  8 gives  the  results. 


As  expected,  the  functional  relationship  predicted  the 
exact  failure  coverage  for  the  20  sets  of  data  used  in  its 
development . However,  the  prediction  uf  the  failure  coverage 
was  quite  poor  for  the  other  tests.  A total  of  100  differ- 
ent combinations  of  data  sets  were  investigated.  The  coeffi- 
cients obtained  from  the  different  data  combinations  were 
very  erratic. 


The  same  procedure  was  also  applied  to  the  dual-wheel 
test  data.  Again  the  developed  functional  relationship 
failed  to  reproduce  the  failure  coverage.  Hence  this  pro- 
cedure was  abandoned . 


Regression  Analysis 


Actual  wheel  loadings  were  applied  to  the  mat-soil  model 
using  Eqns.  (13),  (34),  ana  (15),  and  the  deflection  of  the 

mats  was  computed  for  the  failure  coverage.  In  the  single— 
uheel  tests , the  deflection  of  the  mats  at  the  center  of  the 
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Table  8 


Results  of  Dimensional  Analysis  Study 
of  Single-Wheel  Tests 


Test  Coverages 


Sect . 

Lane 

Item 

Actual 

Failure 

Predicted 

Failure 

2 

3 

1 * ** 

600. 

600. 

2 

3 

2 * 

120. 

120. 

9 

21 

2 * 

300. 

300. 

6 

11A 

1 * 

130. 

130. 

6 

11A 

2 * 

130. 

130. 

10 

23A 

1 * 

32. 

32. 

10 

23A 

2 * 

2. 

2. 

10 

2 3B 

1 * 

4. 

4. 

10 

23B 

2 * 

2. 

2. 

7 

★ * 

94. 

94. 

8 

★ * 

6. 

6. 

9 

** 

94. 

94. 

36 

** 

14. 

14. 

37 

** 

98. 

98. 

38 

** 

197. 

197. 

39 

*★ 

60. 

60. 

40 

** 

30. 

30. 

41 

** 

30. 

30. 

71 

** 

160. 

466. 

72 

** 

675. 

11384. 

73 

** 

473. 

5428. 

74 

** 

74. 

74. 

75 

** 

360. 

7292. 

76 

** 

223. 

738. 

92 

** 

300. 

300. 

93 

** 

23. 

156. 

* Reference  12 

**  Reference  10 
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wheel  was  recorded  and  is  presented  in  Table  9 and  Figure  10. 
For  the  dual-wheel  test  sections,  the  deflection  of  the  mats 
at  the  centerline  of  the  wheel  assembly  was  recorccd  ar.c  is 
presented  in  Tabic  10  and  Figure  11.  It  was  anticipated 
that  a failure  criterion  could  be  established  empirically 
providec  that  the  results,  as  shown  in  Figures  10  and  11, 
could  be  grouped  in  some  manner.  Since  the  predicted  deflec- 
tions were  themselves  functions  of  other  test  parameters,  in 
acoition  to  coverages,  the  predicted  cef lections  for  each 
test  were  taken  as  a product  function  of  these  parameters. 

For  a single-wheel  test,  this  relation  was  expressed  as 

|=  (WLR)Cl  (EIR)C2  (°VC3  (CBR)C4nCs  (26) 

where  A is  the  predicted  deflection  at  the  failure 
coverage  in  inches 

WL_  is  the  wheel  load  relative  to  35  kips 
Eld  is  the  rigidity  per  foot  of  width  relative  to 
* 13680  kip-in 

CA  is  the  contact  area  relative  to  150  sq.  in. 

N is  the  number  of  coverages  to  failure 

Cl  , Cz  i C 3 , C4,  C5  are  exponents  whose  values  are 

to  be  determined. 

The  parameters  of  tire  inflation  pressure  and  tire  ply 
rating  wore  not  included  because  their  effects  were  con- 
sidered to  be  reflected  by  the  parameters  WLD  and  CA„. 

The  data  from  each  single-wheel  test  section  was 
written  in  the  form 


logAf=  c xlcg  (WLr)  +C£  log  (EIR)  +C3  log  (CAR)  -K^  log  (CBR)  +C  logs 

(27) 

A multiple  stepwise  regression  analysis  of  the  data  was 
performed  using  a computer  program  [ Cl]  to  obtain  the  coef- 
ficients, Ci  , C2  , C3,  C 4 , and  C,,  The  failure  condition 
for  all  28  tests  are  shown  in  Figure  12.  Also  indicated 
in  this  figure  is  the  regression  line  for  the  failure  con- 
dition. The  results  of  this  analysis  indicated  that  the 
factor  sought,  failure  coverage  (N) , was  of  only  minor 
importance. 
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Table  9 


Predicted  Deflections  of  the  Mat 


at  Failure 

for  Single 

-Wheel  Tests 

Failure  Deflection 

Sect. 

Lane 

Item  Coverage  in. 

2 

3 

1 

600  2.32 

<•> 

3 

2 

120  1.52 

6 

11A 

1 

130  2.68 

6 

11A 

2 

130  1.74 

9 

21 

2 

300  1.70 

10 

2 3 h 

1 

32  2.69 

10 

2 3A 

2 

2 2.50 

10 

23B 

1 

4 4.24 

10 

23B 

2 

2 4.28 

| 

Failure 

Deflection 

j 

Test 

Coverage 

in. 

^ ■ 

' 7 

94 

1.14 

it  ' 

8 

6 

1.74 

j; 

'1 

o 

> 

94 

1.37 

l|  , 

36 

14 

1.80 

J 

i- 

37 

98 

1.56 

38 

197 

1.35 

i 

i 

39 

60 

1.59 

f 

> 

40 

30 

1.49 

i 

41 

30 

1.63 

l 

i 

t 

71 

160 

1.23 

i • 

72 

675 

0.95 

; 

73 

473 

0.97 

' ;■ 

74 

74 

1.59 

! 

75 

360 

1.16 

* \ 

76 

223 

1.43 

i 

91 

570 

0.77 

j, 

93 

23 

1.75 

? 

108 

160 

1.21 

j 

v.j" 

109 

50 

1.38 

t 

{ 

■ j 

*■  ( . 

ci 

>: 

j 
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FIGURE  10.  DEFLECTION  AT  CENTER  OF  WHEEL  FROM 

MODEL  AT  THE  FAILURE  COVERAGE. 


Table  10 


Predicted  Deflections  of  the  Mat 
at  Failure  for  Dual-v;hee.l  Tests 


Sect . 


Lane  Item 


Failure 

Coverage 


Deflection 

in. 


1 

1 

2 

2 

3 

3 

3 

3 

4 
4 
6 
6 
9 
9 

13 

13 

13 

13 


1 

2 

4 

4 

5 

5 

6 
6 
8 
'8 
12 
12 
22 
22 
28 
28 
29 
29 


3 

V 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 

1 

2 


300. 

1.52 

40. 

1.68 

20. 

3.76 

20. 

1.87 

28. 

3.03 

28. 

1.65 

130. 

2.04 

76. 

1.27 

460. 

0.57 

142. 

0.50 

90. 

1.37 

44. 

0.82 

400. 

3.14 

100. 

2.62 

700. 

2.37 

700. 

0.87 

140. 

3.27 

200. 

1.68 
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TIRE  PRESSURE 
0 50  PSI. 


(EIr)Cz(CAr)C’(C8R)C4 

FIGURE  12.  REGRESSION  ANALYSIS  FOR  SINGLE  - WHEEL  TESTS. 


From  the  prototype  data,  it  was  observed  that  the  sane 
tire  ply  rating  was  not  used  in  all  tests;  thus,  the  ratio 
of  tire  inflation  pressure  to  contact  pressure  was  not  a 
constant.  This  observation  indicated  that  possibly  the  tire 
inflation  pressure  might  indeed  be  a significant  parameter. 
Thus  another  regression  analysis  was  performed  on  the  data 
where  each  test  '..-as  expressed  as 


logA  f=C  j log  (*-;Lr)  +C2  log  (EIr)  +C3  log  (CAr)  +q,  log  (CBR) 


4C5log(N)+C6log(TPR) 


(28) 


where  TP  is  the  tire  pressure  normalized  relative  to  100 
psi.  After  obtaining  the  coefficients,  the  failure  condi- 
tions were  as  shown  in  Figure  13.  This  analysis  also  indi- 
cated that  the  failure  coverage  was  of  minor  importance. 

The  dual-wheel  test  sections  were  investigated  by  a 
similar  regression  anailysis  with  an  additional  parameter 
of  wheel  spacing.  Each  dual-wheel  data  set  was  expressed  as 


logAf  =c  J log  (WLr)  +c2  log  (EIR)  +c3  log  (CAR)  +<^  log  (CBR 

+C5log (N) +C 6log (TPR) +C7log (WSR) 


(29) 


where  WSR  is  the  wheel  spacing  normalized  relative  to  25  in. 
Surprisingly,  the  N parameter  was  again  found  not  to  be  a 
significant  parameter.  In  view  of  the  fact  that  the  N 
parameter  is  of  primary  importance  to  the  establishment  of 
a failure  criterion,  the  regression  analysis  program,  as 
applied  to  the  data  in  the  form  of  Ecn.  (29) , was  also 
discarded. 


Educated  Trial  and  Error  Procedure 

Finally,  an  educated  trial  and  error  procedure  was 
undertaken  to  achieve  the  combination  of  test  parameters 
whose  order  was  more  dependent  upon  the  number  of  coverages, 
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N.  Combinations  of  parameters  vcre  set  up  by  deleting 
certain  para:.. ■-'tors  and  various  forms  of  the  function  '/ere 
obtain  ad  by  fystin^  exponents.  For  each  combination, 
corrpuoatJ  or.c  ’’ore  oirectod  toward  the  development  of  a 
ldr.'C’  - re 2 " t:  orshi„>,  on  ~ la /-log  plot,  between  the  pre- 
dicted' fnil’r-e  refections  (Tables  9 end  10),  the  failure 
covt  rag  , me  the  selected  test  pa  rameters . ( Following 
thi:  pro'  . u ‘a  a r'  *’f  nna'-le  correlation  was  jdvent  rally 
achieve,  for  the  dual-wheel  failure  conditions. 

A linear  approximation  of  the  failure  conditions  was 
established  as  the  failure  criterion  for  both  the  single- 
wheel and  the  dual-. .'heel  tests.  A computer  program  was 
developed  whereby  the  performance  of  a mat  could  be  pre- 
dicted from  the  parameters  of  the  test  section.  This  com- 
puter program  is  given  in  Appendix  G.  It  requires  entries 
of:  modulus  of  elasticity  of  the  mat,  moment  of  inertia  of 

the  mat  per  foot  of  width,  the  length  of  a rectangle  whose 
area  is  equivalent  to  that  of  the  tire  print  and  whose  width 
is  the  same  as  the  actual  tire  print,  the  number  of  wheel 
loads,  the  weight  of  water  per  unit  volume  of  soil,  CBR, 
wheel  load,  contact  area,  tire  inflation  pressure,  the  dis- 
tances to  the  beginning  and  end  of  each  wheel  load  from  an 
arbitrarily  located  origin  and  the  magnitude  of  the  wheel 
load  expressed  as  a uniform  load  over  the  width  of  the  tire 
print. 


The  equation  of  the  failure  criterion  line  can  be 
expressed  as 


Af  = C(N)M  (30) 

where  A is  either  the  deflection  of  the  center  of  the  wheel 
for  the  single-wheel  test  or  the  deflection  at  the  center- 
line  of  the  assembly  for  the  dual-wheel  test  and  C and  M are 
each  constants  which  are  dependent  on  the  type  of  test.  In 
the  computer  program,  the  N parameter  is  incremented  and  the 
deflection  at  the  specified  point  is  computed.  For  values 
of  N less  than  200  the  increment  interval  is  5 and  for 
values  of  N greater  than  200  the  increment  interval  is  20. 

For  a single-wheel  test  as  the  value  of  N increases, 
the  deflection  as  given  by  the  failure  criterion  will  also 
increase;  thus,  the  deflection  will  follow  the  line  "C”  in 
Figure  14.  From  the  simulation  model,  she  deflection  of  the 
same  point  also  increases  as  the  value  of  N increases,  lino 
"D",  Figure  14.  As  N is  incremented  the  model  deflection  is 
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(LOG  SCALE) 


FIGURE  14.  PERFORMANCE  TEST  FOR  A 
SINGLE  WHEEL  TEST 


(LOG  SCALE) 

FIGURE  15.  PERFORMANCE  TEST  FOR  A DUAL 
WHEEL  TEST. 


computed  and  compared  to  the  deflection  required  by  the 
failure  criterion,  '..hen  the  latter  deflection  is  just  equal 
to  or  greater"  than  that  of  the  model,  the  corresponding 
value  of  i;  is  considered  the  failure  coverage. 

The  procedure  given  above  was  also  followed  for  the 
dual-'-'heei  tesrs.  For  this  case  the  failure  criterion 
indicate:'  th.ee  oho  deflections  cecrcased  as  N increased, 
line  "F",  Figure  15.  Thus,  for  small  values  of  N,  the 
deflection  as  given  by  the  criterion  would  be  larger  than 
the  model  deflections,  line  "G",  Figure  15.  Failure  was 
considered  to  be  imminent  when  the  model  deflection  became 
just  equal  io  or  greater  than  the  deflection  given  by  the 
criterion. 


Results 


VThen  the  trial  and  error  procedure,  which  was  biased 
toward  a significant  N parameter,  was  applied  to  the  dual- 
wheel test  data  the  "best"  ordering  of  the  parameters 
resulted  in  a plot  of  the  failure  conditions  as  shown  in 
Figure  16.  All  data  points,  except  three  shown  encircled, 
were  fitted  with  a linear  function  using  a least  squares 
technique.  The  three  points  were  believed  to  be  in  error 
because . 

1)  The  maximum  differential  deformation  for  Section  1, 
Lane  2,  Item  3 at  the  indicated  failure  coverage 
was  2.31  in.  This  was  less  than  the  established 
criterion  for  failure.  To  meet  the  roughness 
criterion  more  actual  coverages  would  be  necessary. 
This  indicated  that  the  failure  condition  should 
place  the  point  to  the  left  of  that  shown  in 
Figure  16. 

2)  The  deflection  pattern  for  Section  4,  Lane  8, 

Item  1 was  pronouncedly  unsymmetrical.  This 
indicated  the  likelihood  that  the  reported  soil 
properties  were  not  representative  of  the  soil 
conditions  that  governed  in  the  process.  If 
this  was  the  case,  then  the  resulting  predicted 
deflection  should  be  larger  than  shown  in  Figure 
16*  I*'  sedition,  the  maximum  differential 

deformation  for  this  section  at  failure  was 
recorded  as  2.63  in.  Consideration  of  these 
two  factors  would  place  the  failure  condition 

up  and  to  the  left  of  that  sho\/n  in  Figure  16. 
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3) 


The  CBR  values  reported  at  the  beginning  of  test, 
Section  4,  Lane  C , Item  2,  shoved  a decrease  of 
strength  with  depth  for  the  upper  lc  in.  of  soil. 
For^all  other  cual-vheel  test  sections  the  irr  erse 
of  this  condition  existed.  It  appears  that  tre 
performance  of  this  test  section  is  governed  by 
weaker  underlying  material.  Accounting  for  the 
weaker  conditions,  the  predicted  deflection  at 
failure  would  have  to  be  larger  than  that  shown 
in  Figure  16. 

The  "calculated"  trial  and  error  procedure  was  also 
applied  to  the  single— vhee3  test  data.  The  failure  condi- 
tions for  one  of  the  "better"  parametr'c  combinations  of 
these  is  shown  in  Figure  17.  In  an  attempt  to  establish  a 
failure  criterion  for  the  single-wheel  tests,  the  data 
shown  were  fitted  by  a linear  function.  The  three  points 
indicated  by  circles  in  Figure  17  were  omitted  in  the  curve 
fitting  procedure.  These  points  represent  tests  which 
failed  at  four  or  less  coverages.  In  view  of  the  extremely 
low  number  of  coverages  to  failure,  it  was  felt  that  these 
tests  should  not  be  included  as  they  did  not  represent 
reasonable  engineering  solutions  for  which  mats  could  be 
employed. 

The  computer  program,  Appendix  G,  developed  to  predict 
the  performance  of  a test  section,  incorporates  the  linear 
relationships  indicated  in  Figures  16  and  17.  The  data  for 
the  given  prototype  tests  [ 10,  12]  were  processed  by  this 
computer  program.  The  predicted  failure  coverages,  along 
with  the  actual  failure  coverages,  for  the  dual-wheel  tests 
are  given  in  Table  11.  Similar  tabulation  for  the  single- 
wheel tests  is  provided  in  Table  12. 


Summer-,  and  Discussion  of  Fai  lure  Gy  iter <_ 

The  roughness  criterion  established  by  *•;,£  Corps  of 
Engineers  [ 12  1 for  the  prototype  tests  was  based  upon  the 
permanent  displacements  of  the  unloaded  mat.  With  the 
developed  elastic  mat-soil  model  displacements  of  the  mat 
could  be  determined  only  when  the  surface  was  subjected  to 
external  loading.  Thus  in  oraer  to  use  the  developed  mat- 
soil  model  as  a means  of  predicting  performance,  the 
behavior  of  the  model,  under  actual  or  fictitious  loading, 
had  to  be  related  to  failure.  Initially,  a fictitious  load 
with  a distribution  identical  to  the  traffic  distribution 
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DEFLECTION  AT  t WHEEL*  IN. 


SINGLE -WHEEL  TESTS. 


Table  11 


Comparison  of  Performance 
for  Dual-Uheel  Tests 


Actual  Failure 

Predicted  Failure 

Section 

Coverage 

Coverage 

1-1-3 

300 

220 

1-2-3 

40  (56)* 

260 

2-4-1 

20 

25 

2-4-2 

20 

25 

3-5-1 

28 

35 

3-5-2 

28 

35 

3-6-1 

130 

150 

3-6-2 

76 

85 

6-12-1 

90 

120 

6-12-2 

44  (54)* 

65 

9-22-1 

400 

320 

9-22-2 

100  (56)* 

30 

13-28-1 

700  (967)* 

880 

13-2C-2 

700 

1560 

13-29-1 

140 

170 

13-29-2 

200  (145)* 

130 

4-8-1 

460 

> 5000 

4-6-2 

142 

1000 

* More  realistic  values,  see  Summary  and  Discussion  of 
Failure  Criterion. 
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Table  12 

i 

Comparison  of  Performance 

for  Single-Wheel  Tests 

Test 

Actual  Failure 

Predicted  Failure 

Coverages 

Coverages  * 

2-3-1 

600 

680 

2-3-2 

120 

120 

6-llA-l 

130 

125 

6-11A-2 

130 

115 

: 9-21-1 

* 600 

400 

9-21-2 

300 

135 

7 

94 

40 

8 

6 

65  ' 

9 

94 

105 

36" 

12-16 

75 

37 

70-126 

140 

38 

170-224 

120 

39 

40-80 

120 

40 

20-40 

75 

: 41 

20-40 

. 135 

71 

150-170 

45 

! 72 

650-700 

80 

\ 73 

400-546 

40 

r 74 

40-108 

50 

; 75 

350-370 

110 

! 76 

170-276 

55 

91 

570 

35 

92 

300 

60 

:i  93 

23 

220 

::  108 

160 

195 

* 109 

50 

170 

V 10-2 3A-1 

32 

580 

10-23A-2 

2 

340 

> 10-23B-1 

4 

> 5000 

10-23B-2 

2 

4940 
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was  assumed  to  act  on  the  mat-soil  model.  The  magnitude  of 
this  loa^  was  determined  from  simulation  of  the  maximum 
deform. cion.  Investigation  of  the  curvature  of  the  result- 
ing \ of  el  v- ailection  pattern  revealed  that  the  roughness 
criterion  could  not  be  satisfied  at  the  defined  failure 
coverage.  Comparisons  of  the  magnituae  of  the  fictitious 
lo'-cs  to  a.iiure  coverages  was  also  unsuccessful. 


An  attempt  was  made  to  develop  a functional  relation- 
ship between  raj  lure  coverages  and  test  parameters  by  dimen- 
sional moddlin^.  The  results,  Table  F,  indicated  that  a 
unicue  relationship  could  not  be  established.  The  inability 
to  develop  a satis  factory  functional  relationship  by  this 
procedure  was  attributed  to:  1)  consideration  of  only  a 

limited  number  of  terms  in  the  series  expansion,  and/or  2) 
omission  of  significant  factors  in  the  selected  dimensional 
model.  For  this  study  wore  terms  were  not  considered  due  to 
the  lack  of  reliable  data  available.  For  example,  to  include 
the  next  term  in  the  series  for  the  single-wheel  tests  would 
require  fifteen  additional  sets  of  data.  Further  examination 
of  Eon.  (25)  failed  to  establish  which  tr  term  was  most  influ- 
ential in  the  functional  relationship. 


The  mat-soil  model  was  subjected  to  loading  identical 
to  that  used  to  obtain  the  "average  deflection"  patterns 
[12].  Under  this  loading  the  mat  deflection  was  computed 
for  each  cp__ific  failure  coverage.  For  consistency,  the 
mat  deflection  was  computed  at  the  center  of  the  wheel  and 
under  the  centerline  of  the  wheel  assembly  respectively  for 
the  single-wheel  ana  dual-wheel  tests.  These  deflections 


were  then  expressed  as  a product  function  of  the  failure 
coverage  and  remaining  test  parameters.  A regression  analy- 
sis of  both  the  single— wheel  and  dual— wheel  test  data  indi- 
cated that  the  failure  coverage  was  an  insignificant  param- 
eter, Figures  12  and  13. 


An  educated . trial  and  error  procedure  was  undertaken  to 
achieve  an  ordering  of  the  test  parameters  wnich  was  depen- 
dent upon  the  failure  coverage.  Numerous  combinations  of 
the  test  parameters  were  investigated.  For  each  combination, 
computations  were  Directed  toward  the  development  of  a linear 
relationship  on  a log— log  plot.  A reasonable  correlation  of 
the  c.ual— wneel  tests  was  achieved  and  the  results  are  pre- 
sented in  Figure  16.  The  results  as  shown  were  fitted* with 
a linear  function  which  was  established  as  the  failure  cri- 
terion. Using  this  criterion  the  performance  of  each  test 
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section  was  predicted  v:ith  the  aid  of  the  computer  program, 
Ajjpesu.;  ix  G.  iw  waits  front  this  program  may  be  compared  with 
th ' actual  ini  lure-  coverage,  Table  11.  The  results  indicate 
that  the  porf ormanca  of  similar  tests  may  be  predicted  with 
reasonable,  confidence  by  this  procedure.  Apparent  discrepan- 
cies for  specific  tests  may  be  attributed  to  the  following  * 
ebso eve tioar  of  the  prototype  data: 


1) 


2) 


3) 


4) 


5) 


Section  1,  Lane  2,  Item  3.  This  test  section  was 
considered  to  have  failed  at  40  coverages  due  to 
roughness  and  mat  deterioration.  However,  the 
maximum  differential  deformation  at  this  coverage 
level  was  2.31  inches.  To  meet  the  established 
failure  criterion  for  roughness,  additional  cover- 
ages would  be  necessary.  Thus  it  is  believed  that 
the  performance  of  this  test  was  governed  by  mat 
deterioration  rather  than  by  surface  roughness 
upon  which  the  prediction  of  performance  is  based. 

Section  6,  Lane  12,  Item  1.  The  "average  deflec- 
tion" patterns  for  the  various  coverages  when  the 
load  assembly  was  at  the  center  of  the  panel  are 
pronouncedly  unsymmetrical.  In  addition  the  deflec- 
tion did  not  increase  consistently  with  the  number 
of  coverages.  Due  to  the  unsymmetrical  character- 
istic of  the  deflection  patterns,  the  soil  proper- 
ties reported  may  not  be  representative  of  the 
governing  subgrade  conditions. 

Section  6,  Lane  12,  Item  2.  The  maximum  differential 
deformation  at  failure  was  reported  as  2.63  in.  Con- 
formance with  the  roughness  criterion  would  necessi- 
tate more  coverages  and  therefore  better  agreement 
with  the  predicted  failure. 

Section  9,  Lane  22,  Item  2.  The  maximum  differential 
deformation  at  failure  was  reported  as  4.39  in.  This 
indicated  that  the  failure  criterion  was  satisfied 
at  a coverage  level  less  than  100.  Applying  a linear 
interpolation  between  the  maximum  differential 
settlement  reported  for  20  coverages  and  100 
coverages,  the  section  would  have  satisfied  the 
failure  criterion  at  56  coverages. 

Section  13,  Lane  29,  Item  2.  The  maximum  differ- 
ential deformation  at  failure  was  4.00  in.  Again 
applying  a linear  interpolation  for  the  maximum 
differential. deformations  reported  for  early 
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coverage  levels,  failure  would  have  occurred  at 
145  coverages. 

6)  Section  4,  Lane  S,  Item  1.  As  noted  earlier  the 
deflection  pattern  was  pronouncedly  unsymmetrical. 
The  soil  properties  reported,  which  were  usee,  zor 
the  nat-scil  moo  el,  may  represent  a stronger  sub- 
grade  than  actually  existed.  If  there  was  an 
appreciable  oiffcrenca  between  the  worst  soil  con- 
ditions a no  the  conditions  as  reported,  one  should 
not  expect  the  performance  of  this  section  to  be 
predicteo  by  the  developed  procedure. 

7)  Section  4,  Lane  8,  Item  2.  Also,  as  noted  earlier, 
the  strength  of  the  subgrade  as  reflected  by  the 
CBR  values  decreased  with  depth.  The  performance 
of  this  section  appeared  to  have  been  governed  by 
underlying  weaker  material. 

8)  Section  13,  Lane  28,  Item  1.  This  section  failed 
at  700  coverages  due  to  elastic  deflections;  that 
is,  excessive  displacements  were  in  evidence  under 
che  wheel  loads.  This  failure,  therefore,  could 
not  be  attributed  to  roughness.  Extrapolation  of 
the  given  deformation  data  indicated  that  the  sec- 
tion would  satisfy  the  roughness  criterion  at  S67 
coverages . 

9)  Section  13,  Lane  28,  Item  2.  No  apparent  reason 
can  be  found  from  the  available  data  for  the  dis- 
crepancy shown  for  this  section  between  observed 
and  predicted  performance. 

Apparently  as  incicatea  in  Table  12  and  again  in  Figure 
17,  the  developed  procedure  demonstrated  less  success  ’■’hen 
predicting  the  performance  of  single-wheel  tests.  The 
paucity  of  available  information  and  test  results  for  the 
sections  showing  the  poorest  comparison  of  performance 
(sections  reported  in  reference  [ 10])  rencer  the  correlation 
for  the  single-wheel  tests  academic.  This  is  particularly 
true  since  the  three  inch  roughness  criterion  was  not  specif- 
ically employed;  further,  the  magnitudes  of  the  differential 
deformations  and  the  "deformation"  patterns  were  also  lack- 
ing. If  one  can  speculate,  a prime  factor  which  might 
effect  this  correlation  is  the  variability  of  the  number  of 
passes  necessary  to  complete  one  coverage.  For  the  tests 
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erteo  in  reference  ClO  ) these  ranged  from  a minimum  of  7 
* nioni'mum  of  29  passes.  Thus  the  input  energy  level  per 
''  coulo  be  rmito  different  for  seemingly  similar 
t ion:  . bonce,  it  was  thought  advisable  to  turn  to  an 
covcc  nodel  whereby  actual  loadings  could  be  simulated. 
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ASYMMETRIC  LOAD  MODEL 


Introduction 


The  "deformation"  patterns  of  the  prototype  landing 
mat  tests  conducted  by  the  Corps  of  Engineers  [ 1 2)  indicate 
that  the  time  rate  of  deformation  decreases  with  coverages. 
Typical  mat  deformations  at  the  centerline  of  the  traffic 
lane  when  plotted  as  a function  of  coverages,  shown  by  the 
circles  in  Figure  1C,  suggest  that  a continuous  plot  of  this 
relationship  might  follow  the  dash-line.  Larew  and  Leonards 
[ 33]  obtained  the  same  general  shape  in  their  study  of  the 
deformation  characteristics  of  fine-grained  soils  subject  to 
repetitive  loads.  They  concluded  that  for  repetitive  loads, 
wherein  the  ratio  of  the  repeatedly  applied  deviator  stress 
to  the  static  deviator  stress  causing  shear  failure  was  in 
excess  of  0.91,  the  behavior  might  best  be  represented  by  a 
viscoelastic  model. 

Most  works  to  date  which  employ  viscoelastic  models 
assume  for  mathematical  expediency,  the  existence  of  a con- 
tinuously applied  constant  load  or  a single  step  load  in 
time  [20,  55].  Generally,  the  first  type  of  load  is  utilized 
in  the  study  of  creep  phenomena  while  the  latter  is  employed 
in  the  study  of  stress  relaxation. 

The  time-dependent  response  of  a soil  may  assume  a 
variety  of  forms  depending  on  such  factors  as  soil  type, 
soil  structure,  stress  history,  type  of  loading,  and  other 
factors  [58  ].  In  most  studies  these  factors  are  classified 
into  two  croups  of  effects:  a creep  effect  and  a relaxation 

effect.  It  has  been  shown,  Singh  and  Mitchell  [58]  and 
Konaer  and  Krizek  [31]  , that  in  general  soils  exhibit-both 
linear  and  non— linear  behavior.  It  is  also  commonly  assumed 
that  the  magnitude  of.  the  strain  is  stress  dependent.  Thus 
any  strain— time  function  which  is  to  describe  the  response 
of  such  materials  must:  1)  be  applicable  to  the  range  of 

stress  levels  encountered;  2)  account  for  both  linear  and 
non-linear  relationships  between  strain  and  time;  and  3)  be 
relatively  simple  in  form  for  ease  of  application  and  deter- 
mination of  parameters. 
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COVERAGES 

FIGURE  18.  TYPICAL  DEFORMATION  AT  THE  CENTERLINE  OF 

TRAFFIC  LANE. 


In  this  section  a viscoelastic  model  is  developed  for 
predicting  the  response  of  landing  mats  to  a simulation  of 
the  actual  sequence  of  traffic  loads.  Although  only  the 
subgrecc  will  he  oaken  to  possess  viscoelastic  properties, 
tie  rocc-1  oarav'.eters  which  are  established  from  the  response 
of  actual  tests  will  also  indirectly  reflect  the  viscoelastic 
behavior  of  the  lanaing  mats. 


/ nalysis  of  Model  Behavior 

Assume  that  a load  applied  to  a mat-soil  system  pro- 
duces a constant  stress  Oj  , as  shown  in  Figure  19(a).  The 
response  of  the  system  due  to  this  stress  will  be  as  shown 
in  Figure  19(b)  wherein  e represents  the  immediate  elastic 
strain  and  e(t)  represents  the  strain  for  times  t > 0.  If 
transient  effects  are  considered  negligible  in  the  elastic 
range,  e will  be  a constant  and  e (t)  may  be  expressed  as 


c(t)  = cQ  + Ae  (t) 


(31) 


where  Ac(t)  contains  all  the  time  effects,  Figure  19(b).  An 
increase  in  the  magnitude  of  the  constant  stress  will  simply 
displace  the  response  curve  upward,  Figure  19(c).  Assuming 
strain  is  time  dependent  and  stress  is  independent  of  time, 
r.r.  e.eeccirtcd  ctrccc-etroin  rclcticrchip  will  take  the  form 


e(t)  = Du(t)o  (32) 

where  cr  represents  the  magnitude  of  the  (constant)  stress 
and  D (t)  is  a "compliance"  function.  According  to  Schapery 
[55]  this  function  can  be  expressed  as 


Du(t)  = Do  + D(t)  (33) 

where  D is  a test  dependent  constant 

D?t)  is  a function  of  time  so  defined  that  D(0)  = 0. 

The  assumption  is  commonly  made  [55,  58  ] that  time  effects 
can  be  expressed  by  a power  law  of  the  form 

D(t)  = Di^~ j (0<  n < 1)  (34) 
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TIME  (t) 

o)  APPLIED  STRESS 


b)  RESPONSE  DUE  TO  Oi 


TIME  (t) 

c)  EFFECT  OF  STRESS  LEVEL 


FIGURE  19.  STRESS  - STRAIN  - TIME  BEHAVIOR. 
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where  Dx  and  n are  parameters  to  be  determined  from  response 
curves;  an  is  the  shift  factor  which  is  stress  dependent  and 
oxpressed~'as 

Vn 


and  o.  is  a constant.  For  constant  stress  conditions  a is 
seen  £0  be  equal  to  unity.  Hence  from  Eqns.  (32)  and  (-53), 
the  strain  at  any  time  t > 0 under  a constant  stress  will  be 
given  by 

€(t)  = (D0  + q tn)0  (36) 

If  a step  stress  of  time  duration  t*  is  applied,  the 
corresponding  strain  will  be  as  shown  in  Figure  20.  For  any 
time  t < 1 1 the  strain  will  be  defined  by  Eqn.  (36).  For 
times  t > t the  residual  strain  in  the  system  can  be 
expressed  b$ 


aD  = aD(o)  - 


0/0. 
sinh ,/ 


®r(t)  = b0  + Dj  (t)n  10  - l Dq  + Dj  (t  - tx  )% 

or  (37) 

er(t)  - [ D,  tn  - Dj(t  - tj  )r*]0 

Development  of  Landing  Mat-Soil  Model 

Under  the  action  of  a.  single  pass  of  a load  vehicle,  a 
point  on  the  surface  of  the  subgrade  experiences  a normal 
stress  variation  similar  to  that  shown  in  Figure  21.  This 
stress  variation  can  be  represented  bv  a step  stress  provided 
that  the  time  duration  for  stress  builc-up  and  release  i= 
small.  The  effective  time  duration  of  the  step  stress,  say 
tj  = t.  - t^,  will  be  a fraction  of  the  time  required  to 
complete  one  pass. 

It  will  be  assumed  that  trafficking  of  the  section  was 
continuous  and  performed  at  uniform  speed.  The  actual  time 
to  complete  one  pass  will  then  be  a constant  and  one  ’~nss 
may  serve  as  a basic  time  unit.  Also,  the  duration  tj  of 
the  step  stress  will  be  a constant  independent  of  past . The 
magnituce  of  the  step  stress  will  be  cependent  upon  the 
assembly  load  and  the  position  of  the  point  under  considera- 
tion relative  to  the  path  followed  during  the  pass. 
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a)  APPLIED  STEP  STRESS 
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TIME  (t) 

b)  RESPONSE  DUE  TO  STEP  STRESS 
FIGURE  20.  BEHAVIOR  OF  SYSTEM  Ti 
STEP  STRESS. 
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STRESS 


FIGURE  21.  LOAD  DISTRIBUTION  AT  A POINT  DURING  A 
SINGLE  PASS. 
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b)  MAXIMUM  STRESS  AT  POINT  I 

FIGURE  22.  VARIATION  OF  STRESS  WITH  LOAD 
POSITION . 


Assuming  the  surface  displacement  of  the  subgrade  for 
time  t < tx  can  be  expressed  by  Ecn.  (36),  the  residual 
displacement  at  the  end  of  one  pass  vd.ll  be  given  by 


&r(l)  = [D,(l)n  - D^l-t^o  (1>  ^ ) (38) 

where  tj  is  the  time  duration  of  the  step  stress 
a is  the  magnitude  of  the  step  stress 
and  n are  constants. 

The  superposition  principle  nay  be  used  to  extend  the 
applicability  of  Ecn.  (3o)  to  more  than  one  pass  of  the  load 
vehicle.  For  examp_u,  if  two  passes  of  the  load  vehicle  are 
made  along  the  same  path,  the  residual  deformation  at  the 
end  of  the  second  pass  will  be 

*.( 2)  = Id,  (2)n  - Dj  (2-t j ) n + Dj  (l)n  - D (1-d  )n)a  (39) 

Generalizing  Eqn.  (39)  for  any  number  of  passes  N,  along  the 
same  path,  one  obtains 


m=N 

A (N)  = Dj  or  ( mn  - (rn-tjrl  (40) 

r m=  l 

Under  normal  traffic  conditions  the  path  of  the  load 
vehicle  changes  with  passes.  Thus,  the  magnitude  of  the 
step  stress  is  also  altered  as  the  path  of  the  load  vehicle 
changes,  see  Figure  22.  Therefore  the  stress  is  space 
dependent  and  must  remain  in  the  summation  operation,  Eqn. 
(40)  . The  residual  deformation  at  any  nu:nber  of  passes 
regardless  of  the  path  followed  would  then  be  given  by 

m=N 

A (N)  = Dj  l (mn  - (m-tj  ) ]<j  m (41)* 

m=  l 

where  o is  the  magnitude  of  the  stress  at  the  point  for 
the  specific  position  of  the  load  vehicle. 


* It  is  noted  that  Dj  ^ K [rnn-  (m-t^  )n]  is  analogous  to  the 

m=  1 

reciprocal  of  the  sufcgrade  modulus  used  in  the  Winkler 
hypothesis  assuming  it  to  be  time  dependent. 
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A quasi-elastic  analysis  was  developed  to  predict  the 
residual  deformation  of  a landing  mat-subgrade  system  (com- 
pute r program  given  in  Appendix  H) . The  mocel  consisted  ot 
an  infinite  beam  with  rigidity  identical  to  that  of  the  mat 
and  continuously  supported  by  a subgraoe  whose  response  to 
loading  could  be  expressed  by  Ecn.  (41).  The  loading  on 
the  beam  was  c csignco  to  approximate  the  loading  seouence 
On  ployed  in  the  procotype  rests  (from  reference  l 12J  ) where 
uniform  distribution  of  traffic  was  applied.  The  loading 
was  made  Loth  time  and  space  dependent.  The  sequence  of 
loading  for  a single-wheel  test  is  given  in  Figure  23  and 
for  a dual-wheel  test  in  Figure  24. 

To  be  consistent  with  the  prototype  tests,  the  basic 
time  unit  was  taken  as  one  coverage.  Due  to  lack  of  infor- 
mation as  to  the  actual  time  to  conduct  the  prototype  tests, 
the  duration  of  load  application  was  assumed  equal  to  one- 
half  the  time  recruired  to  make  one  pass.  This  assumption 
may  be  unrealistic  but  lacking  more  information  it  was 
expedient  and  in  addition  it  does  not  effect  the  validity 
of  the  model  or  the  method  of  application.  Only  the  numeri- 
cal values  of  the  parameters  Dj  and  n are  influenced  by  this 
assumption.  Not  having  adequate  test  information  available, 
these  parameters  could  only  be  roughly  approximated. 


The  developed  model  was  applied  t-o  a dual  wheel  test. 
The  values  of  the  parameters , Dj  and  n,  were  obtained  by 
comparing  the  model  response  for  a point  located  on  the 
centerline  ot.  s-he  traffic  lane  to  the  actual  response  of 
the  same  poi  t luring  the  prototype  test.  In  Figure  25, 
the  fitted  mocai  deformation  of  a point  on  the  centerline 
is  indicate  by  the  curve  labeled  "theoretical1'  and  the 
actual  deformation  is  given  (by  crosses)  at  three  coverage 
levels.  The  actual  deformations  shown  are  the  average  of 
the  centerline  deformations  at  three  coverage  levels 
obtained  from  the  two  "deformation"  patterns  given  in 
reference  [12].  Some  predicted  deformation  patterns  at 
various  coverage  levels  are  shown  in  Figure  26. 

The  model  was  also  employed  for  a single  wheel  test. 
Some  results  are  given  in  Figure  27  and  Figure  26. 
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NOTE;  NUMBER  ON  CURVE  • NUMBER  OF  COVERAGES 

FIGURE  26.  PREDICTED  DEFORMATION  PATTERNS 
FOR  SECTION  I,  LANE  I,  ITEM  3. 
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NOTE : NUMBER  ON  CURVE  « NUMBER  OF  COVERAGES 

FIGURE  27.  PREDICTED  DEFORMATION  PATTERNS 
FOR  SECTION  2,  LANE  3,  ITEM  I. 


ITEM 


FIGURE  28.  COMPARISON  OF  DEFORMATIONS  (section  2,  lane  3,  item  i>. 


Discussion 


It  should  be  recoiled  that  the  parameters,  Dj  and  n 
shown  in  Figures  25  and  2b,  were  selected  so  as  to  provide 
the  simulation  of  the  deformation  at  a single  point.  i:o 
attempt  v e been  m.a-e  to  attach  any  physical  significance 
to  the  .'es  of  these  parameters.  In  addition  to  the  usual 
test  variables,  wheel  load,  number  of  wheels , wheel  spacing, 
contact  area,  mat  rigidity,  and  soil  conditions,  these 
parameters  are  influenced  by  the  width  of  the  traffic  lane 
and  the  number  of  positions  of  the  load  vehicle  required  to 
complete  one  coverage.  More  detailed  information  of  the 
deformation-coverage  response  of  sections  and  loading 
sequences  are  necessary  to  develop  a procedure  of  identify- 
ing these  parameters  with  conditions  of  either  past  or  con- 
templated future  tests. 

It  is  noted  that  the  shape  of  the  deflection  pattern 
for  the  dual-wheel  test,  Figure  26,  is  very  different  than 
that  for  the  single-wheel  test,  Figure  27.  This  is  due 
primarily  to  the  assemblages  used  and  the  magnitude  of  the 
total  load  applied  to  the  subgrade.  Another  contributing 
factor  is  the  test  procedure  itself.  For  the  dual-wheel 
test  six  different  assembly  locations  were  used  to  complete 
one  coverage?  whereas,  for  the  single-wheel  test,  eleven 
assembly  locations  were  necessary  to  complete  one  coverage. 
To  maintain  a correspondence  between  coverages,  it  was 
assumed  that  the  load  for  the  dual-wheel  test  was  applied 
for  a time  1.835  longer  than  that  for  the  single-wheel  test. 
For  these  conditions,  the  effect  of  the  preceding  position 
of  loading  is  more  pronounced  for  the  dual-wheel  test  than 
for  the  single-wheel  test. 

Ecruation  (41)  may  be  generalized  further.  Under  the 
assumption  of  a steady,  continuous  trafficking,  the  value 
of  t j may  be  taken  as  a constant  and  Eqn.  (41)  can  be 
expressed  as  the  difference  equation 


Ei  t pM-  ; >1  (42) 

m=  i L '-j  J m 

Recognizing  that  the  value  of  tj  is  small  in  comparison  to 


the  time  required  to  complete  one  pass, 
brackets  may  be  taken  as  the  derivative 


respect  to  t 
the  integral1  of 
f(n).  Under  an 
ex Dressed  as 


the  terms  in  the 
of  the  f (m)  with 


The  summation  of  this  quantity  approximates 
the  derivative  which  is  the  function  itself, 
average  stress  condition  Eqn.  (42)  may  be 
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arC» 

Hr—  = Vave£<K>  <43> 


If  the  residual  deformation 
at  two  specific  number  of  passes, 
f (:;)  may  be  established  by 


at  a point  k can  be  measured 

K.  and  N . , the  function 
i 3 


Ar(i;.)  _ f(N.) 

" fTiry 


(44) 


Thus  the  function  f (W)  could  be  taken  as  any  form  which 
satisfies  Eqn.  (44) . 


In  summary,  a mechanistic  model  has  been  developed 
which  is  capable  of  simulating  the  actual  loading  sequence 
and  which  also  provides  a means  whereby  residual  deforma- 
tions of  the  mat  surface  may  be  predicted.  The  potentiali- 
ties of  this  model  have  bean  demonstrated  for  a single-wheel 
and  dual-wheel  prototype  test.  To  fruitfully  pursue  the 
capabilities  of  this  model,  more  detailed  test  data  are  neces- 
sary to  identity  or  attach  physical  significance  to  the 
model  parameters.  A continuous  record  of  the  deflection  as 
a function  of  the  number  of  passes  is  necessary  along  with 
a record  of  the  position  of  the  load  vehicle  during  each 
pass.  Also,  information  concerning  the  speed  of  the  load 
vehicle  and  the  time  interval  (expressed  as  a function  of 
the  time  required  to  complete  one  pass)  required  to  shift 
the  position  of  the  load  vehicle  is  necessary.  With  thxS 
additional  data  and  by  selecting  a single  pass  rather  than 
a single  coverage  as  the  basic  time  unit,  better  correspon- 
dence between  deformation  patterns  can  be  achieved  thus 
allowing  the  possible  establishment  of  a limiting  deformation 
as  the  failure  criterion. 
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MOMENT  TRANSFER  INVESTIGATION 


In  most  o£  the  test  sections  reported  in  the  Corps  of 
Engineer's  report  (12}  it  was  observed  that  the  distress  in 
the  traffic  lane  was  quite  pronounced  in  the  vicinity  of 
the  end  joints  of  the  mat  elements.  This  is  why  in  all 
simulations  the  load  wat  laced  in  the  center  of  the  mat 
element.  In  view  of  t'  is  observation,  a study  was  under- 
taken to  investigate  tne  effectiveness  of  the  end  joints 
to  transfer  moment  from  one  mat  element  to  another  under 
simulated  test  conditions. 

Due  to  the  masonry  type  of  placement  of  the  mats, 

Figure  29,  it  was  not  possible  to  determine  the  amount  of 
moment  transferred  solely  across  the  end  joints.  A signifi- 
cant amount  of  moment  could  be  transferred  around  the  end 
joint  between  mat  elements  A and  B,  Figure  29,  due  to  the 
transfer  of  shear  across  the  transverse  joints.  An  inves- 
tigation that  considers  the  presence  of  elements  C and  D 
under  field  conditions  is  further  complicate,  by  the  contact 
length  of  the  tires  under  test  loading  being  greater  than 
the  width  of  the  mat  element.  Thus  it  was  decided  in  this 
study  to  restrict  the  investigation  to  the  overall  influence 
that  the  presence  of  an  end  joint  has  on  the  transverse 
continuity  of  the  traffic  lane.  Emphasis  will  be  placed 
upon  the  performance  of  the  Til  aluminum  and  the  MS  steel 
mats  as  reported  by  the  Corps  of  Engineers ( 10,  12)  . 

The  type  of  end  connector  employed  for  the  Til  mat  was 
described  in  some  detail  in  the  Corps  of  Engineer's  report 
111].  The  method  of  connecting  the  MB  mat,  which  is  a 
self-locking  mat,  was  also  described.  A modified  or 
strengthened  joint  connection  was  given  in  a later  report 
l 9].  An  examination  of  the  performance  of  these  connectors 
suggested  efficient  shear  transfer  across  the  enc  joints; 
however,  the  capacity  of  the  end  joint  to  transfer  moment 
appeared  to  be  less  than  that  of  the  mat  itself.  The  latter 
is  attested  to  by  the  enc  joints  shown  in  the  photographs, 
Figure  22,  Part  III  and  Figure  8,  Part  XI  of  ( 12 ] . As  a 
consequence  of  the  above,  the  assumption  was  mace  that  in 
the  simulation  model  there  exists  complete  transfer  of 
shear  across  the  end  joint. 
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FIGURE  29.  ARRANGEMENT  OF  MAT  ELEMENTS 
WITHIN  THE  TRAFFIC  LANE. 
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FIGURE  30.  SUPERPOSITION  OF  LOAD  CONDITIONS 


Since  the  subgrade  in  the  prototype  tests  was  con- 
structed under  tight  control,  it  was  further  assumed  that 
the  strength  of  the  subgrace,  as  expressed  by  the  k param- 
eter, was  uniform  throughout  the  traffic  lane.  In  line 
with  this  assumption,  any  variation  in  the  "average  deflec- 
tion" patterns  (between  that  at  the  end  joint  and  that  at 
the  center  of  panel)  might  then  be  taken  as  a measure  of 
the  d'.grec  to  which  the  end  joint  provioed  moment  transfer. 

The  mat-soil  model  developed  earlier  was  employed  to 
simulate  the  "average,  deflection"  patterns  at  the  end  joints. 
The  k parameter  was  assigned  the  values,  Table  13,  Column  5, 
which  resulted  in  a minimum  of  the  error  functional  for  the 
previous  simulation  of  the  "average  deflection"  patterns  at 
the  center  of  the  panel.  These  k values  were  used  to  exclude 
any  inherent  errors  in  Eqns.  (13)  and  (14).  The  parameter  y 
was  permitted  to  vary  with  coverages  in  accordance  with 
Eqn.  (15) . 

Postulating  the  validity  of  superposition,  Figure  30, 
two  distinct  steps  were  followed  to  simulate  the  "average 
deflection"  pattern  reported  in  the  vicinity  of  an  end  joint. 
Initially,  the  mat  was  considered  to  be  of  infinite  extent 
and  devoid  of  any  joints.  Actual  wheel  loads  were  applied 
symmetrically  about  what  would  be  the  location  of  an  end 
joint.  The  moment  developed  in  the  mat  at  the  location  of 
the  fict.it end  joint  was  then  computed  in  addition  to 
the  deflections  of  the  nine  or  more  discrete  points  used  to 
characterise  the  "average  deflection"  pattern.  Next,  the 
mat  was  considered  to  be  hinged  and  a concentrated  moment 
was  applied  at  the  end  joint,  Figure  30(c).  The  mathematical 
development  for  this  condition  is  given  in  Appendix  E.  The 
deflections  of  t.he  mat  under  the  action  of  the  concentrated 
moment  loading  were  determined  and  algebraically  added  to 
the  deflections  of  the  mat  wherein  complete  continuity  was 
assumed.  Simulation  of  the  "average  deflection"  pattern 
was  achieved  by  increasing  incrementally  the  magnitude  of 
the  concentrated  moment  until  the  error  functional,  Eqn. 

(D2) , was  a minimum.  The  increment  size  was  taken  as  one 
percent  of  the  moment  which  existed  at  the  ficticious  joint 
in  the  infinite  beam.  The  deficiency  of  moment  transfer 
across  the  end  joint  was  taken  as  the  ratio  of  the  magnitude 
of  the  concentrated  moment  to  the  moment  in  the  assumed 
infinite  beam.  These  deficiency  percentages  are  listed  in 
the  last  column  of  Table  13. 
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Table  13 


End  Joint  Iloment  Transfer  Deficiency  Based 
Upon  Uniform  Subgrade  Conditions 


Sect.  Lane  Item 

1 1 3 

1 2 3 

2 3 1 

2 3 2 

2 4 1 

2 4 2 

3 5 1 

3 5 2 

3 6 1 

3 6 2 

6 11A  1 


Cov. 

k-pci 

Deficiency 

0 

54.1 

1. 

20 

41.4 

1. 

200 

41.4 

1. 

300 

46.0 

1. 

0 

52.5 

337. 

20 

43.5 

158. 

40 

39.8 

197. 

0 

38.4 

107. 

200 

28.8 

187. 

600 

28.7 

271. 

0 

42.6 

531. 

120 

57.4 

2896. 

200 

30.3 

711. 

0 

19.9 

72. 

20 

16.5 

122. 

0 

33.0 

246. 

20 

32.2 

155. 

0 

17.7 

1427. 

30 

13.2 

504. 

0 

23.5 

104. 

30 

19.2 

1. 

0 

12.3 

1. 

20 

10.5 

1. 

76 

15.0 

1. 

156 

13.9 

1. 

0 

21.8 

78. 

20 

15.9 

32. 

76 

18.4 

41. 

0 

32.1 

73. 

20 

27.0 

59. 

130 

28.4 

64. 
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1 

Table 

13  (cont'd. 

) 

1 Sect. 

Lane 

Item 

Cov. 

k-pci 

Deficiency 

JJ  6 

IIP. 

2 

0 

46.8 

1. 

20 

38.0 

68. 

1 

130 

32.6 

47. 

1 6 

12 

1 

0 

21.6 

48. 

20 

IS. 5 

28. 

| 

90 

21.6 

1. 

1 6 

12 

2 

0 

60.0 

134. 

20 

34.1 

53. 

Is 

44 

34.3 

79. 

||  9 

21 

1 

0 

18.9 

85. 

20 

15.6 

16. 

II 

200 

12.6 

8. 

1 

300 

11.3 

16. 

I 9 

21 

2 

0 

29.1 

61. 

20 

26.7 

106. 

i 

200 

25.9 

175. 

i 

300 

21.8 

186. 

1 9 

22 

1 

0 

13.9 

1. 

20 

15.3 

1. 

1 

400 

9.6 

1310. 

' q 

n 3 

f 

22 

2 

0 

16.8 

105. 

20 

16.8 

71. 

i| 

i i 

100 

16.8 

1. 

f! 

r t 

;'!  10 

23A 

1 

0 

42.9 

1. 

32 

32.9 

1376. 

10 

23B 

1 

4 

32.1 

1. 
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For  many  coverage  levels  the  percentage  of  moment 
deficiency  was  found  to  be  greater  than  100;  obviously,  the 
selected  assumptions  were  in  error,  it  vas  apparent  for 
these  cases  that  even  a perfect  hinge  could  not  account 
mathematically  for  variations  in  the  "average  deflection" 
patterns . 


To  examine  the  assumption  of  uniform  support  (as 
expressed  by  a unique  k parameter),  the  "average  deflection" 
patterns  for  the  end  joints  v/ere  again  simulated  by  the 
mat-soil  model  assuming  an  infinite  beam  and  that  no  joints 
existed.  A new  value  of  the  k parameter  was  obtained  from 


this  simulation.  A comparison  of  these  k values,  denoted 
by  k , with  those  obtained  from  the  simulation  of  the 
"average  deflection"  patterns  at  the  center  of  the  mat  ele- 
ment can  be  had  from  Table  14,  Columns  5 and  7.  It  may  be 
observed  that,  in  general,  the  difference  in  the  k values 


is  not  large;  however,  they  are  large  enough  to  influence 
the  magnitude  of  deflections  because  of  the  great  sensi- 
tivity of  the  model  to  this  parameter,  particularly  at 
relatively  low  values  of  k.  The  respective  simulation 
errors  for  the  continuous  mat  analysis  are  also  given  in 
Table  14.  A comparison  of  these  errors  indicates  that 
considerably  better  simulation  of  the  end  joint  "average 

deflection"  patterns  were  achieved  by  using  k . 

© 


Values  of  k , like  those  of  k,  were  found  to  decrease 
with  the  number  of  coverages.  More  importantly,  the  value 
of  k at  any  coverage  level  was  found,  with  only  very  few 
exceptions,  to  be  less  than  that  of  the  corresponding  k. 
This  tends  to  indicate,  as  was  noted  previously  in  Eqn. 

(13) , that  the  mat  and  subgrade  in  the  vicinity  of  the  end 
joints  effectively  experienced  more  coverages  (in  the  sense 
that  more  energy  was  applied) . 


One  might  speculate  on  the  nature  of  the  behavior  at  a 
joint.  The  difference  in  the  values  of  k and  k^  at  zero 
coverages  could  be  attributed,  in  part,  to  the  inherent 
slack  generally  found  in  end  joint  connections  necessary  to 
permit  such  connections  to  be  made  under  field  conditions. 
Also,  high  stress  concentrations  are  undoubtedly  developed 
in  the  connectors  which  could  produce  local  yielding  of  the 
material.  A strong  case  can  be  mace  for  the  latter  point 
by  examining  the  distress  exhibited  by  the  end  joints  in 
the  photographs  accompanying  the  Corps  of  Engineer's  reports 
( 9,  10,  12  ].* 
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Table  14 

Comparison  of  k Values 


Error 

Error 

Sect. 

Lane 

Item 

Cov. 

k-pci 

k-pci 

in./pt. 

in./pt 

1 

1 

3 

0 

54.1 

.103 

48.1 

.072 

20 

41.4 

.150 

37.4 

.085 

200 

41.4 

.340 

32.4 

.133 

300 

46.0 

.231 

40.0 

.151 

1 

2 

3 

0 

52.5 

.157 

47.5 

.114 

20 

43.5 

.266 

37.5 

.184 

40 

39.8 

.440 

31.8 

.287 

2 

3 

1 

0 

38.4 

.158 

33.4 

.137 

200 

28.8 

.338 

20.8 

.204 

600 

28.7 

.539 

15.7 

.196 

2 

3 

2 

0 

42.6 

.214 

41.6 

.214 

120 

57.4 

.342 

37.4 

.245 

200 

30.3 

.507 

24.3 

.466 

2 

4 

J 

0 

19.9 

.200 

19.9 

.200 

20 

16.5 

.591 

14.5 

.525 

2 

4 

2 

0 

33.0 

.244 

29.0 

.164 

20 

32.2 

.552 

30.2 

.543 

3 

5 

1 

0 

17.7 

.640 

13.7 

.227 

30 

13.2 

1.030 

9.2 

.482 

3 

5 

2 

0 

23.5 

.206 

22.5 

.180 

30 

19.2 

.252 

17.2 

.157 

3 

6 

1 

0 

12.3 

.225 

li.3 

.168 

20 

10.5 

.116 

10.5 

.116 

76 

15.0 

.758 

11.0 

.400 

156 

13.9 

.690 

10.9 

.254 

3 

6 

2 

0 

21.8 

.225 

20.8 

.211 

20 

15.9 

.324 

15.9 

.324 

76 

18.4 

.873 

14.4 

.709 
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Table  14  (Cont'd.) 


Error 

Vpci 

Error 

ct . 

Lane 

IteT* 

Cov. 

k-pci 

in . /pt . 

in./pt. 

6 

11A 

l 

0 

20 

130 

32.1 
27.0 
28  .4 

.211 

.210 

.261 

27.1 

24.0 

23.4 

.152 

.181 

.175 

6 

HA 

2 

0 

20 

130 

48.8 

36.0 

32.6 

.167 

.273 

.210 

48.8 

38.0 

32.6 

.167 

.273 

.210 

e 

12 

1 

0 

20 

90 

21.6 

18.5 

21.6 

.427 

.395 

.362 

21.6 

18.5 

17.6 

.427 

.395 

.137 

6 

12 

2 

0 

20 

44 

60.0 

34.1 

34.3 

.216 

.240 

.232 

53.0 

30.1 
34.3 

.202 

.204 

.232 

9 

21 

1 

0 

20 

200 

300 

18.9 

15.6 

12.6 
11.3 

.259 

.085 

.152 

.201 

13.9 

14.6 

12.6 
11.3 

.089 

.059 

.152 

.201 

9 

21 

2 

0 

20 

200 

300 

29.1 

26.7 
25.9 

21.8 

.119 

.223 

.411 

.422 

25.1 

23.7 
18.9 

18.8 

.096 

.211 

.348 

.409 

9 

22 

1 

0 

20 

400 

13.9 

15.3 

9.6 

.439 

1.034 

2.160 

10.9 

8.3 

4.6 

.199 

.290 

.734 

9 

22 

2 

0 

20 

100 

16.8 

16.8 

16.8 

.213 

.190 

.395 

15.8 

15.8 

12.8 

.200 

.172 

.173 

10 

23A 

1 

0 

32 

42.9 

32.9 

.252 

.850 

36.9 

24.9 

.206 

.761 

10 

23B 

.1 

4 

32.1 

.668 

32.1 

.668 

95 


In  general,  ki  decreased  more  rapidly  with - coverages 
thi.n  k (Table  15,  Columns  5 and  6).  This  behavior  is  not 
surprising  sinct.  in  all  test  sections  the  mots  in  the 
vicinity  of  uhe  end  joint  experienced  progressive  deteriora- 
tion uuc  to  o concentration  of  mat  breaks  (rivet  failure, 
cracking,  ana  curling)  [32].  The  weakened  mat  evidently 
losses  its  effectiveness  and  the  underlying  subgrade  is 
subjc.cce.  to  higher  energy  levels.  Thus  in  the  vicinity 
of  tho  end  joint,  the  subgrade  experiences  larger  shear 
deformations  anc  considerably  more  remolding  than  in  other 
locations . 

Regardless  of  the  exact  nature  of  the  behavior  of  the 
end  joints,  a measure  of  their  effectiveness  might  follow 
from  an  examination  of  the  ratios  of  k to  k.  Table  15, 
Column  7,  provides  a listing  of  these  ffatios.  Discounting 
the  differences  in  the  error  functionals,  average  values 
for  these  ratios  are  0.84  and  0.90  for  the  Til  and  M8  mats 
respectively.  If  credence  can  be  placed  upon  the  above 
supposition,  it  appears  that  the  effectiveness  of  the  end 
joint  connections  is  approximately  10  to  16  percent  less 
than  that  of  the  mat  elements.  The  foregoing  results  are 
not  in  conflict  with  the  observations  of  the  Corps  of 
Engineers  noted  earlier  [9].  Indications  suggest  that  some 
increase  in  serviceability  can  be  achieved  by  further 
strengthening  of  the  end  joint  connector. 
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Lane 

1 


Table  15 

Ratio  of  k values 
Item  Cov.  kNAIKT 


^e**/^ke*INT 


0 

1.00 

1.00 

20 

.77 

.78 

200 

.77 

.67 

300 

.85 

.83 

0 

1.00 

1.00 

20 

.82 

.79 

40 

.76 

.67 

0 

1.00 

1.00 

200 

.75 

.62 

600 

.75 

.47 

0 

1.00 

1.00 

120 

1.35 

.90 

200 

.71 

.59 

0 

1.00 

1.00 

20 

.83 

.73 

0 

1.00 

1.00 

20 

.95 

1.07 

0 

1.00 

1.00 

30 

.75 

.67 

0 

1.C0 

1.00 

30 

.82 

.76 

0 

1.00 

1.00 

20 

.85 

.93 

76 

1.22 

.97 

156 

1.33 

.96 

0 

1.00 

1.00 

20 

.73 

.76 

76 

.84 

.69 

(k  )M/kM 
' e N N 


1.00 

.68 


.92 

1.00 

.74 

.79 

.96 

1.00 

.79 
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Table  15  (Cont'd.) 


Sect . 

Lane 

Item 

Cov. 

V^INT 

(ko>,.As 

6 

11A 

1 

0 

1.00 

1.00 

.65 

20 

.84 

.69 

.69 

130 

.89 

.66 

.83 

6 

llA 

2 

0 

1.00 

1.00 

1.00 

20 

.78 

.78 

1.00 

130 

.67 

.67 

1.00 

6 

12 

1 

0 

1.00 

1.00 

1.00 

20 

.86 

.86 

1.00 

90 

1.00 

.82 

.82 

6 

12 

A 

2 

0 

1.00 

1.00 

.89 

20 

.57 

.57 

.89 

44 

.57 

.65 

1.00 

9 

21 

1 

0 

1.00 

1.00 

.74 

20 

.83 

1.05 

.94 

200 

.67 

.91 

1.00 

300 

.60 

.81 

1.00 

9 

21 

2 

0 

1.00 

1.00 

.86 

20 

.92 

.95 

.89 

200 

.89 

.75 

.73 

300 

.75 

.75 

.86 

9 

22 

1 

0 

1.00 

1.C0 

.79 

20 

1.10 

.76 

.54 

400 

.69 

.42 

.48 

9 

22 

2 

0 

1.00 

1.00 

.94 

20 

1.00 

1.00 

.94 

100 

1.00 

.61 

.76 

10 

23A 

1 

0 

1.00 

1.00 

.66 

32 

.77 

.68 

.76 

10 

23B 

1 

4 

1.00 

1.00 

1.00 
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SUGGESTIONS  FOR  FUTHER  RESEARCH 


Based  upon  the  results  of  this  study,  the  following 
suggestions  for  future  study  are  presented: 

1.  A study  should  be  undertaken  to  establish  a basic 

unit  for  defining  a single  coverage:  one  which 

is  not  test  dependent. 

2.  It  is  urged  that  a comprehensive  study  be  conducted 
to  investigate  the  full  capabilities  of  the  asym- 
metric load  model  described  herein.  Although 
laboratory  testing  of  isolated  subgrade  samples 
subjected  to  repetitive  loads  would  be  of  value  for 
such  a study,  a test  series,  smaller  in  size  than 
the  conventional  prototype  tests,  could  provide 
meaningful  data  for  the  determination  of  the  model 
parameters.  Necessary  information  to  establish  the 
time  duration  for  the  applied  load,  could  be  ob- 
tained by  installing  a continuously  monitored  load 
cell  at  the  subgrade  surface  of  a conventional 
prototype  test. 
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SUMMARY  AND  CONCLUSIONS 


Two  mechanistic  models  have  been  developed  to  simulate 
the  behavior  of  landing  mat  systems.  The  first  model#  which 
is  based  upon  elastic  theory,  was  shown  to  be  capable  of  dup- 
licating the  action  of  such  systems  under  static  loads  with 
a reasonably  low  degree  of  error.  Numerical  values  for 
parameters  entering  the  model  were  obtained  from  the  simu- 
lation of  prototype  deflection  patterns.  Values  of  the  k 
parameter#  the  subgrade  modulus#  were  found  to  be  less  than 
51  pci.  at  the  start  of  the  prototype  tests.  The  magnitude 
of  this  parameter  was  then  found  to  decrease  as  trafficking 
progressed.  It  is  particularly  noteworthy  that  contrary  to 
prevailing  opinion#  model  behavior  is  extremely  sensitive 
to  the  magnitude  of  the  subgrade  modulus. 

Empirical  relationships  were  developed  which  related 
test  properties  to  all  model  parameters.  Use  of  these  re- 
lationships along  with  the  failure  criterion  established 
herein  indicated  that  the  performance  of  dual-wheel  tests 
can  be  predicted  with  a reasonable  degree  of  confidence. 
Apparently,  the  developed  procedure  is  less  satisfactory 
when  employed  for  the  prediction  of  single-wheel  tests. 

There  is  a strong  indication  that  the  variability  of  the 
number  of  passes  necessary  to  complete  one  coverage  in  the 
latter  case  may  be  a prime  factor  influencing  performance. 
This  is  particularly  true  as  input  energy  level  per  cover- 
age can  be  quite  different  for  seemingly  similar  situations. 

It  appears  that  the  effectiveness  of  the  end  joint 
connections  is  approximately  10  to  16  percent  less  than  that 
of  the  mat  elements.  Indications  suggest  that  some  increase 
in  the  serviceability  of  a landing  mat  system  can  be  achieved 
by  strengthening  the  end  joint  connectors. 

The  second  mechanistic  model  developed  is  capable  cf 
simulating  actual  loading  sequences  and  can  also  provide  a 
means  of  estimating  residual  deformations  of  the  mat  surface. 
The  potentialities  of  this  model  were  demonstrated  for  both 
single-wheel  and  dual-wheel  prototype  tests.  Lacking  ade- 
quate experimental  information  the  full  capabilities  of 
this  model  could  not  be  examined  in  this  study. 
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APPENDIX  A 


VARIATION  -.L  ITITHCD  OF  ANALYSIS  FOR  ELASTIC  FOUNDATIONS 


The  variational  method  of  analysis  for  elastic  founda- 
tions presented  here  was  developed  originally  by  Vlasov  and 
Leont'ev  [62].  Due  to  the  limited  circulation  of  this 
important  publication,  a derailed  formulation  is  given  for 
a single-layer  foundation  under  plane  strain  conditions. 

Consider  a homogeneous,  elastic  compressible  layer  of 
soil  of  thickness  H resting  on  a rigid  base  subjected  to 
imposed  loadings  which  produce  conditions  of  plane  strain, 
Figure  Al(a).  Under  such  conditions  if  the  horizontal  and 
vertical  displacements,  respectively,  u(x,y)  and  v(x,y),  at 
all  points  are  known,  the  stresses  and  strains  ir.  the  soil 
layer  may  be  obtained  from  established  stress-strain  rela- 
tionships. Following  Vlasov  and  Leont'ev,  the  unknown  dis- 
placements will  be  taken  as  the  finite  series 


n 


u(x,y)  » 

ii  i ui  <x>  yy> 

(i  = 1,  2,  3,...m) 

(Ala) 

v (x,y)  = 

* Vk(x)*  (y) 
k=  1 K K 

(k  » 1,  2,  3, . . ,n) 

(Alb) 

Here  it  is  assumed  that  the  functions  (y)  and  (v) , which 
represent,  respectively,  the  distributiin  of  the  displacements 
u (x,y)  and  v(x,y)  over  any  section  x = constant,  are  either 
known  or  may  be  assigned.  With  ^ (y)  and  * (y)  established 

as  dimensionless  functions,  the  functions  U.  (x)  and  VT  (x) 
will  have  dimensions  of  length.  The  lattez^are  seen  to 
represent  the  magnitude  of  the  displacements  u(x,y)  and 
v(x,y)  at  the  section  x = constant. 

Once  the  form  of  the  functions  (y)  and  (y)  are 
established  for  a given  physical  problem  the  functions 
U. (x)  and  V^ (x)  can  be  determined  from  equilibrium  condi- 
tions relative  to  a differential  strip  of  length  dx,  Figure 
Al(b),  and  the  stress-strain  relationships  developed  from 
the  conventional  theory  of  elasticity. 
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From  the  Hookean  stress-strain  relationship,  the 

stresses  may  be  expressed  as 


O =s 
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(A2a) 


(A2b) 


(A  3b) 


The  equilibrium  conditions  are  obtained  by  ecruating  the  work 
of  all  internal  forces  to  the  work  of  all  external  forces 
acting  on  the  strip  under  any  virtual  displacement. 

Consider  the  differential  strip  when  given  a virtual 
displacement  in  the  x direction,  u.  =♦  (y)  for  U.  = 1 where 
i can  have  rn  different  values,  andxa  virtual  displacement 
in  the  y direction,  v,  = (y)  for  V.  = 1 where  k can  be 

any  of  n virtual  displacements.  TheKdif ferential  strip 
thus  possesses  (m  + n)  degrees  of  freedom  in  the  xy  plane. 
The  external  forces  acting  on  this  strip  result  from  the 

do 


normal  stresses, 

3t 


I 0 + 

r x 


x 


3* 


dx,  the  shearing  stresses, 


Tyx' 


T + T— -1 

yx  3x 


dx,  and  from  the  distributed  applied  load 


whose  x and  y components  (per  unit  of  length)  are  p(x,y) 
and  n(x,y),  respectively.  The  internal  forces  are  caused 
by  the  normal  stresses  o and  the  shearing  stresses  t , . 

An  examination  of  Figured  A2(a)  and  . J (b)  provides  th!yfol- 
lowing  relationships: 


in  the  x direction,  Figure  A2(a), 

y y «o 
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no 


in  the  y direction,  Figure  A2(b), 
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Substitution  of  Eqns.  (A2)  into  Eqns.  (A3)  and  (A4) 
leads  to  a system  of  ordinary  differential  equations  in 
U.  (x)  and  V-,  (x) . This  system  consists  of  (m  + n)  equa- 

tions and  can  be  expressed  as 
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(h  * 1*  2,  3,...n) 

where  the  coefficients  are  obtained  from  the  following 
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acting'in  lit P Jitivl  “regain’  of  Ihof  oSrdLSa  axis 
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FIGURE  A2.  DIFFERENTIAL  ELEMENTS 
OF  SOIL  MASS. 


FIGURE  A3.  DISPLACEMENT  OF  SOIL 
SURFACE. 


The  applied  loads  p(x,y)  and  q(x,y)  are  eissumed  to  be 
distributed  over  the  plate  height  as  arbitrary  functions  of 
v.  However,  if  only  external  normal  ana  shearing  forces, 
q(x,0)  a.x  p(x,0) , act  on  the  upper  surface  of  the  layer, 


the  coefficients  p.  and 

3 


become 


Pj  = P(^)^j(0)  ; qh  = q(x)*h(0)  (A7) 


The  system  of  ordinary  differential  equations  in  U.  (x)  and 
V,r  (x)  obtained  from  Eqns.  (A5)  can  be  solved  to  obtain  the 
unknown  displacements  for  the  plane  strain  problem. 

If  it  is  assumed  that  the  horizontal  displacements  in 
the  compressible  layer  are  either  zero  or  are  negligible 
(y)  = 0#  the  displacements  u(x,y)  and  v(v,y)  can  be 
expressed  as 


and  the  function  VI (x)  represents  the  deflection  of  the 
foundation  layer  surface  as  illustrated  in  Figure  A3.  This 
last  assumption  docs  not  produce  a model  with  a unique  solu- 
tion as  the  c i iuc.'.sio. ileus  function  ♦ (y)  may  be  selected  in 

any  arbitrary  manner  provided  that  <|r  (H)  - 0. 

Applying  the  precious  assumption  of  a fixed  base,  Eqn. 
(A9)  reduces  to  th single  differential  ecruation 


1-W, 


■a — r V"(x)  - s V (x)  + 
2 111  li  l 


1-  ••  2 


“0 


E 


g - o 


(A12) 


where 


11 


*i6dy  * Bu“  |0  ♦! 


(ai:i 


Eqn.  (Al2)Emay  be  further  simplified  by  multiplying  each 


term  by 


(lmo  ) 


to  yield, 


where 


2tVj"  (x)  - kV  (x)  + q(x) 


(A14) 


Eor, 


li 


‘ " ?nv; ) ’ k 


Vn 
i-  *i 


(f.15) 


and  r^  and  s ^ are  defined  by  Egns.  (A13) . 

Differential  equation  (A14)  defines  a model  whereby  the 
vertical  displacements  of  the  surface  of  the  soil  are  related 
to  the  load  applied  on  the  surface.  This  model  not  only  con- 
siders the  compressive  strains  in  the  elastic  foundation  but 
also  the  shearing  strains  within  the  foundations. 

With  the  displacements  of  all  points  expressed  by  Ecn. 
(All)  and  employing  the  stress-strain  relationships  of  the 
classical  theory  of  elasticity  the  normal  and  shearing 
stresses  may  be  obtained  from  Eons.  (A5)  as 


0 a 

X 


EoyO 
1- WO 


v,  (x)*!(y) 


(Al6a) 
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0 = 


E0 

— T V (x)  (y) 

1-M0  1 1 


(Al6b) 


<*> 


(Al6c) 


Selecting  any  section  x = constant,  the  total  force  T in  the 
x direction  is  given  by  T.  “ o <dy  (i  « 1,  2,  3,...m) 

1 | XI 


and  the  total  force  in  the  y direction  is  given  by 


1 f 

e y dil: 


L 

a JoV*h6dy 


(h  * 1,  2,  3,...n).  Pleading  the  assump- 


tion of  negligible  horizontal  displacements  and  applying 
Eqn.  (A16c)  there  results 


n 

f c ♦ . 

JO  x l 


idy  «*  0 


(Al7a) 


8 j B , Tyx*j8dy  = vi(x)  f (y)dy  " 2tVi  (x>  (Al7b) 

J Jo 

The  foundation  model  is  now  completely  defined  once  the 
dimensionless  function  y.  (y)  has  been  established. 

Tl 
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APPENDIX  B 


i::n:;n::  beam  ok  a:;  elastic  foundation 


Superimpose  on  the  surface  of  the  soil  model,  as 
defined  in  Aopcnuix  A,  an  infinitely  long  beam  which  obeys 
Eulerian  conditions.  Let  p(x)  represent  the  applieo  normal 
load  on  the  beam,  q(x)  represent  the  reaction  due  to  the 
elastic  foundation  and  V (x)  represent  the  beam  deflection. 
The  bending  of  the  beam  will  be  described  by  C 63  i 

4 

El  A-V-tgl  » p (x)  - q(x)  (Bl) 

dx 

In  Appendix  A the  relationship 

4 

- 2t  2-2-lp-  + kV (x)  = q(x)  (B2) 

dx 

was  developed  to  relate  surface  deflections  of  the  soil  to 
superficial  loads.  Assuming  that  the  beam  and  the  soil  are 
always  in  contact,  Eqns.  (Bl)  and  (B2)  can  be  combined  by 
elimination  of  the  interfocial  normal  force  cr(x)  , to  yield 

El  -■  - 2t  ^ + kV (x)  * p(x)  (B3) 

dx  dx 


Implicitly,  this  last  assumption  requires  the  development 
of  normal  tensile  stresses  across  the  beam-soil  interface. 
Fortunately,  any  development  of  these  stresses  generally  will 
occur  in  a region  somewhat  removed  from  the  load  and  will 
be  of  small  magnituce.  ‘ 


For  computational  expediency,  Vlasov  §nd  Leont'ev 
defined  a non-dimensional  coordinate,  n = 7-  , as 

li 


(B4) 
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APPENDIX  3 


7”  *>T  “•  «**fT  *■' 

Xi-A.  J.J.*  J-  J.  — 


beam  on  an  elastic 


FOUNDATION 


Superimpose  on  the  surface  of  the  soil  model,  as 
defined  in  Appendix  A,  an  infinitely  long  beam  which  obeys 
Eulerian  conditions.  Let  p(x)  represent  the  applied  normal 
load  on  the  beam,  q(x)  represent  the  reaction  due  to  the 
elastic  foundation  and  v(x)  represent  the  beam  deflection. 
The  bending  of  the  beam  will  be  described  by  C 63  J 

4 

El  a p(x)  - q(x)  (Bl) 

dx 

In  Appendix  A the  relationship 

4 

- 2t  + kV (x)  ■ q (x)  (B2) 

dx 

was  developed  to  relate  surface  deflections  of  the  soil  to 
superficial  loads.  Assuming  that  the  beam  and  the  soil  are 
always  in  contact,  Eqns.  (Bl)  and  (B2)  can  be  combined  by 
elimination  of  the  interfocial  normal  force  cr(x),  to  yield 

El  - 2t  - + kV  (x)  = p (x)  (B3) 

dx  dx 


Implicitly,  this  last  assumption  requires  the  development 
of  normal  tensile  stresses  across  the  beam-soil  interface. 
Fortunately,  any  development  of  those  stresses  generally  'rill 
occur  in  a region  somewhat  removed  from  the  load  and  will 
be  of  small  magnitude.  ‘ 


For  computational  expediency,  Vlasov  $nd 

defined  a non-dimensional  coordinate,  n = y , 

L 


Leont'ev 

as 


(B4) 
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where  E is  Young's  modulus  of  the  beam 

I is  the  moment  of  inertia  of  the  beam 
£0  is  Young's  modulus  of  the  soil 
Hg  is  Poisson's  ratio  of  the  soil 
4 is  the  width  of  the  beam 

Employing  the  r.ono imensional  coordinate  and  dividing  Eon. 
(D3)  by  the  beam  rigidity,  the  governing  differential  ecua- 
tion  is  established  as 

Oijl  . 2r2  — V-l^  + s"v(  J (B5) 


where 


2 tL 

r air 


(B6a) 


4 pj 

S2  = ||-  = 2L  r dy  (B6b) 

If  a concentrated  forced,  Figure  Bl,  acts  at  the  origin 
of  coordinates,  the  following  homogeneous  differential  equa- 
tion is  applicable  for  all  points  except  tie  origin 


2 * 


.fLvj.iL  _ 2r  + s4v{n)  *=  0 


The  corresponding  characteristic  equation  is 


m4  - 2r2m2  + s4  = 0 (E8) 

As  r and  s are  positive,  the  ratio  — must  be  positive  and 
the  solution  of  the  characteristic  equation  may  take  one  of 
three  forms  depending  on  the  relative  magnitudes  of  r and  s. 
The  sc'  hion  of  prime  importance  for  the  present  study  is 
obtains,  for  the  condition  when  s > r.  Physically,  this 
condition  implies  that  the  majority  of  the  supporting  capac- 
ity of  the  soil  is  developed  by  normal  stresses  rather  than 
by  the  shearing  stresses  within  the  foundation.  Close  inves- 
tigation of  Eqn.  (37)  reveals  that  for  the  condition  r = 0, 
the  model  becomes  identical  to  the  one  which  employs  the 
Winkler  hypothesis.  At  the  other  extreme  as  s approaches 
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FIGURE  Bl.  CONCENTRATED  LOAD  ON  MODEL. 


FIGURE  B2.  UNIFORM  LOAD  ON  MODEL. 


zero,  the  model  becomes  one  for  a beam  supported  by  a mem- 
brane. ;J.l  further  formulation  is  restricted  to  the  condi- 
tion • hen  s > r.  Tor  this  condition  the  roots  of  the  charac- 
teristic equation  are 


m = + o + i?  (BS) 

where  S and  I are  real  and  positive 

2 2 2 2 

; . M-tL.  , | . S-=£-  (bio) 

The  solution  of  Eqn.  (B7)  then  has  the  general  form 

«■  m m 

VtiJ-Cje’^sin^+Cje-  cos  ^’VCje  8nsii5n +C^e®T,cos?r>  (Bll) 

Due  to  symmetry  only  that  portion  of  the  beam  for  n > 0 
need  be  considered.  Since  the  deflection  of  the  beam  must 
approach  zero  as  the  nondimens ional  coordinate  becomes 
large  ( rr—  ) , the  following  must  occur 


C3  » C4  ' 0 (B12) 

The  deflection  equation,  Eqn.  (Bll) , then  becomes 


v<n> 

where 

* C F + CP 

(B13) 

11  2 2 

F. 

F 

o 

= e~8nsinBn 

(Bl4a) 

« e”8ncos5n 

(Bl4b) 

C 

The  remaining  constants  C 

and  C can  be 

determined  from  the 

conditions  at  the  origin 
Thus  for  n - 0 

ielatec?  to  the 

slope  and  the  shear. 

e(o)  = 

. 1 = o 

L ^ 

(B15) 
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(BIG) 


K(0)  = - 


FI 

T 

L 


I Id  y 

7L~ 


2r 


M±n i]=  _ 

dti  J 


£ 

2 


It  should  bo  noted  that  the  shear  at  any  section  is 
cor posed  of  two  distinct  components.  The  first  tern  Eqn. 
(BIG) , results  fro-1  the  shear  force  developed  in  the  bear 
whereas  the  fr'or.r  tern  arises  from  the  shear  force  developed 
in  the  soil.  The  latter  term  is  expressed  in  integral  form 
in  Eqn.  (/a 7b) . 

Substituting  the  appropriate  derivatives  of  Ecn.  (B13) 
into  Eons.  (Bl5)  and  (BIG)  and  solving  for  the  constants 
produces 


C1  = 


PL' 


4||s2EI 

PL3 
-«2. 


(B17a) 


2 43s2  El 


(Bl7b) 


Thus  the  general  equation  for  the  deflection  of  the  beam 
subject  to  a concentrated  load  becomes 


v<n>  - -—i — t!  Fi(”>  + 5 F,  <i>) 

%0  s El  1 2 


(B18) 


Modification  of  Eqn.  (B18)  is  necessary  when  uniformly 
distributed  loads  are  encountered.  Referring  to  Figure  B2, 
three  expressions  are  developed  for  the  beam  deflection 
which  depend  upon  the  location  of  the  uniform  load  relative 
to  the  point  at  which  the  deflection  is  sought.  In  the 
development  of  each  expression  a differential  element  is 
used  in  conjunction  with  Eqn,  (BIS)  which  is  then  integrated 
over  the  length  of  the  load.  Thus  for  0 < n < - 

- L 

V<")  >■',(!  -")]+  25f/f  -n,W! -4 


(Bl9a) 


f°r  I < V L 
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v ( n)  = 


L3n{x) 

4 S]f4LI 

-2S[F2("rf)  + ps(t-0]}  (B19b> 

b 

an..  -O v 7 < n, 

~ 3 


v(n) 


= L3c(x) 

455s 


(Bl9c) 

With  the  use  of  Eqns.  (B19) , the  deflection  of  the  beam  at 
any  point  can  be  determined  for  a given  uniform  load.  After 
employing  superposition,  any  number  of  uniform  loads  may  be 
handled.  The  parameters  5,  g,  s 4,  and  r 2 are  all  functions 
of  the  displacement  distribution  function  ^ (y) . 


APPENDIX  C 


STATE  VARIABLE  FILTER  FORMULATION 


The  procedure  for  identification  of  the  model  parameters 
presented  below  follows  Kohr  [30] . 

From  the  data  available  from  prototype  tests  [ 12  1 a 
system  (relationship)  can  be  established  between  the  known 
input  and  the  resulting  output.  The  form  of  this  system  has 
been  postulated  as  Eqn.  (B3) . For  this  system  to  reflect 
adequately  the  behavior  of  landing  mats,  the  model  parameters 
must  be  representative  of  the  test  conditions. 

By  writing  Eqn.  (B3)  as 


A . a,  div  + btv 


El 


where 


a*  - § ••  b*  ■ fee 


separation  of  the  unknown  quantities  can  be  achieved . The 
term  on  the  right  is  exclusively  input  data;  that  is,  load 
configuration,  expressed  by  step  functions,  and  mat  rigidity. 

Let  the  "model  residue",  c,  at  a point  be  defined  as 

c = £-?  - a*  ^-4p  + b*V  - (C 2 ) 

dx  dx  *x 

and  the  error  functional,  f ( c) , for  all  points  be  defined 
as 
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If  V, 


can  he  obtained,  then  the  "model  residue 


d 

and 

dx 

and  the  error 


proper  sclocti 
are  called  rh.e 


functional  can  be  made  a minimum  through  the 
on  of  o*  and  b*.  T:hen  this  occurs,  a*  and  b* 
"best''  cctirr.atcr  of  the  model  parameters. 


The  derivatives, 


d2y 

2 


d 

•nd  — r required  in  Eqn.  (C2)  may 

dx* 


be  chtnfricJ  from  the  actual  deflections  by  use  of  a state 
variable  filter  [20].  For  this  operation,  it  was  convenient 
to  represent  the  "average  deflection"  patterns  by  a continu- 
ous function  rather  than  by  a series  of  discrete  points.  To 
obtain  deflection  as  a function  of  position,  the  selected 
points  from  the  "average  deflection"  patterns  were  approxi- 
mated by  a polynomial  using  a least  squares  technique.  The 
criterion  established  for  the  polynomial  was  that  it  approxi- 
mate the  discrete  points  witnin  a selected  standard  devia- 
tion of  0.05.  For  the  three  deflection  patterns  investigated 
by  this  procedure,  a sixth  degree  polynomial  was  found  to 
satisfy  the  standard  deviation  criterion.  The  polynomial 
thus  obtained,  denoted  by  V,  was  established  as  the  expres- 
sion for  the  test  output  and  served  as  the  forcing  function 
for  the  state  variable  filter.  Following  Kohr,  the  state 
variable  filter  employed  was  a differentiel  eouation  of  one 
order  higher  than  the  highest  derivative  of  the  model.  Thus 
the  filter  selected  was  taken  as 


d 5V 
dx5 


dV 


+ q. 


dx 


T + <J 


+ q 


^5 

dx 


+ Vc  = V 


(C4) 


where  the  coefficients,  q,  , q # q , snd  q , for  this  filter 

, 12  3 4 

were 


q = 3.29710 

1 

q = 4.89532 

2 

q = 4.53948 

3 


(C5) 


4 = 2.15678 

These  coefficients,  called  tfte  Paynter  filter  coefficients, 
were  also  given  by  Kohr  £ 30  ) . Solution  of  Eqn.  (C4)  yielded 
d^c  d*Vc 

the  values  of  — , • 

dx*»  dx* 
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Once  these  values  were  known,  they  were  substituted  into 
Eon.  (C2)  and  the  error  functional,  Eqn.  (C 3 ) , was  minimized 
by  using  the  "descent  equations"  [30,  66].  These  equations 
for  a two  parameter  system  are 


= « K JlIUI 

dx  a 3 a* 


2K  r 

ae 


JSu. 

ga* 


db* 

dx 


iliil 

3 b* 


3c 

= - 2Kb£  3^7 


(C6) 


where  K and  K.  are  arbitrary  positive  constants.  These 
constants,  called  the  "identification  gains",  stipulate  the 
speed  at  which  the  minimum  is  approached. 


Noting  from  Eqn.  (C2)  that 


3c 
3 a* 


3c 
3 b* 


(C7) 


and  making  the  appropriate  substitutions,  the  values  of  a* 
and  b*  are  obtained  by 


a* 


dx 


b*  - - 2KbJLeVcdx 


(C8) 


Solution  of  Eqns.  (C8)  was  obtained  by  evaluating 
gral  between  the  points  of  zero  deflection  and  by 
iently  setting  K and  eaual  to  0.5. 


the  inte- 
conven- 
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APPENDIX  D 


STEEP  DESCENT  METHOD 


As  demonstrated  in  Appendix  B/  the  behavior  of  the  mat- 
soil  model  is  governed  by  an  expression  of  the  form 


vtheoretical  = G(VT> 


(Dl) 


where  Eo  and  y are  unknown  parameters  and  G is  the  differen- 
tial equation,  Eqn.  (B5) . The  behavior  of  the  prototype 
tests  can  be  simulated  by  the  model  only  if  the  parameters 
E o and  y are  representative  of  the  test  subgrade  conditions. 
To  determine  E and  y , a criterion  is  established  as  to 
minimize  the  s&m  of  the  squares  of  the  deviations  of  the 
theoretical  deflections  from  the  experimental  deflections. 

As  the  "average  deflection"  patterns  [12]  have  been  defined 
by  discrete  points,  let  the  error  functional  to  be  minimized 
be  defined  as 


n , , 

f(c)  =*  Z Cj  (l  * 1#  2,  3,...n)  (D2) 

i= ) 

where  t represents  the  deflection  deviation  at  a discrete 
point  and  n is  the  number  of  discrete  points.  Equation  (D2) 
can  be  minimized  by  a "trial  and  error"  procedure;  however, 
a more  direct  procedure,  "steep  descent  methoo  l 40,  66  J, 
which  is  based  upon  the  behavior  of  the  error  functional  can 

be  employed. 


Five  combinations  of  E and  y are  selected  in  a manner 
as  to  form  a five  point  grid,  Figure ^1 ,_ and  the  error 
functional  is  computed  for  each  combination,  from  this 
grid  a spatial  surface  is  developed  which  reflects  the 
local  behavior  of  the  error  functional.  Assuming  that  the 
error  functional  possesses  a bowl-like  surface  in  the 
vicinity  of  the  minimum  a second  order  surface  of  the  form 
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(D3) 


Z = ay  + 2hyE0  + + 2gy  + 2fEQ  + c 

is  fitted  to  the  grid.  The  influence  of  the  quadratic  terns, 
in  iqn.  (L'3)  , '/ill  become  negligible  in  a region  removed 
fro.r.  the  minimum  and  a plane  containing  the  central  grid  can 
adequately  describe  the  local  behavior  of  the  error  functional. 
The  equation  of  this  plane  is  taken  as 


Z » 2gy  + 2fEQ  + c (D4) 

where  g and  f are  obtained  from  Eqn.  (D3)  and  c is  the  value 
of  the  error  functional  at  the  central  grid  point.  The  coef- 
ficients, g and  f are  determined  in  this  manner  so  as  to 
obtain  the  most  representative  plane.  The  direction  of  the 
line  of  maximum  slope,  steepest  descent,  in  this  plane  is 
given  by 


tan8«l£  (D5) 

2g 

A minimum  of  the  error  functional  will  be  rapidly  attained 
by  proceeding  in  the  8 direction  from  the  central  grid 
point. 


Reasonably  small  increments  should  be  used  when  pro- 
ceeding along  the  steep  descent  line  to  avoid  bypassing  the 
minimum.  Also,  it  is  advisable  to  initiate  the  procedure 
from  various  locations  on  the  error  functional  surface  to 
determine  whether  the  minimum  obtained  is  global  rather 
than  local. 


APPENDIX  E 


’O^ilULATIO”  FOR  KCllEKT  TRANSFER  INVESTIGATION 


To  facilitate  the  computation  of  the  deflections  of  a 
hinged  beam,  two  semi-infinite  beams,  one  on  each  sice  of 
the  hinge,  are  loaded  with  a concentrated  moment  as  shown 
in  Figure  El.  Two  distinct  expressions  are  developed,  one 
for  the  deflection  of  the  semi-infinite  beam  and  the  other 
for  the  exposed  ground  surface.  The  latter  is  necessitated 
to  insure  the  continuity  of  the  soil  media. 

With  reference  to  Figure  El,  it  is  apparent  that  Egn. 
(B7)  is  applicable  for  r\  > 0 and  that  the  deflection  for 
this  region  can  be  expressed  by  Eqn.  (Bl3) . The  governing 
differential  equation  for  the  exposed  foundation  surface 
from  Eqn.  (A14)  for  x < 0 is 

2 

2t  - kV  (x)  - 0 (El) 

dx  2 

of  which  the  general  integral  is 


V (x)  » D e"8*  + D e0X  (E2) 

l 2 

where 


Applying  the  boundary  condition  that  the  deflection  V(x) 
approaches  cero  as  x approaches  minus  infinity  yields 


= 0 


120 


FIGURE  El.  MOMENT  LOADING  ON  END  OF 
SEMI  - INFINITE  BEAM. 


The  deflection  of  the  foundation  surface  for  x < 0 is 


V(x)  = D e (E4) 

2 

To  determine  the  constants  Cj  , Cj  , of  Eqns.  (B13) , and  D2  , 
the  following  boundary  conditions  are  used  at  the  origin 


V(-0)  = V (0) 


M = - — V (0)  (E5) 

L* 

N(-0)  " N (0) 

where  N(-0)  = 2tV’  = ?a  tD  is  the  generalized  shear  at  minus 
4ero.  Applying  the  first2of  Eqns.  (E5)  produces 


D2  - C2  (E6) 

From  the  second  of  Eqns.  (E5)  after  substituting  the  second 
deviatives  of  the  F functions,  Eqns.  (Bl4) , at  the  origin 
there  results 


ML—  — r 

El + r C; 


where  S«  S»  and  r2  are  as  defined  in  Eqns. 
The  third  boundary  condition  giyes 


(E7) 

(B6a)  and  (BlO) . 


2o  tD  2 
3 


d3y(o) 
d "3 


2r 


dy  (0) 


Multiplying  by  -L  /El  and  making  the  appropriate  substitu- 
tions for  the  derivatives  of  the  F functions,  Eqns.  (B14) , 
one  obtains 


-2  cfcL  D2 
— 


,-2  *2 


= [5  (3S‘  -r)-2r"eJ  C +[5  (3?z  -J2)  +2r2,S]  C2  (E8) 


I 
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From  the  simultaneous  solution  of  Eqns.  (E6)  , (E7),  and  (E6)  , 
the  intv.  jre  _i  e.i  constants  are  found  to  be 


D 


2 


[B  s2j 


where 


D = - 


(£9) 


2 

and  s is  as  defined  by  Eqns.  (B6b) . Thus  the  deflection 
for  the  semi-infinite  beam  is  given  by 


v(n)  = C F + C F 
11  2 2 

and  the  deflection  of  the  ground  surface  by 


V(x) 


2 

ML 

EID 


?s2e 


ax 


(E10) 


(Ell) 


If  the  semi-infinite  beam  extends  to  the  left  of  the 
origin,  the  coordinate  system  is  reversed  and  Eqns.  (E10) 
and  (Ell)  are  then  applicable. 
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APPENDIX  F 


COMPUTER  PROGRAM  FDR  PARAMETER  IDENTIFICATION 


******************************* ***• 

PARAMETER  IDENTIFICATION  PROGRAM 


INFINITE  3EAM  ON  A VLASOV  AND  LEONTE V MODEL 
DOWNWARD  DEFLECTION  IS  TAKEN  AS  POSITIVE 


DEFINITION  OF  VARIABLES  USED 

KNK  * NUMBER  OF  TEST  SECTIONS  TO  HE  INVESTIGATED 
KNN  ■ INTEGER  COUNTER 

E ■ YCUNCS  MODULUS  OF  THE  MAT  MATERIAL  I Kl PS/SO. I N. I 
I » MOMENT  OF  INERTIA  OF  THE  MAT  PER  FOOT  OF  WIDTH  UN.  **A) 
CELTA  * WIDTH  OF  THE  MAT  - TAKEN  AS  THE  LENGTH  OF  A «ECTA\3LC 
WHOSE  AREA  IS  EQU4L  TO  THAT  OF  THE  TIRE  PRINT  AND  WHOSE  wIDTH 
lb  THE  MAXIMUM  WIDTH  OF  THE  TIRE  PRINT  UN.) 

KKK  * NUMBER  OF  SETS  OF  DATA  PER  SECTION  - A SET  IS  A 
CEFLECTION  PATTERN  FOR  A PARTICULAR  COVERAGE  WITHIN  THE  TEST 
NN5  • NUMBER  OF  STATIC  UNIFORM  LOADS  - WHEEL  LOADS 
POIS  * POISSONS  RATIO  OF  THE  SOIL 

H * THICKNESS  OF  SOIL  LAYER  - 1CCO  USE3  FOR  H«  INFINITY 
SAI  « DISTANCE  TO  THE  BEGINNING  OF  STATIC  UNIFORM  LOAD  FROM 
TIE  ORIGIN  (IN.) 

SQG  « MAGNITUDE  OF  STATIC  UNIFORM  LOAD  (KIPS/INCH) 

sei  * cistance  to  the  end  OF  STATIC  uniform  loao  from  the 

ORIGIN  (IN.) 

NN  * INTEGER  COUNTER 

IH  « INTEGER  VARIABLE  ASSOCIATED  WITH  THE  TYPE  OF  STRESS 
DISTRIBUTION 

1 » L INEAR 

? « HYPERBOLIC  - FINITE  DEPTH 
3 « HYPERBOLIC  - INFINITE  DEPTH 
NOP  « NUMBER  OF  DEFLECTION  POINTS  PER  SET 
Il_*  SECTION  NUMBER 

12  * LANE  NUMBER  " ' ' 

13  * ITEM  NUMBER 

14  « IOC AT  ION  NUMBER 

1 * CFNTER  OF  PANEL'  ' ~ 

2 * QUARTER  POINT  OF  PANEL 

_ 3_^_  JOINT 

ICOVER  * NUMBER  OF  COVERAGES  FOR  THE  SET  “ 

XS  « DISTANCE  FROM  ORIGIN  TO  POINT  AT  WHICH  DEFLECTION  IS  GIVEN 

I In.  ) 


A 1 
A 2 
A 3 
A A 
A 5 
A 6 
A 7 
A 8 
A 9 
A 13 
A 11 
A 1? 
A 13 
A 1 A 
A 15 
A 15 
A 17 
A 18 
A 19 
A 20 
A 21 
A 22 
A 23 
A 2A 
A 25 
A 25 
A 27 
A 28 
A 29 
A 30 
A 31 
A 32 
A 33 
A 3A 
A 35 
A 35 
A 37 
A 38 
A 39 
A AO 
A A 1 


w‘Onr»t  j o or>  ono  o o o w o o o noo  non  or»o 


c 

c 


4 


OS  * Mt*SJRFD  OEFLFCTION  AT  tHE  POINT  XS 

CAMMA  * PATE  OF  DECREASE  OF  STRESS  WITH  DEPTH 

EC  * YOUNGS  MODULUS  OF  THE  SOIL 

********* ********* ****************** *********••**»•••*# ****** 

PEAL  I,L,kS,mOME,MSMAI,MOmES 

riMFNSH'N  SAKS),  SQO(S),  SB  1(5}*  XS  f 20)  * 05(20),  X0EFK20) 

COMMON  L, A, e,A„),S2,S6,R2,E, I *N5» XI 60) ,AII5) ,Q0I5) »BI (5» ,DEFLI60),$ 
IL0PEI60  ),K0MF(6C),SHEAR16C) 

RF  AT  IS, 1S»  KNR 
KNN-0 
KNN-KNN ♦ 1 

IF  (KNN.EC.IKNK*l))  STOP 

RE  AO  INFORMATION  ABOUT  BEAM  PROPERTIES 

READ  IS, 16)  E, I, DELTA 
1*  !*DELTA/12. 

READ  INFORMATION  ABOUT  FOUNDATION  PROPERTIES 

RF AC  (5,20)  KK<,NN5,P0IS,H 

READ  LOADING  INFORMATION 

RF  AD  IS, 16)  I S AI I < ),  SOOK  ) , S8I l K ) ,K* 1 ,NN5) 

NN*  0 
NN* KN+1 

IF  (NN.EC.UKK*!))  GO  TO  1 

READ  INFORMATION  ABOUT  DISTRIBUTION  OF  STRESS 


READ  (5,17)  I H, NOP 

REA  D TEST  S EC  T ION  IOENTIFIC  A T ION  ...  _ 

REAC  (5,18)  11,12, 13, I A, ICO VER 


REAITdAI  A POINTS 

REAC  (S,  19  ) ( XSU  ) , DSI  J ),J » 1,NDP  )_ 

WRITE  I 6, 21 ) 

WRITE  (6,22)  II,  12,  13, H, (COVER 

WR I TE  16,23) 

W^IT6(6,2A*) 

WRITE  (6,25) 

GAMMA- l, 


ERRS  ■ STORED  VALUE  OF  ERROR  - SET  INITIALLY  TO  10000. 
Errs*T66oo. 

L- ( 1 2 . * E ♦ I * I l.-P0IS»*2.))/I EO*OELTA) )»» 0.33333333 


A A 2 
A A3 
A AA 
A A5 
A A6 
A A7 
A AS 
A A9 
A SO 
A SI 
A 52 
A 53 
A 5A 
A 55 
A 56 
A 57 
A SB 
A 59 
A 60 
A 61 
A 62 
A 63 
A 6A 
A 65 
A 66 
A 67 
A 68 
A 69 
A 70 
A 71 

A n 
A 73 
A 76 
A 75 
A 76 
A 77 
A 78 
A 79 
A 80 
A 81 
A 82 
A 83 
A 86 
A 85 
A 86 
A 87 
A 88 
A 81 
A 9C 
A 91 
A 92 
A 93 
A 96 
‘95 


o 0!0  ® 


C SELECT  EQUATION  FOR  PROPER  STRESS  OISIRI3UTION 
C 

IF  ( IH.fC.l  ) 30  TO  5 
IF  ( 1H. EC.? ) CJ  TJ  6 
IF  ( IH.EQ.3)  GO  TO  7 
*>  31 1 * CEL  T A*H/  3 . 

Sll-CELTA/H 

T * ( £0*R 1 1 )/ ( 4 . ♦( l.+POISI) 

KS»IEO*SU)/(  1.-P01S**2.) 

R2* ( T *L  **2.  )/< E*l) 

S6* <KS*L*»6. )/( EM ) 

GO  TO  8 

c 

6 GHL*(GAHMA*H)/L 

PSIK«0.5*GHL*I (SINHIGHL >*C0 SHIGHL J+GHL I /I  Si NH (GHt ) **2. )) 
PSIT-1.5/GHLM  (SlNHl  GHL  )*C0  SHIGHL  )-GHU  /(  SI  NH  (GHU  »*2.  )I 
T* ( EO*OELT A*H*PS IT  I/I 12. M l.+POI SJ ) 

KS«(E0*CELTA*PS!<  I/O-*  I l.-P0IS«*2.)) 

R2* I H*I l.-P0IS)*PSITI/I6.*L) 

S6* (2,»L*PS  IK l/H  ' - 

GO  TO  8 
C 

I ' KS*  ( E0*‘0ELT  A*GAMHA  |/(  2.*LM  1 ,-P0Ts"**2.  ) ) 

T»( E0*0ELTA#L )/< 8.*GAMMA*t l.4P0l$|) 

R2« I l.-POIS  )/ 1 6.*GAHMA I 
$6*gahma  '■ 

S2-S6M0.5 

t€S T T0~6 EfTfRM  lNE*!?TOL\OWlN~G“tQuSTlWS  A rF"  A PP U C A BLE 
IF  IR2.GT.S2j  GO  TO  1* 

*■ I IS2*R2  )7 2. )**6.5  — 

6«  HS2-R2 )/  2.  ,5 

AB»2.*A*B 

N5*  NN5  — 

00  9 K«  I, NS 

AKK  )»SA  I(  K ) 

QOIK  )*SQO( K j * " 

9 Bl  1 K )*S  B 1 1 K ) 

.00  10  J*  1 * NDP 

x(j)»xsTJ)  ~~ — — 

CALL  ULOAD  IJ) 

XOEFLIJ )»OEFLI J ) 

10  CGNT  INU  E — — 

ERR0R*0  .0 

00  11  J»1 ,NOP 

II  ERROR*  ERR0R4A8SI  XOEFLIJ  l-DSIJ)  )**2. “ 

IF  IERROR.GT.ERRSI  GO  TO  13 

ERRS«ERROR 

IF  I ERRS  .L T . 1 . 1 1 GO  TO  12  — 

ECS*  EO 
SKS«KS 

TS*T  

EO*  E040 .05 


A 96 
A 97 
A 98 
A 99 
A 100 
A 101 
A 102 
A 103 
A 106 
A 105 
A 106 
A 107 
A 108 
A 109 
A 110 
A 111 
A 112 
A 113 
A 116 
A US 
A 116 
A 117 
A 118 
A 119 
A PO 
A 121 
A 122 
A 123 
A 126 
A 12S 
A 126 
A 127 
A 128 
A 129 
A 130 
A 131 
A 132 
A 133 
A 136 
A 1 35 
A 136 
A 137 
A 138 
A 139 
A 160 
A 1 M 
A 162 
A 163 
A U6 
A 165 
A 166 
A 167 
A 168 
A 169 


GO  TO  b 

eos-eo 

SKS 

TS*T 

E0*£0»0 .01 
GC  TO  b 

S'J8fl»SKS/CELTA 
Tf  * TS/ CELT  A 

WtUTE  t 6,  ?6  ) EOS«  GAMMA,  SUBR  , TT,u ,ERR  S 
IF  1GAMHA.GE.8.) L GO_T0_2 

gamha*gamma«i. 

GO  TO  3 

wafrcTfcWTl 


_ FORMAT  U10) 

FORMAT  < 3F20.51  * 

FORMA!  <21151 
FORMAT  <51101 

FORMAT  < 2F20.5 1 " “ " 

FORMAT  <2115, 2F15.51 

_ FORMAT  J//10X,  INSECT  ION,  10X,  AHLANE , 1 OX,  AHl  TE  M,10X,8HL0CATl  0N,nx,9 
lHCOWERAGEST  ' ' 

FORMAT  <115,116,114,116,1191 

FORMAT  </5X,9HE  OF  SQ1L.9X, 5HG4PMA, 5 X , RHSUB3RA0E ,1 2 X ,IHT ,1 3 X ,IML, 1 
12X , SHERROR  1 ~ 

FORMAT  < 37X , 7HMOOULUS 1 

FORMAT ( 3X,  1 1HK  IPS/SO.  IN.,  20X,  15MKIPS/S0.  I N. /l_N.  ,5X,8HK1 PS/lN. 21 

FORMAT  < 6 F 1 5 , 5 1 

FORMAT  I/10X.34HR2  IS  GREATER  THAN  S2  NO  SOLUTION! 


APPENDIX  6 


COMPUTE*  POGRAM  FOR  PERFORMANCE  PREDICTION 


PERFORMANCE  PREDICTION  PROGRAM 

C INFINITE  BEAM  OH  A VLASOV  AND  LfcONTEV  MODEL 

C HYPERBOLIC  STRESS  01STR l BunON 

C DOWNWARD  "DEFLECTION  IS  TAKEN  AS  POSITIVE 


C DEFINITION  OF  VARIABLES  USED 

C KNK  ”•  NUMBER  OF  TEST  SECTIONS  TO  BE  fNVE  ST^TTfC 

C KNN  - INTEGER  COUNTER 

C II  » SECTION  NUMBER 

C 12  » LANE  NUMBER ” 

C 13  ■ ITEM  NUMBER 

C. _ I NCA_f^  ACTUAL  NUMBER  OF  COVERAGES  AT  WHICH  THE  SECTION  FAILED. 

C E • YOUNGS  MODULUS  OF  THE  MAT  MATERIAL  ( KIPS/SO. I N. I 

C I * MOMENT  OF  INERTIA  OF  THE  MAT  PER  FOOT  OF  WIDTH  UN.***) 

C OELT A « WIOTH  OF  THE  MAT  - TAKEN  AS  THE  LENGTH  OF  A RECTANGLE 

C WHOSE"  AREA  IS  E3UAI  TO  THAT  OF"  THE  TIRE  PRINT  AND  WHOSE"  WIDTH 

C IS  THE  MAXIMUM  WIOTH  OF  THE  TIRE  PRINT  (IN.) 

C NN5  » NUMBER  OF  STATIC  UNIFORM  LOADS  - WHEEL  LOADS 

C WGAMMA  • THE  PRODUCT  OF  THE  AVERAGE  WATER  CONTENT  (PERCENT) 

C AND  THE  AVERAGE  ORY  UNIT  WEIGHT  OF  THE  TOP  18  INCHES  CF  SClL. 

C C.B.R.  VALUE  FOR  THE  TOP  IB  INCHES  CF  SCI  L 

C WL  • MAGNITUDE  OF  THE  LOAD  - FOR  DUAL* WHEELS”  USE  TOTAL  LCAD  ON 

C THE  ASSEMBLY  - FOR  SINGLE  WHEEL  USE  WHEEL  LOAD.  (IN  KIPS) 

_C CA  • CONTACT  AREA  (IN  S3. IN.)  - FOR  BOTH  SINGLE  AND  DUAL 

C ASSEMBLY  USE  CONTACT  AREA  OF  ONE  WHrEL.  

C TP  « TIRE  INFLATION  PRESSURE.  (IN  LBS./SO.IN.) 

C SAI  » DISTANCE  TO  THE  BEGINNING  OF  STATIC  UNIFORM  LOAD  FRCM 

C THE  ORIGIN  ( IN.)  

C SCO  » MAGNITUDE  OF  STATIC  UNIFORM  LOAD  (KIPS/INCH) 

c SB]  _■  DISTANCE  TO  THE_ENO  OF  STATIC  UNIFORM  LOAD  FROM  THE 

C ORIGIN  (IN.)  — 

C POIS  • POISSONS  RATIO  OF  THE  SOIL 

_C G/!MM_A_«  RATEOF  DECREASE  OF  STRESS  WITH  DEPTH 

C EO  « YOUNGS  MODULUS"  OF”  THE  SOIL  * 

C **********«*»*********•••«***••••**»•********«»•****„«•«„«, 

PCAL  I.KS.Lt  INTKS*  II.MOME 

DIMENSION  SAU  5),S00(  5),SBI< S)  

COMMON  X(60  )• All  5), QO(  5),BI ( 5),L iA,B ,AB ,S2, S*,R2 ,E ,1 *N5  »0EFL(6D) , 

1 SLOPE]  60),  HOME]  60 ), SHEAR ( 6C) 


■ 4 


REAL'  (5,420  ) KNK 

420  FORMAT  (110)  

C KNK=NUMBER  UF  SECTIONS  INVESTIGATED 
KNN*0 

410  KNN*KNN  ♦ 1 

IF  ( KNN  .EQ. ( KNK ♦ lY)  STOP 
C RE  AC  SECTION  IDENTIFICATION 

RE  AC  (3.101)  1 1,  12,  13.  INCA  _ 

C READ  PROPERTIES  OF  MAT 
READ  (5,  100)  E,  I, DELTA 

C READ  NUMBER  OF  WHEELS  AND  SOIL  PROPERTIES.  

READ  (3,102)  NN5.WGAMMA.CBR 
C RC AC  WHEEL  LOAD,  CONTACT  AREA, TIRE  PRESSURE. 

RF AC  (5,  100)  WL » CA, TP  

C READ  LOAD  PLACEMENT  AND  MAGNITUDE. 

READ  (5,100)  (SAl(K),SQU(K),S8I(K),K*l,NN5) 

JlK*l  

II*  l 

!« l*0ELTA/12. 

P0IS-0.4 

INTKS*164.+3.*CBR-5.45*WGA‘MA 
ICOVER*  1 

IF  ( J IK  .EQ.l)  GO  TO  69  _ ..  - 

? 3 ICOVER* 1 

IF  (JIK.EQ.2)  GO  TO  70 

69  ICOV  Eft*  INCA  

70  EIM*E*I  I 
NNNK* 1 

501  SU Rft*  I NTKS/ ( FLO AT(  ICOVER )**0.C48455)  . 

GAHMA«FLOAT( ICOVER )**0. 25 
GAMMA«GAMMA+(  13680. -EIM )/EIM 
KS  * SUBR^DELTA/IOOO. 

EQ* ( ( 2.  *(  l . -PO  IS** 2. 1 ) /OEL  TA }*((E*l*K$**3.) /(G AMMA**3. 1 1 **0. 25 
L* ( (2.*E*l*(  1 ,-PO  IS** 2. ))/( EO*DEL TA ) )*►  0.  33333333 

T*( EU*OFLTA*L  )/( B.*GAMMA*< I.+POIS))  _ . 

TT*T/ DELTA 

R2* (l.-POIS  )/( 4 .*GAMM A ) 

S4*  GAMMA  - - 

S2=S4**0.5 

A* ( (S2*R2)/2. )**0,5 

E*( (S2-R2)/2. >**0.5  

AB*  2 . * A*B 
N5*NN5 

00  111  K-I.N5  --  - 

A! ( K)*S  AI(K  ) 

C0(  K )*S  Q0( K ) 

III  BI(K)«S8I(K  ) 

IF  (NN5.EQ.2)  XI  J )«  A I(  1 ) + ( B T(  2)-A  I ( 1 )) /2. 

I F ( NN5  . EO . 1 ) X(  J )*A  1(  1 )♦( B I ( 1 )-A  j < 1 ) I./.2* 

CALI  ULOADlJ) 

IF  (JIK.EO.ll  GO  TO  71 

IF  (NN5.F0.1)  GO  TO  78 

EfC« (£♦!!✓  13680.  )**2»5 


CAC*CA*  *0.5 
WLCMWl/  79.  )**l. 75 

T r>  c = 100  ./TP 

f a i ln- float < icuver  )**.325 

CR1  m=  ( EIC*COR*CAC*T»C  )/<wLC*FAILN) 

c«ire-t=<c;tiTe-t«'*o.9A59^A>/i  ic.c^^c.scspa) 

IE  INNNr.E3.lJ  I COVER  = ICO VFR- 1 
IF  <C-tim.tE.DEFi<  1))  GO  TQ  80 
GO  TO  79 

73  EIC-lc*  11/1368  0. 1**0, 5 
WLC= (WL/ 35.  )**  l . 3 
CAC»CA**0.7? 

CPRC=C0R**0  .9 

F.MLN=F|rOAT  ( ICOVER )**0.3 

C.R  IT  E-t - HLOFAILN  )/(  FIC*CAC*CBRC  ) 

CR  I!  tK  = ( 10.  <*1.  156626)*  ICR  I TER**  0. 559264) 

IF  (NNNK.Eg.il  ICOVER*ICOVEK-l 
IF  (C*l TE3.GF.0CFL( 1 ) ! GO  TO  8C 

79  IF  ( ICOVER. CT .5000)  GO  TO  81 

IF  ( ICOVER. GE. 200)  ICOVER  «ICOVER*20 
IT  ( ICDVER.LT.200)  lC0VER»lCQVER+5 
NNNK  * ? 

GO  TO  501 

80  WRITE  (6,  10) 

10  Ft:  KMC  I//10K,  7HSECTI0N,  10X , 4HLANE , l 0X.4HI  TEN, 10X.8HCCVFRASE  ) 
WRITE  ( 6,  129)  I 1,  12,  13,  INCA 
129  FORMAT  ( 114,  116,  114,  119) 

WRITE  (6,82)  ICOVER 

62  FORMAT  ( / 10X , 25HTH I S SECTION  WILL  FAIL  AT  ,14,11H  COVERAGES.  ) 

GO  TO  410 

81  WRITE  (6,  10) 

WRITE  ( 6,  129)  II,  12,  13,  INCA 
WRITE  (6,83) 

83  FORMAT  (/ 10X,48HPRE01CTEO  FAILURE  IS  SREATER  THAN  5000  COVERAGES  ) 
GO  T 0 4 10 

100  FORMAT  (3F20.5) 

102  FORMAT  ( 1 20 , 2F 20 . 5 ) 

110  FORMAT  ( 2F20.5 ) 

101  FORMAT  (4110) 

120  FORMAT  (4F20.5) 

71  WRITE  (6,72)  ICOVER, UCVH  1) 

7?  FORMAT  ( 110, 5X, 21HFAILURE  DEFLECTION  IS 

JIM2 
GO  TO  73 
END 


,F 10. 5 ) 


OUO  UUUU  O OO  OUO  O O OO  O O O OO  0,0  0 
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APPENCIX  H | 

COMPUTER  PROGRAM  FOR  JSYMM  ETR  IC  LOAO  MODEL  1 


ASYMMETRIC  LOAC  MODEL 

REAL  l 

Cl  MENS  t CN  #Xl(20),AX2(  20),  BXl(  2C  I,  3X  i(  2C  ) , XDEFH  25 1 .ADDEFL  (25) 
COMMOM  X(60),0EFL(fa0)«AI(5),BI(5),QQ(5),N5 
HE AO  SECTICN  IDENTIFICATION 
I l«SECT I ON  NUMBER 
12«  LANE  NJMEER 
I3«  ITEM  NUMBER 
READ  (5.101)  11,12.13 

Rt AO  INFORMATION  ABOUT  REAM  PROPERTIES 
E»  MODULUS  CF  ELASTICITY  OF  BEAM  - MPS/SQ.IN. 
l»  MOMENT  CF  INERTIA  Of  BEAM  PER  F03T  OF  WIDTH  - IN.**4 
RFAO  (5,100)  E , I 

NLOAO*  NUMBER  CF  WHEELS. 

HHEFlS>WHEFL  SPACING  IN  INCHES. 

WIDL AN»W  1 OT  CF  THE  TRAFFIC  LANE  IN  INCHES. 

READ  (5,2)  NLO^O, WHEELS, WICLAN  • 

NSHI FT«  NUMBER  CF  POSITIONS  OF  THE  ASSEMBLY  BEFORE  IT  IS 
MOVE C LATERALLY  THE  WlCTHOF  THE  ASSEMBLY.  THIS  VALLE  IS  L'SFO 
ONLY  FOR  MULTIPLE-WHEEL  TESTS.  FOR  THE  LOADING  SECLENCE  ShUwN 
IN  FIGURE  24  CF  THIS  REPORT  - NSHFT«3.0.  THIS  VALLE  FOR 
SINGlE-wHEEL  TESTS  IS  SFT  EQUAL  TO  UNITY. 

CONT  A.CCNT  ACT  AREA  OF  TIRE  IN  SQ  . IN  . 

©•TOTAL  LOAD  ON  ONE  WHEEL  IN  KIPS. 

RFAO  (5.2)  NSHIFT  ,CONT  A,  Q 

NFOS 1 T* NUMBER  CF  POSITIONS  OF  WHEEL  ASSEMBLAGES  TO  COVER  HE 
TRAFFIC  LANE  WICTH. 

READ  (5,421)  NPOSIT 

TIREW-WIDIH  OF  THE  TIRE  PRINT  IN  INCHES. 

TIREl«LENGTH  CF  THE  TIRE  PRINT  IN  INCHES. 

ABOVE  TWO  VALUES  BASEC  UPON  A RECTANGULAR  TIRE  PRINT. 

IF  (NLCAC.EO.l)  TIREW-WICLAN/(FL0AT(NP3SIT)) 

IF  (NL0AC.EQ.2)  T I REW»W ICL AN/ ( 2 ,*FLO AT( NPD S I T ) ) 

TI  REL*C  ONT  A/T 1 REW 
I*  I ATI  REt/ 1 2. 
fclK«0 
Q*Q/TIREW 

AT-LENGTH  CF  TIME  THAT  TFE  LOAC  IS  APPLIEO. 

BT-LENGTH  CF  TIME  THAT  THE  LOAC  IS  OFF  BEFORE  THE  LOAD  IS 


i 


1 


I 


C Hulf.  A GAIN. 

*•*  at  (r>(nj)  «i,ei 

■]  c « m:  t t * v 

j C 4T  AND  RT  MUST  BE  SO  SELkCTEC  Thai 

i C (41  « BT  MNRPS  IT  =1 

« C THt  AlCVF  PeCUTREMENT  IS  NECESSARY  FOR  l CUVER4GE  10  BE  T.iE 

i c I'iic  m*'  iMr. 

C REA'J  l.VFCRMAT  ICN  Cf  V6RI4TIUN  OF  *31) ULUS  OF  SOIL 

C i'  IMS  UCCF'H  CMY  THE  RESIDUAL  DEFORMATIONS  ARE  OF 

'j  C INTEREST,  THIS  THE  VALUE  OF  C CAN  BE  ARBITRARY. 

C IN  Ti-li  FRCCFAH  AC=l.O  ( SEE  REPORT- SECTION  ON  ASVMMFTrIC 

i C LOAD  MCDLL). 

RT An  (5  ,1)  c . C I ,EN,AC 

C XC*niSTANCF  FROM  CRICIN  TO  ECGE  OF  TRAFFIC  LANE, IN  T-»  I S 

C FkCC.MF  THF  CPICIN  FAS  BEEN  SELECTri)  AS  A POINT  1C  FEET  LEFT 

C OF  TFF  ECCF  fF  TFE  TRAFFIC  LANE. 

>0*1/0. 

IF  (MCAC.EC.1  ) CL  TC  650 

fC  651  J- I , NPCS  IT 

IF  (J.GT .NSFIFT ) CO  TC  720 

A X l < J ) = XO  + FlCAT  (J-IMT  IREW 

AX?  (J)*  AXl  ( J ) ♦ Ft C AT  ( NS  H I FT  ) *T  ItR  E W 

PX1 ( J ) * *X 1 ( J ) *T JRFW 

9X?  (J)«  AX2I Jl  + T 1«£W 

GO  TO  651 

7?0  AXl  (J)» XG*2.*FL0AT (NSHIFT  ) *T IREW hFIO A II J- 1-NSHIF 1 )* TIRE  V 
AX?  (J)* AXl (JI  + FLCATI NSHIFT  )*TIREW 
Pxl  ( JJ*  AXl ( U)hT IREW 
fix? (J)» AX2( J ) ♦T  IREW 

651  CONTINUE 
C>  C I C (j  > 2 

65 C OC  653  J* l, NPCS IT 

AXl  (J)*  X0+FL0AT1  J-1W  IREW 
653  I XI  ( J)* AXl ( J)*T  I H Em 

652  N5  * M C A C 

OC  65A  XK*l,N5 
..  659  CC  ( K x ) * C 

C ESTABLISH.  OISTANLES  TO  WHERE  CFFL  ECT  ION  S ARE  TO  BE  COMPLIED. 

C LEI  THIS  ARRAY  eF  XU). 

C COMPUTED  DEFLECT  ICN  AT  ONE-TENTH  POINTS  ACROSS  TRAFFIC  LANE. 

00  303  „J*I  ,11 

3 03  X(J  J)=XO»FLCAT(.JJ-l)*rt  ICLAN/1C. 

WRITE  (6,10) 

WRITE  (6,129)  11,12,13 
C NOMCOV*  NUMBER  OF  COVEPACES. 

C NFASS*  NUMBER  CF  SPECIFIC  PASS  WIIFM  A COVERAGE  AT  THE  END  OF 

C WHICH  THE  CEFLFCTICN  IS  CESlREC. 

C NFASS  MUST  BE  LESS  THAN  OR  ECUAL  TO  NP3SIT. 

C NCFASS*  TOTAL  NUMBER  OF  PASSES  NECESSARY  TO  COMPLETE  THE  GIVEN 

C NUMBED  CF  CCV  ERAGES  . 

810  RF A C (S.A20)  NJMCOV  , NPASS 
NCPASS*  M'MC  CV  HNFCS  IT 

C LUP*  COMPLET  t NUMBER  OF  PASSES  FOR  ALL  FULL  AND  PARTIAL  COVER. 

IUF= NO PASS* NPASS 


c 

c 


IS  DESIREO 

BE  REPRE  JENTED 
INITIALLY 


BY 


LI,cfn!rMN<'L  T,M£  TU  hHRE  ELECTION 
FU  M£*F  LQAT  (LUP)  *AT  ♦ FLOAT  (LUP  )*ET 

LET  THE  TCTAL  CEFLECT ICN  AT  EAO  POINT 
XCEFKJ),  ZERC  ALL  OF  TFESE  LOCATIONS 
00  AAA  J*l,  U 
AAA  XQEFl ( J )*0, 

^COMPUTE  THE  CEFLECT IONS  FOR  EACF  LOAD  APPLICA 1ICN  AND  RELEASE. 

m*NFOSIT 
N(jMCC2»  NUMC0V42 
DO  30V  NN«2  ,NUNC02,? 

IF  1NN.CC.NUMCC2)  M3»NPASS 
IF  (M3.EC.0)  GC  TO  301 
DC  302  **»!,*} 

IF  INLCAO.EC> 1 ) CO  TC  660 
AI  (D-ATUMH) 

A I C2)« AX2(HM) 

BI  <1)«3X1(NM) 

61  (?)*BX2(MM) 

KHM-O 
GC  TO  661 
AI ( 1 )• AX1 ( MM) 

81  ( 1 )•  BX1  ( MM) 

XHM«0 

NMl  «MN*NFCS IT -I 2*NP0S IT-2*MM«2-KMM ) 

T ST  APT»FLOAT(  NMII<AT4FL0AT<NHD<BT 
TlME-FTlME-TSTART 
Tl»AT 

C0Mf-04  01*(Tl*<  l./AC-l.l4TIMn«EN-IC40  l*  t TIME-  T 1 )**£  N) 

OXS»TI REl/COMP 
WLAMCAa  <eXS/(A.*F*in**0.25 
CO  309  J»l,U 
CALL  UULOAO  < J.WLAMCA, BKS ) 

AOqEFLI  JI-OEFKJI 
00  663  J«1,U 

X0EFKJI*X0EFLIJ)4ACCEFL<JI 

I !■  I !♦  1 


660 

661 


309 

663 


KMM«KMM4l 

IF  (KMM.EO.l)  GC  TO  661 
302  CONTINUE 
301  CONTINUE 
K 1 1 » II - 1 
JNN«NN-? 

WRITE  (6,130)  JNN 
WRITE  (6,131)  MM.K1I 
WRITE  (6,132)  Cl , EN 
WRITE  (6,133) 

WRITE  (6,13A) 

WRITE  (6,135)  ( XOEFL(J),X(J),J-l,D  ) 

GO  TO  8 10 

1 FORMAT  (AF10.5) 

2 FORMAT  U5.2F10.5) 

10  FORMAT  (//  lOX.  7HSECT  ICN,  10X,  AH.ANE,  1C  <• AhlTEM 
100  FORMAT  (2F10.5 ) 


101  FORMAT  (3110) 

120  CORMAT  (114,116,114)  „ . 

130  . OR  HAT  (6X,4lHTHE  FOLLOWING  RESULTS  OCCUR  AT  THE  ENO  OF  / 

* 17X,I4,IX,9HCC*E«ACES) 

131  FORMAT  (6X,33HTME  LAST  POSITION  OF  TFE  LOAD  WAS  , 1 5 , 1 X,  3H  AND/ 
*2X,11HA  TOTAL  Cf  , 15,  IX,  38H.0ACS  FAVF  BEEN  APPLIED  TO  THE  SUTEM. ) 

132  FORMAT  USX,  19HTHE  VALUE  OF  Cl  IS  * F6 • 3/ I5X, 1SH  THE  VALLE  CF  EN  IS  , 

* F6.3) 

133  FORMAT  (SX, 37HTHI  CEFORMAT ICN  PATTERN  IS  AS  FOLLOWS) 

134  FORMAT  (4X,20HRES10UAL  CEFORMAT ION  , 4X, 2CHD I STANCE  FROM  ORIGIN/ 

I!  • 10X.6HIACHES  , 1NX,6HNCHES  ) 

i 135  FORMAT  I 10 X, F8 .5, 15X , F 10 .5 > 

'i  420  FORMAT  (2110) 

!i  421  FORMAT  (110) 

i END 


fi  , 

E-0-»J 


o o 0,0 o o o1  o j . i i no;  noi  o | jo  ooooo 


APPENDIX  I 


C CHFUT E R FR03PAW  FOR  NOP ENT  TRANSFER  INVESTIGATION 


MOMENT  TRANSFER  ANALYSIS 


VLASOV  AFC  LEONTEV  MODEL 
HYFEneOLIC  STRESS  OISTRIEUTION 
REAL  ( » KS  1 1 FTKS « KUEFl, MOEFLLf MOMES*MQM£»L * II 

01  MENS  I CN  $AI15),SCC(5  ),  SB  II  5 >♦  XS1  2C  ),DSI  20  tXQEFU  20  *MO£Fl  (2  0) 
I FDEFLIifO  ) »KEFLL(  20 >•  CCE FL ( SCI 
COMMON  X(60),A|f5>,CC(  5),  B l(  5 ),l » A,  e , AB, S2» $6,R2,F  , I ,N5,DEFl(60) 
I SLOPE (60) » HOME (60 )» SHEAR! 60  )»CM(5)»01ST(5)*N2*AA*I 

FNK» NUMBER  CF  SECTIONS  INVESTIGATED 
REAC  (5,620)  KNK 
RNN-0 

610  RNN-RNN+l 

IF  (XNN.EQ. (KNK«1) ) STOP 

READ  I NF CRM AT  ION  ABOUT  BEAM  PROPERTIES 
RE  AC  (5,100)  E,  I, CELT  A 

CELT  A* WIDTH  OF  BEAM  IN  INCHES. 

READ  INFORMATION  ABOUT  FOUNDATION  PROPERTIES. 

PF  AC  (5  ,102)  KKK,NN5,P0IS 

FKK» NUMBER  CFSETS  OF  CATA  PER  SECTION. 

NN5« NUMBER  CF  STATIC  UNIFORM  LOACS. 

READ  (5,100)  (S  A!  ( K ) , S C0(  K )» $6  t(K  ),K«l,NN5) 


1 1*  I 

l« 1 *UfcLT A/ 12. 

EIM«E*I 1 
NN»0 

602  NN*NN»1 

IF  INN. C 0. ( KKK+! ) ) GO  TO  610 
REAO  (5,6)  N2.N0P 

N2*  NUMBE  F CF  MOMENT  RELEASES. 

REAO  (5,101)  tl, 12, 13,  16,  ICOVER 

1 NTKS*SU EGR  *( E MOCULUS  IN  L eS  ./CU  . IN  . 


REAO  15  ,999)  INTKS 

RCASAI«c1s?1nCES TO^eViNFi’iNG^^STATIC  UNIFORM  LOAD  IN  IFCHFJ. 
SB1*CISTANCE  TO  ENC  OF  STATIC  UNIFORM  LOAD  IN  INCHES. 

1 1* SECT  t CN  NUMBER 


I2-LAFE  NJM8ER 
13*  I TEN  NUMBER 
16* LOC AT  JCN  NUMBER 


, 


, 


I 


r»nj 


r- 


C !COVfR« NUMBER  OF  COV ERAGES 

C ■ LET  l*CENTER  OF  PANEL  2»QUARTER  PONT  3F  PANEL  3*JQINT 

WRITE  (6,3) 

WRITE  (6,129)11,12, 13, U,  ICOVER 
. SUBP*!N7K$ 

IF  (ICOVER.GT.O)  GO  TO  SOI 
GAMMA*  1 ,0 

GAMMA* GAMMA* (13660. -EIM)/EtM 
GO  TO  502 

501  GAMMA* FLOAT (ICCVER)*»0.25 
GIMMAaG  AMMA*( 13680. -EIMJ/EIM 

502  RS*SU8R*CELTA/1000. 

A03  feO*  ( (2. ♦(l.-PCIS**2.) )/C£LTA)*(  ( E*I*K$**3. )/( GAMMA**  3. ) )**0.25 
l* ( (2. *E* I * ( 1. -FC IS  **2 • I )/( EQ»CELTA))**C. 33332331 
T«(E0*DELTA*L)/(6.*GAMMA*(  1.4P0IS)) 

TT«T/OELTA 

R2*  (l.-FOIS )/ (A.AGAMMA ) 

. SA*GAMMA 
S2*  S4**0«5 
AA-SORT  (KS/ (2.*T)| 

A*  ALPHA 
8-BETA 

A*(  (S2*P2)/2.)M0.5 
B*(  (S2-R2)/2.)*«0.5 
AB-2.«A«6 
M5-NN5 

CO  111  Ml, NS 
Al  Ml-SAHK) 
o)*SCO(K) 

in  ei « Ki-sei (k » 

00  1 14  J*  1 , NOP 
X( J)*XS ( J ) 

CALL  ULOAO  (Jl 
>OEFL(J)*CEFL(J) 

1 I A CCNT1NUF 

C CALCULATE  THE  MOMENT  AT  TFE  JOINT  LOCATION  ASSUMING  THE 

C BEAM  TC  eE  CONTINUOUS. 

J*1 

X(J)*AI  (l)*(BI(2)-AI(l))/2, 

CALL  ULOAO(J) 

ACT  MON* HOME ( J) 

ERRS* 10000. 

JI  J-l 
N2*l 

OIST ( 1 ) * X ( J ) 

CM( l )■-  ( ACTMONAO.Ol ) 

117  CO  115  J* 1 , NCP  ' 

X(J)*XS(J) 

CALL  SEMIM(j) 

MOEFl ( J )*  CE  F , ( J ) 

115  CONTINUE 

CO  116  J* 1 , NOP 
X( J)*0IST(1 )♦( CIST ( 1J-XSI J ) I 
CALL  SEMIM(J) 


moeflu  j)*reruj) 

C’.TEF  1 ( J )»  MCEFL  ( J )♦  MCE  FLU  J I 
116  FCCFLU)«XCEFLCJ1*CCEFUJ  ) 
FERC«  CMU)/ACrMON*lOO. 
£RRO«*0 .0 


UO  ?00  J*  l * NOP 

POO  i'R^03*EfF0«  + «eS(FDEFLl  J l-CSU  ))*•?. 
IF  (JlJ.EQ.2)  GC  TO  401 
IF  (ERRCR,GT.E«RS>  CO  TO  400 

errs*er»cr 

CM(1)*CM(1)-  (ACTMCN*0.C1) 

GC  TC  117 

<,00  CHIlKKdH  I ACTMCNM0.C1I 


Jl  J*2 
GO  TO  117 

401  WRITE  (6,130)  ACTMQN 
WRITE  (6,131)  Cell) 

WRITE  (6,132)  PERC 
WRITE  (6,133)  ERROR 
WRITE  (6,134) 

WRITE  (6,135) 

WRITE  (6,13b) 

WRITE  16,137) 

WRITE  (6*,  113)  (XS(J),CS(J),XCEfLU),WC£FUJ),MOEFLL(JI,COEK(JI  » 

I FDEFl  ( <J) , J*  1 , NOP ) 

3 FORMAT (// 10 X»7HSECT  ION,  10X,  4KANE,  1CX,4HTEM,  ICX.8HLCCATICN, 

1 1DX.9HCCVERAGES  ) 

4 FORMAT  (2113) 
ino  FORMAT  (3F20.5) 

101  FORMAT  (5110) 

102  FORMAT  (2115,  F15.5) 

110  FORMAT  (PF20.5) 

113  FOPMAT  (7F15.5) 

130  FOR  MAT  < / / 5 X *65FTH  C^MCK ENT  ^ AT  THE  LOCATION  OF  THE  HINGE  IN  THE  CONT 

i3i*}o(iXn/55^e£m^  concentrated  moment  applied  to 

♦THE  END  OF  THE  SEMI— INFINITE  BEAM  IS  ,F 5 , 2 , IX, 8HKI P* IN »•  I 
132  FORMAT (/5X.78HTFE  RATIO  OF  THE  APPLIED  CONCENTRATED  MOMENT  TC  THE 

♦ MOMENT  AT  THE  HINGE  IN  THE  /5X,  18FC0NTTNU3US  3EAM  IS  ,F?.2,l*» 

♦8HPEPCEM.  ) ,,  ,, 

.,1  FORMAT (5X.23HT HE  SIMULATION  ERROR  IS  «F9.5) 

. FORMAT (/3X. 13HCISTANCE  FRCM  , 6X,  fcHACTUAL  , 7X,  I CHOEFLEC  TICN  ,3X, 

• JSSSeFLCCT  CN  .5X.I0HCEFLICT1C.N  , 5X,  1CHDEFLECT  ION  ,7X,5H  TOTAL  ) 

l„  ;7,.10«,IUECT». I .M.IIHOF  l*IM* 

♦rh OF  SEMI-  ,7X , 8H0F  SEMI-  ,fcX,9hOF  HINGED  ,eX«5HM00EL  ) 

U6  JoSJuTl5!x,4HBEAM  ,6X  , 13F  INF  IN  IT  E BEAM  ,2X,  13HINF1NI  TE  BEAM  ,6X, 
♦4HBEAM  ,9X  , 10HCEFLECT  ION  ) 

137  FOR  MAI  (51X,7H(RlGHT  ) , 9X.6HILEFT  ) ) 

136  FORMAT (7X.6HINCHES  ,6( 9X,  6HINCHES )// ) 

420  FORMAT  ( 1 10 ) 

999  FORMAT  (F10.0) 

END 


E-O-J 


APPENDIX  J 


subroutines  for  compute^  programs 


SUBROUT INE  ULOAD( J ) 

REAL  ItMOMES,MSMAl,HOME,l 

COMMON  X{60),An5),Q0l5J,BI(SifLtA,B,ABtS2»S4»R2*EtI  »NS,OEFUn>, 
l SLOPE!  60). MORE!  60  It  SHE AR  ( 6C  ) *CM  ( 5)  ,01  ST ( 5 » ,N?  »AA , T 

OEFL(J)«0.0  _ 

SLOPE! J )*0.0 
HOME! J ) *0 .0 
SHEAR!  J)*0.0 
CO  600  K» 1, N5 

IF  <X(  J J.GT.AKK  ))  GO  TO  601 
AIETA*! Al(K )-X(  J ) )/L 
8LETA»(  BIU  )-X(  J ))/L 
Cl*  U**A.)/!  A.*A*B*S4*E*l)*R2 
C2*  ! L ♦♦<»«)/!  A** A*B*SA*E*  l )*AB 
GO  TO  603 

601  IF  (X! J ).GT  • B I ( K ) ) GO  TO  602 

ALET  A* ( X ( J 1-A  I! K ))/L 

BL  E T A*  ( B 1 1 K )-X  ( J ) ) /L 
Cl*-!L**A.)/! A.*A6B*SA*E*l  )*R2 

C2*-(L»»A.)/! A.«A«B»SA»£»1)»AB . 

Cn  yc 

602  ALETA*|X!J)-AI!K  )J/L 

BLET  A* ! X ( J )-Bl<  K 1 >/L 

Cl-  !L**A.)/<  A.*A*B*$A*E*I)*R2 
C2«  (L**A.)/!A.*A*B*SA*E*I)*AB 

603  AAX»A*ALETA 

BAX»B*ALET A 

ABX-A*BLETA 

BBX- B*BLETA 

FA1*SIN!BAX  )/ ( EXP ( A AX  ) ) 

FA2«C0S!BAX I/!EXP!AAX1» 

FA1I«B»FA?-A»FA1 ' 

FA2  l»*(  A*FA2+B*FAl) 

FAl?«R2*FAi-AB*FA2 
FA22»R2*F A2*AB*FAl 

FA13»8*!3.*A**2.-B**2.  )*FA2»A*! 3.*B**?.-A**2. )*FAl 
FA23*  B*( B*-2.-3.*A**2. )*FA l*A*< 3.*B**2.-A**2. ) *F A2 
F81*SIN!EBX)/(EXP!ABX>)  • 

FB2*C0S ( EBX J/T EXP! ABX  I I 
F81 l»B*FB2-A*FB  1 

FB 2 1»-(  A *FB2+B*F B_IJ 


o o 


V 


FR12*R2*F81«AB*FP.2 

F«22=«2*FB24AB*FB1 

Ff>  l 3*3*1  3.*A**2.-ft**2,  )*FB2*A*(  3. *3**2. -A**2.  ) *FB1 

FB?3=0*( fi**2.-3.*A**2.  )*FB t*A*( 3. *3**2. -A* *2, ) *FB’ 

IP  (X(J } - GT  . A I ( K ) ) GO  TO  fcOA 

VSMAL=Cl*(  FAI-FB1)+C2*(  FA2-FB2I 

Lf>  M AL=  C 1 * ( F B 1 l-F A 1 1 ) +C  2*(  FB  21-FA 2 1 3 - 

PSMAt  «Cl*l  FA12-FB12)*C?*(FA22-F8  22) 

CS  F ALS  C 1*( FB 13-F A13 )+C2*( FB23-FA23) 
osmal*-o$mal 

CS^Af-CS^AL 
GC  TC  606 

6GA  IF(X(J).GT.eiUn  GO  TO  6C5 
CC«U**A.)/(SA*E*I) 

VSMAL*C1*( FAl*F8ll*C2*(FA2*FB2)*CC 
CSHAL*C1*(  FAU-FBllltC  2*1  FA21-FB21) 

FSFAL«Cl*( FAl?*FB12)*C2*(FA22+FB22l 
G$MAt*Cl*<  FA13-FB13)*C  2* ( FA 23-FB  23 1 
CSMAL--OSMAL 

csmal«-osmal 

GO  TO  606 

605  VSMAL«C1*(  FB1-FA  l ) *C2*(  FB2- FA 2 ) 

OSFAL*-(Cl*(  FBIl-FAlU*C2*(FB2l-FA2in 

KSM  AL*C  1*(  FB  12-FA  12 ) *C 2*1 FB 22-FA22)  ' " 

qS/*Al»-(Cl*(F813-FAI3)*C2*(FB23-FA23)) 

606  DEFLS»QOU)*VSMAL 

S LOPES* QOIK )/L*OSMAl 

MOMES*-( QO(K)*F*l)/(L*L )*MSMAL 

SHEARS»-(QO(K )*E*I )/(L**3.)*QiMAL_ 

CFFL( J j *DEFl( J ) ♦OEFt S 
SLOPEIJ  )*SLQPE( J ) ♦SLOPES 

M0ME1J )»MOME( J )+MOMES  _ 

600  SHEARU  )*SHEAR(  J )+SHEARS  ' “ 

RETURN 

6N0_  __ 


SUBROUTINE  WULOAO  ( J t WL AMOA, BK S ) 
'"YNFINITE  3EAM  ON  WINKLER  HCOEL 
UNIFORM  LO AOS 

CO FMCN  X ( 60  ) , 0 EFU  60  I , A I ( 5 ) ,B 1 ( 5 1 ,00 (6) ,N5 

CEFLI J ) *0.0 
DO  800  K»l,N5 

IF_«XIJ_».GT.AI(K)J  GO  TO  801 

AW*  At(K)-X(J) 
ew«  BI(K)-XU) 

GO  TO  803 

801  IF  (xm.GT.bllK))  GO  TO  802 
AW*  X(J)-AHK) 

BW*  BHK)-X(J> 

" GO  TO  003 

802  AW*X(J)-Al(K) 
eW*X(J)-BI(K) 

803  ALAMCA*  AW*WlAMOA‘ 


BLAPOAs  BW*WLAPCA 

CIA*  COS  ( AL  AMOA )/(  EXP  ( A|.  AMOA  ) ) 

CLB*  COS(BLAMOA)/(EXP(BLAMDA)) 

IF  (X(  J ) .gt  . A I ( K ) ) GO  TO  eOA 
CEFLS={ cn(K  )♦( 0lA-0t8  ) )/{  2.»BKS) 

GC  TO  800 

804  IF  (X(j|.r,r  .3IIK))  GO  TO  805 

CFFLS*!  QO(K  )*(  2.-DLA-0LB)  )/(  2.*SKS) 
GO  TO  800 

80  5 CEFLS*-!  00 ( K )*(  DLA-OLBH/I  2.*BKS) 

800  CEIHJ  j*CEFL!J  )»OEFLS 

RET  JRN 
END 


SUC ROUT INE  SEM 1M( J 1 
REAL  I»  POMES, M3ME,  l 

CCM MON  X I 60 ) , A 1 < 5 ) , QO ( 5) , 8 1 ( 5 ) ,L ,A  «B  ,AB  » S2  » S4 , R2 ,E  , l , N5»  DE  F L (5  0 ) 
1 SLOPE! 60), MONEI 60 ), SHE AR I 60) ,CM (5) ,DI ST  15) ,N2 ,AA ,T 

DEFLI J )«0.0  _ 

SLOPE!  j )*0 .0  ~ ~ " 

MOME! J )»0.0 
SHEAR! J )*0.0 
00  300  K*  1,N2 

D*-B*(S  2*! 3.*A**2.-B**2.)*( A.*A*AA* T*L**3. ) /!E*l ) ) 

Cl»-!!CM!K)*L**2.  )/!E*l*D))*«A*S2+(2.*\A*T*L**3,  )/!E*I)) 

C2* (CM!  K )*S 2*e*L**2. 1/1 E+  1*0) 

IF  (X!  J J.GE.01STU  ) ) GO  TO  301 

X0IST«0  IS T < K ) -X (J J 

AANEG«AA*xblST 
CEFLS«C2/( EXP!  AANEG)) 

JCEFL(J)*DEFL(J  )*OEFLS__ 

GO  TO  300 

301  xniST«X(J)-CIST(K ) 

AX*  A*X 0 IST/L 
BX*B*XD  IST/L 
F1*S IN! OX )/ ( EXP! AX ) ) 

F2*C0SJ  BX  )/  ( EXP!  AX  ) ) 

FI  1*  B*F  2-A*F l ' 

F21«-( A*F2*B*F 1 ) 

_ F12*R2*f  i-AB*F2 

F22*R2*F2+A6*Fi  ~ 

F13*B*( 3 • * A * * 2 .-B**  2 • )*F2*A*( 3,*B**2.-A**2. )*F l 

F23* B*«  B**2.-3.*A**2.i*Fl+A *J  3.*6**2.-A**2. )*F2 

QEFLS*C l*Fl*C2*F2 
SL0PES«Cl*FU+C2*F21 

PCMES*CI*F1  2+C2*  F22 

SHEAR$*Cl*F 13*C2*F23 
DEFL<J)*OEFK  J HOEFLS 

SLOPF!  J )=S10PE( J )+SLOPE_S_ 

MOME(J)*MUHE(J  MMOMES 
SHEAR! J )«SHEAR(J )+SHEARS 
300  CONTINUE 
RETURN 
ENO 


